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Einleitung und Debersictit über die Hilfsmittel der modernen 
Reclieflteclinik, insbesondere die mattiematisctien Tabellen. 



§ 1. Die Entwicklung der Rechentechnik kann auf eine lange 
Geschichte zurückblicken ; denn die Wahrnehmung der Tatsache, 
dasa öfters Dinge gleicher Art und in grosser Zahl — Menschen, 
Tiere, Pflanzen u. s. w. — vorhanden waren, musste schon früh 
den Wunsch erwecken, die Anzahl derselben zu bestimmen und 
derartig anzugeben, dass auch anwesende Genossen des schauen- 
den Forschers hiervon Kunde erhalten konnten. Eine solche Mit- 
teilung konnte anfänglich nur in einfachster Weise geschehen, 
z. B. durch Emporheben von Fingern, durch Nebeneinanderlegen 
von Steinchen, Körnern, Früchten, in spätem Zeiten durch eine 
Reihe von Strichen, Punkten und dergleichen auf Steinplatten, 
Tontafeln, Papyrus u. s. w. 

Mit der Entwicklung der Sprache stellten sich dann später 
die Zahlwörter ein, d. h. Ausdrücke, durch welche die Menge der 
vorhandenen gleichartigen Einheiten zu einem Begriff zusammen- 
gefasst und demgemäss in kürzester Weise bezeichnet wurde. 
Indessen hatte die UnmögUchkeit, alle Zahlen bis zu einer einiger- 
massen beträchtlichen Höhe durch neue Zahlwörter zu bezeichnen, 
die Bildung systematischer Anordnungen zur Folge, welche die 
Benennung aller Zahlen durch geeignete Verknüpfung weniger 
Wörter ermöglichten. Am verbreitet sten war damals schon das 
dekadische Zahlensystem, das je 10 Einheiten einer Stufe zu einer 
Einheit, der nächst höheren Stufe zusammenfasste. Seine Ent- 
stehung verdankt es der von altersher gebräuchlichen Abzahlung 
an den Fingern. Ebenso entstand das in vielen Sprachen hervor- 
tretende, wenn auch niemals konsequent durchgeführte quinäre 
System durch Abzählen an den Fingern einer Hand und das bei 
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den Azteken und Kelten gebräuchliche vlgesimale System durch 
Abzählen an Fingern und Zehen. 

Wenn die Menschen bei Aufstellung dieser Systeme mehr 
unbewusst der Anleitung der Natur folgten, so waltete dagegen 
der überlegende Geist bei der Schöpfung des dodekadischen und 
sexagesimalen Zahlensystems, die — entsprechend den praktischen 
Bedürfnissen — bequeme Teilungen gestatteten. Ihre innere Be- 
rechtigung wird durch den zähen Widerstand bezeugt, den sie 
dem Vordringen des Zehnersystems entgegenzusetzen vermochten. 
Insbesondere behauptet das babylonische Sechzigersystem in der 
Zeit-, Kreis- und Winkelteilung noch heute seine mindestens 
4000-jährige Herrschaft. 

Das Zusammenfassen mehrerer Einheiten einer Stufe zu einer 
Einheit der nächst höheren Stufe veranlasste die Unterscheidung 
zwischen Zahlzeichen für die Stufenzahlen und solchen fQr die 
Einerzahlen. Dabei hatten nach dem multiplikatorischen Prinzip 
die Einerzahlen die Rolle von Koeffizienten — 583 wurde sche- 
matisch geschrieben: 5 CSX 31 — . Es ist bei einer derartigen 
Zahlenschreibung naheliegend, die Zeichen für die Stufenzahlen 
ganz wegzulassen. Solche Versuche sind denn auch wirklich 
gemacht worden, indessen scheiterte die allgemeine Durchführung 
dieser Abkürzung daran, dass man Zahlen wie 2C3 (203) nicht 
zu bezeichnen wusste. Erst die Erfindung des Nullzeichens durch 
die Brahmanen Vorderindiens im 4. Jahrhundert n. Chr., das eine 
der grössten Wohltaten darstellt, die jemals ein Erfinder der 
Menschheit erwiesen hat, gestattete eine konsequente Durch- 
führung der vorgenannten Abkürzung. Erst mit dem Gedanken^ 
das Fehlen von Einheiten einer Stufe durch ein Zeichen, — die 
„0" — ersichtlich zu machen, war das Positionssystem erfunden 
und damit ein System geschaffen, welches in wissenschaftlicher 
Hinsicht das Beispiel — wohl das einzige — eines absolut voll- 
kommenen Zahlensystems ist. Was die praktische Bedeutung 
betrifft, so genügt der Hinweis, dass es gestattet, mit denselben 

zehn Zeichen 0, 1, 2 9 jede Zahl, so gross sie auch sein mag, 

zu bezeichnen, — daher es auch Gemeingut des Volksunter- 
richts aller gebildeten Völker geworden ist. Die logische Weiter- 
entwicklung des Stellenprinzips brachte dann am Ende des 
16. Jahrhunderts eine für das praktische Rechnen höchst wichtige 



^Vereinfachung: das Rechnen mit Dezimalbrüchen, Besonders ver- 
dient um die Einführung derselben machten sich die Mathematiker 
Viete, Stevin und Bürgi. Von diesen war es namentUch Stevin, 
welcher in seiner Abhandlung „La Disme" (1585) die Nützlichkeit 
der Dezimalbrüche als Hilfsmittel der Rechentechnik hervorhob 
und klar den Gedanken aussprach, alle Rechnungen des Geschäfts- 
lebens ohne Brüche, nur mittels ganzer Zahlen auszuführen. Von 
ihm unabhängig gelangte der Schweizer Mathematiker Jost Bürgi,* 
der sich zuerst eines Punktes zur Abgrenzung von Dezimalstellen 
bediente, zu der gleichen Vereinfachung. Er lehrte auch schon 
die abgekürzte Multiplikation. 

§ 2. Neben der besprochenen Entwicklung, die mit der Ein- 
führung der Dezimalbrüche ihren Abschluss findet, setzt schon 
sehr früh und unabhängig davon die Einführung eines zweiten 
Hilfsmittels der Rechentechnik ein, nämhch die Darstellung des 
numerischen Verlaufs einer Funktion durch Tafeln oder die 
yfTabuUerung^^ . 

Der erste, der sich mit diesem Zweig des praktischen Rech- 
nens beschäftigte, war der grosse Hipparch, der zwischen 161 
und 146 vor unserer Zeitrechnung astronomische Beobachtungen 
anstellte. In ihm erkennen wir den eigentlichen Erfinder der 
sphärischen Astronomie und der Sehnenrechnung, für welch 
letztere er auch die sogenannte Sehnentafel berechnete. Nicht 
nur bezeugt dies Theon von Alexandrien in seinem Kommentare 
zum Almagest, sondern es geht auch aus der einzigen von 
Hipparch selbst erhaltenen Schrift sicher hervor, dass er bereits 
mit Hilfe geometrisch abgeleiteter Rechnungsregeln unter An- 
wendung dieser Sehnentafel manche Aufgabe der sogen, sphäri- 
schen Astronomie zu lösen wusste. Indessen ist sowohl die 
Tafel selbst, als auch eine ihre Berechnung behandelnde Schrift 
von Menelaus spurlos verschwunden, so dass man heute auch 
nicht das leiseste über die Einrichtung jener Tafel auszusagen 
im Stande ist. Dagegen übten die Arbeiten dieser beiden hervor- 
ragenden Mathematiker grossen Einfluss aus auf Claudius Pto- 
lomäiis (140 n. Chr.), der das vollendete, was Hipparch und Me- 



* Jost Bürgi oder Justus Byrgius stammt aus dem Toggenburg und lebte in 
der Zeit von 1552 — 1632. ursprünglich Uhrmacher und Mechaniker ; später Astro- 
nom unter Kepler. 
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nelaus begonnen hatten. Er schuf für den astronomischen Gebrauch 
eine Trigonometrie von so vollendeter Form, dasM sie weit über 
ein Jahrtausend nicht überboten wurde. Sie findet sich vereinigt 
mit seinem astronomischen Lehrgebäude in den 13 Büchern des 
berühmten Abnagest.* Im 1. Kapitel des 1. Buches finden wir 
die griechische Goniometrie mit der Sehnentafel, welch letztere 
die in Teilen des Radius ausgedrückten Zahlenwerte der Sehnen 
für die Winkel von bis 180 <> (von 30 zu 30') wiedergibt. 
Aehnliche Tabellen, wie sie Ptolomäus für die Sehnen geschaffen, 
finden wir etwas später auch bei den Indiern und Arabern, im- 
merhin mit dem Unterschiede, dass sich diese beiden Völker bei 
ihren Rechnungen nicht der Sehnen-, sondern der Äwms- Tafeln 
bedienten, indem sie an Stelle der Sehnen die Halbsehnen als 
Funktionen des halben Winkels einführten. Unter den Arabern 
war es besonders Äbul-Wefa (930—998), welcher sich als Rechner 
hervortat. Er erfand unter anderem eine neue Methode zur Be- 
rechnung von Sinustafeln, die den Sinus von 1/2 ^ mit einer Ge- 
nauigkeit lieferte, wie sie eine auf 9 Dezimalstellen gerechnete 
Sinustafel aufweisen würde. Ferner führte er die trigonometrische 
Funktion Tangente ein, indem er die Werte des Quotienten 
sina : cosa berechnete und in einer Tafel vereinigte. Gleichzeitig 
entstanden auch die Hakemitischen Sinustafeln, auf Veranlassung 
des ägyptischen Kalifen Ihn- Junis** (ca. 950 — 1008) angefertigt. 
Als später die mathematischen Kenntnisse der Griechen und 
Araber nach und nach ins Abendland gelangten, beschäftigte sich 
dasselbe ebenfalls zunächst mit den praktisch wichtigsten Gebieten 
derselben und es gehört zu den grossen Verdiensten von Johmines 
Begiomontanus*** (1436—1476), dass er die überlieferten Rechnungs- 
vorschriften sammelte und zum Teil ergänzte, so dass er mit 
Recht als der eigentliche Schöpfer der modernen Trigonometrie 
angesehen werden kann. In seinen trigonometrischen Tafeln voll- 
zog Regiomontanus den Fortschritt zur vollständigen Durchfüh- 
rung des dezimalen Systems. Während seine ersten Sinustafeln 
gleich denen seines Lehrers Peuerbach noch eine Vermischung 



* KigeDtlich jieyd^rj (TUVzacc^ oder magna constructio. 

** Er bediente sich, wie es scheint, zuerst der indischen (arabischen) Ziffern. 

*** Eigentlich Johann Müller aus Königsberg in Franken ; kommt schon ala 
12jähriger Knabe an die Universität Leipzig. 



des sexagesimalen und dezimalen Zahlensystems enthalten, ist 
seine dritte Sinustafel auf djen Radius 10^ berechnet. Die Ge- 
nauigkeit ist also dieselbe wie in einer siebenstelligen Tafel, ohne 
dass jedoch eigentliche Dezimalbrüche zur Anwendung kämen. 
Auch in seiner Tangententafel ist die dezimale Teilung durch- 
geführt, indem die „numeri" in Teilen des Radius 100000 an- 
gegeben werden. 

Ein anderes nicht zu unterschätzendes Verdienst von Regio- 
montanus besteht in der Erstausgabe* der sogenannten Ephe- 
merlderiy unter dem Titel „Ephemerides quas vulgo dicunt Al- 
manach" ums Jahr 1474. Es sind dies ebenfalls Tafeln und 
enthalten die Stellungen der Sonne, des Mondes, der Planeten 
und Kometen von Tag zu Tag oder auch in kleineren Intervallen, 
unterscheiden sich somit von jenen naturgemäss dadurch, dass 
der Inhalt derselben veränderlich und deshalb von Zeit zu Zeit 
neu berechnet werden muss. 

Von ganz besonderer Wichtigkeit für die angewandte Mathe- 
matik sind auch die Arbeiten von Georg Joachim Rhäticus (1514 
bis 1576), eines Schülers von Kopernikus, ist er doch der Ver- 
fasser des im Jahre 1596 von Valentin Otho als „Opus Palatinum 
de Trianguhs" herausgegebenen grossartigsten trigonometrischen 
Tabellenwerkes, in dem Sinus, Tangente und Sekante der um 
10" zunehmenden Winkel auf 10 Dezimalstellen berechnet waren. 
Noch genauere -von Rhäticus verfasste Tabellen, den „grossen 
Kanon" gab Bartholomäus Pitiscus, der Urheber des Terminus 
„Trigonometrie", 1613 unter dem Titel „Thesaurus mathematicus" 
heraus. 

Eine weitere Erleichterung erfuhr das praktische Rechnen im 
Anfang des 7. Jahrhunderts durch die Einführung der Logarithmen, 
deren erste Erfindung ohne Zweifel Joost Bürgi zugeschrieben 
werden muss, während der Ruhm der Verbreitung und allge- 
meinen Nutzbarmachung derselben unbestritten dem Lord Neper 
zukommt. 



* Die ersten Ephemeriden stammen allerdings aus einer viel weiter zurück- 
liegenden Zeit, sollen doch nach Ausweis der Keilschriften funde solche schon von 
den babylonischen Astronomen angelegt worden sein, um auf Grund der vorliegen- 
den Beobachtungsreihen in einer für die damalige Zeit durchaus rationellen Weise 
die Orter von Sonne und Mond vorauszuberechnen. 
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Ausgehend von zwei zusammengehörigen Reihen, einer arith- 
metischen und einer geometrischen, verfertigte Bürgi zwischen 
1603 und 1611 eine Tafel der roten und schwarzen Zahlen, deren 
erstere die Zahlen der arithmetischen Reihe — Logarithmen — , 
letztere die der geometrischen Reihe — Numeri — darstellten 
und zwar berechnete er die beiden korrespondierenden Zahlen- 
reihen nach 

X = 10 w und 2/ = 10 • 1,0001 *. 
Trotz Keplers dringender Aufforderung kam Bürgi nicht dazu, 
seine Tafel in Druck zu geben. Dies geschah erst 1620 unter 
dem Titel: „Arithmetische und Geometrische Progress-Tabulen, 
sambt gründlichen Unterricht, wie solche nützlich in allerley Rech- 
nuiigen zu gebrauchen und verstanden werden sol." Beim Er- 
scheinen des Werkes fehlte nun aber seltsamer Weise der „gründ- 
liche Unterricht", so dass die „Tabul" für die meisten Besitzer 
unbrauchbar geblieben ist. 

Inzwischen erschien 1614 die „Descriptio miriflci logarith- 
morum Canonis" des oben erwähnten John Neper (1550 — 1617), 
der 1619, also nach dem Tode des Verfassers, die schon vor der 
„Descriptio" ausgearbeitete „Constructio" folgte. 

Neper ging aus von der Bewegung zweier Punkte.** Der 
erste bewegt sich gleichförmig und die von ihm zurückgelegten 
Wege bilden die arithmetische Reihe der Logarithmen; der zweite 
bewegt sich gleichzeitig, aber so, dass er in der ersten Zeiteinheit 
^ des ganzen Weges, in der zweiten Zeiteinheit ^ des noch übrigen 

Weges, also ^ (1 — ;^) u. s. w. durchläuft. So ergibt sich die fal- 
lende geometrische Reihe der Zahlen mit dem allgemeinen Gliede 

Hierin setzte er nun m = 10^, multiplizierte, um ganze Zahlen 
zu erhalten, die rechte Seite zugleich mit 10''' und erhielt ähnlich 
wie Bürgi zwei sich entsprechende Zahlenreihen, nämlich: 

Xa = n) IJn = 10^ (1 — 4)^ 

wobei er die Zahlen der ersten Reihe Ärtificiales und später 
Logarithmen, die der zweiten Naturales nannte.*** Er besass somit 



* Vgl. R. Wolf: Handbuch der Astronomie pag. 69 T. I. 
** Vgl. A. Sturm : Geschichte der Mathematik (Sammlung Göschen 226). 
*** Siehe Wolf: Handbuch der Astronomie pag. 70 und 71 T. I. 
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in der aus diesen Zahlen zusammengestellten, den Hauptteil seines 
oben erwähnten Canons bildenden „Radicalis Tabula" ein ganz 
entsprechendes Hilfsmittel, wie Bürgi in seiner Progress-Tabul. 

Die praktische Verwendbarkeit der Rechnungszahlen von 
Bürgi und Neper wurde wesenthch durch den Umstand beein- 
trächtigt, dass ihre Grundzahlen nicht mit der Basis des übUchen 
Zahlensystems übereinstimmten.* Diesen Uebelstand fühlte denn 
auch Henry Briggs (1556—1630), als er die Artificiales seines 
Landsmannes kennen lernte, sofort heraus. Die Folge davon 
war, dass er selbst die Berechnung der Logarithmen für die 
Basis 10 unternahm. Die von Briggs aufgelegte „Arithnretica 
logarithmica. Londini 1624" gibt sodann die 14 stelligen Logarith- 
men aller Zahlen von 1 — 20000 und von 90000—100000. 

Die durch Adrien Vlack als „Editio secunda aucta" gegebene 
Neuausgabe vom Jahre 1628 gibt dagegen die Logarithmen aller 
Zahlen von 1—100000, jedoch mit Reduktion auf 10 Stellen und 
überdies auch noch für jede Minute des Quadranten (0 — 90®) 
die Logarithmen der Sinus, Tangens und Secans. 

Mit der Zeit, da der hervorragende Rechner Henry Briggs 
die nach ihm benannten Logarithmen schuf, erreichte eine lange 
Reihe wichtiger Erfindungen ihren Abschluss und es folgte nun 
die Zeit des inneren Ausbaues derselben, welcher hauptsächlich 
darin bestand, dass man die schon bestehenden Tabellenwerke 
den verschiedenen Zwecken, denen sie zu dienen hatten, soviel 
wie möglich anzupassen versuchte. So sind z. B. in der neuern 
Zeit Tafeln beigefügt worden, welche die Möglichkeit herbeiführen, 
aus Log a und Log b direkt, d. h. ohne auf die Zahlen a und b 
zurückzugehen, Log (a+b) zu finden.** Was die Vorgeschichte 
dieser Tabellen anbetrifft, so ist zu erwähnen, dass schon Zecchini 
Leonelli (1776 — 1847) in seinem „Supplement logarithmique. Bor- 
deaux 1803" auf die Theorie und Nützlichkeit derselben auf- 
merksam machte. Auch Gauss pubUzierte im Jahre 1812 eine 
betreffende 5stellig6 Tafel: daher wohl der Name „Gauss'sche 
Logarithmen". 

•Während die Basis der Bürgi'schen Logarithmen 2,7] 81 4593, also bei- 
nahe gleich e ist, findet man für diejenige der Neper'schen Logarithmen annähernd 
1 : c; e = 2,7182 8183... 

** Vgl. das von Dr. J. A. Hülsse herausgegebene Logarithmisch-Trigonometrische 
Vega'sche Handbuch. 
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Die hier mitgeteilten Tataachen über die Entstehung der 
nach und nach zum Gemeingut aller Rechner gewordenen mathe- 
matischen Tabellen beweisen zur Genüge, dass man je und je 
die grösste Sorgfalt auf die Herstellung derselben verwendete. 
Vor keinen Schwierigkeiten scheute man zurück, wenn es galt, 
die Rechnungsvorteile zu vermehren und damit den zu den Rech- 
nungen erforderlichen Zeitaufwand auf ein Minimum zu reduzieren. 

Die Wichtigkeit der Tabellen erhellt wohl am deutlichsten 
aus dem Umfange der darauf bezüglichen Literatur. Würde man 
die gedruckten Tafelausgaben an einem Orte vereinigen, so ent- 
stünde eine ganz ansehnliche Bibliothek, in der die drei grossen 
Ephemeridenreihen des „Berliner Jahrbuch" (1774 — 1906), der 
„Connaissance des temps" (1678 — 1906) und des „Nautical Al- 
manac" (1767 — 1906) für sich allein zusammengenommen schon 
502 Bände zählen würden. 

. § 3. Als drittes Hilfsmittel der Rechentechnik bleiben mir 
noch die mechanischen zu erwähnen übrig. Das Rechnen kostet 
Zeit und Aufmerksamkeit und muss daher mit irgend einem 
Aequivalent bezahlt werden. Es ist deshalb klar, dass man 
von jeher bestrebt war, ausser den mathematischen Tabellen 
Einrichtungen zu ersinnen, die einerseits Zeit gewinnen lassen, 
andei'seits die durch Nachlassen der Aufmerksamkeit erzeugte 
Gefahr von Rechnungsfehlern vermindern. Zu diesen Einrich- 
tungen zählt man die sogenannten Rechnungsmaschinen, und ver- 
steht darunter Apparate, die, wenn auch nicht selbsttätig, so 
doch in der Hand des geübten Besitzers die vom letzteren ge- 
wünschten Resultate auf mechanischem Wege hervorzubringen 
im Stande sind. Die ersten derartigen Einrichtungen waren die 
jetzt so ziemlich wieder in Vergessenheit geratenen Rechenstäbe 
von Neper. Dieselben entsprechen der früheren Uebung, bei 
Multiphkationen sämtliche Teilprodukte einzeln aufzuschreiben. 
Eine andere Art war die nach ihrem Erfinder benannte Gunter- 
Scale, aus welchen nach Vorschlag von Wingate (1627) der jetzt 
noch bei den Praktikern behebte Rechenschieber hervorging. 

Die erste eigenthche Rechenmaschine ist indessen erst die- 
jenige von Pascal. Hierauf folgten ziemhch rasch aufeinander 
solche von Leibnitz, Charles Babbage, Gebrüder Scheutz, Thomas u. a. 
Vor allen diesen verdient die Thomas-Maschine ganz besonderer 
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Erwähnung. Dieselbe gestattet nämlich den Wert eines Aggre- 
gates von der Form aÄ-j-bB-^cC-] , also eines Ausdrucks, 

wie wir ihn bei der Berechnung von Funktions werten aus Reihen- 
entwicklungen und bei den Interpolationen so häufig treffen, 
gewissermassen „vom Blatt" herunterzuspielen, sobald man die 
Faktoren a, Aj b,B, c, C ... vor sich hat. 

Es kann wohl nicht meine Aufgabe sein, hier noch weiter 
auf das Wesen und besonders auf die Einrichtung dieser Ma- 
schinen einzugehen, kommen doch voraussichtlich die wenigsten 
Leser des vorhegenden Schriftchens je in den Fall, dieselben 
wirklich bei ihren Rechnungen verwenden zu müssen. Wer 
dennoch mehr darüber zu vernehmen wünscht, findet eingehende 
Belehrung einerseits in dem von Teignmouth um die Mitte des 
vorigen Jahrhunderts geschriebenen Artikel „On machinery for 
calculating and printing mathematical tables", andererseits in 
dem Artikel „Rechenmaschinen" des Wörterbuchs von Karmarsch 
& Heeren. Was speziell die Rechenmaschine von Thomas be- 
trifft, so verweise ich auf „Franz Releaux: die Thomas'sche 
Rechenmaschine. Freiberg 1862 in 8^". 

Der hier gegebene geschichtliche Rückblick über die Ent- 
stehung und Entwicklung der Hilfsmittel der Rechentechnik lässt 
erkennen, welch grosses Interesse von jeher namentlich den 
mathematischen Tabellen entgegengebracht wurde. Diese Tat- 
sache, dann aber auch die Rücksicht auf alle diejenigen, die sich 
aus irgend einem Grunde mit der Technik des Zahlenrechnens 
— insbesondere mit Interpolation bei mathematischen Tabellen — 
näher zu befassen haben, Hessen in mir den Entschluss reifen, 
diesem höchst wichtigen, aber leider selbst den Studierenden der 
Mathematik oft nur dem Namen nach bekannten Teil des wissen- 
schaftlichen Rechnens, nachstehende Abhandlung zu widmen. 

Den Anfang bilden, um spätem, den Zusammenhang stören- 
den Unterbrechungen vorzubeugen, gewisse Eigenschaften der 
endlichen Differenzen. 





L Endliche Differenzen verschiedener Ordnungen, 



§ 4. Es seien... «1, «2» «s» «4,... irgendwelche gegebene, z. B. 
einer mathematischen Tabelle entnommene, aufeinander folgende 
Funktionswerte. Subtrahiert man jedes a von dem folgenden und 
setzt die so erhaltenen Differenzen ...6i, &2j&3--- in die zweite 
Spalte auf die Mitte zwischen Subtrahend und Minuend, so 
entsteht die sogenannte erste DiflFerenzenreihe, welche nun in 
genau derselben Weise wie die Reihe der gegebenen Funktions- 
werte behandelt werden kann. So erhält man nach und nach die 

Tabelle L* 



«1 


Öl 






«2 




Ci 






&2 




dl 


Ö3 




C2 






&3 




^2 


«4 


&4 


^3 




«5 









in welcher die Grössen ft, c, rf... den Namen: Differenzen erster ^ 
zweiter, dritter u. s. w. Ordnung, oder kürzer erste, zweite, dritte 
u. s. w. Differenzen der Grössen a tragen. Die ganze Tabelle 
bildet somit ein Differenzenschema, wobei man der Gleichförmig- 
keit halber die Ausgangsgrössen a als Differenzen nuUter Ord- 
nung bezeichnen kann. 

§ 5. Es sei nun F (x) irgend eine ganze rationale Funktion 
n^ Grades von x, also 

F(X) = AqX''-}' Äi X'"^ + Ä2X'"^ + + Än-l X + An. 



^ H. Bmns : Grundlinien des wissenschaftlichen Bechnens 1903. 
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Dann ist, wenn man mit J F (x) die erste Differenz von F (x) 
und mit co den konstanten Zuwachs der Argumente oder das 
Tafelintervall bezeichnet 

+ A^ [(»-Dl w ic""- + (w-l)2 («^ a;""^ -I 1- (M-l)„.2 «o""^ « 

+ ^g [(n-2)i <o x"-^H 1- (w-2)„^ «"'^ «+ {«-2)„.2 <«""'] 

+ 

+ ^„., [(2)i o,x^ (2), 6>^ 

+ ^„.1 (Dl '" 

= («)i ^ CO x"'^ 

+ [(«)j Jo <^^ + (*i-l)i -4i «] a;"-^ 

+ [(n)3 A <t>^ + («-1)2 ^1 w^ + (w-2)i ^8 w] ic""^ 

+ • 



+ [(%)„.j A w""' +(«-l)„.2 A ro''-'+(n-2)„.3 J,a>-3 H K2)i ^„.2'"] 



X 



+ (w)„ A 'o" + («-l)„.i A <o""* 4- (n-2),..2 A '"""' H 1- (2)« ^„-2 'o' 

+ (l)i A-i o»" 

d. h. die erste Differenz einer ganzen Funktion n'*" Grades von 
X mit dem höchsten Gliede 

reduziert sich auf eine ganze Funktion (w — 1)'«" Grades von a; 
mit dem höchsten Gliede 

Geht man weiter von i F {x) zu den höheren Differenzen 

/ J' (ic) = A AF (x), 
/ F{x) = A I F (X), 

A"' F(x) ^ A r-^ F (x\ 



* Hienn bedeutet aluremem («),. == -r-^^ — ^ ^— ; ' 

° ^ "" 1-2 '3 k 



i 
) 
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so überzeugt raan sich leicht von der Richtigkeit des folgenden 
Satzes: Wenn F (x) eine ganze rationale Funktion n'*'* Grades 
von X mit dem höchsten Gliede Aq x'^ bedeutet, so ist J"* F (x) bei 
m < n eins ganze rationale Funktion {n — m)^^"* Grades von x mit 
dem höchsten Gliede:* 

Ä^ n (n — 1) (n — 2) (n— 3) (n — m-\-V) r/' x"'"*. 

Hieraus ergibt sich unmittelbar die für die folgenden Betrach- 
tungen wichtige Tatsache, dass 

' j» F (X) = w! A (o\ 
also konstant ist, woraus dann weiter folgt, dass die Differenzen 
von höherer als der n'^'* Ordnung verschwinden. 

Als Beispiel diene die Funktion 3'*'" Grades F (x) = 2 x^ + 
b x^ -\- 7 X -\- l. Berechnet man aus derselben für eine Reihe 
von Argumenten die entsprechenden Funktionswerte, bildet nach 
der in § 4 gegebenen Anleitung die Differenzen der verschiedenen 
Ordnungen und ordnet sie in der in Tabelle I angedeuteten Weise, 
so entsteht die 

Tabelle IL 



1 
X 


F{x) 


I 


II 


III 


1 


15.000 








1,1 


17,412 


2,412 
2,644 


0,232 


0,012 


1,2 
1,3 


20,056 
22,944 


2,888 
3,144 


244 
256 


12 
12 


1,4 


26,088 




268 




1,5 


29,500 


3,412 
3,692 


280 


12 
12 


1,6 


33,192 




292 




1,7 


37,176 


3,984 


0,304 


0,012 


1,8 


41,464 

1 


4,288 







Vgl. Markoff: Differenzenrechnung (deutsch von Friesendorf f u. Prümm) 1896. 
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Der oben ausgesprochene Satz findet sich hier in der Tat be- 
stätigt; denn es ist 

J3 i^ (;:c) == 3! • 2 . (0,1)3 == 0,012. 

§ 6. Betrachtet man die Tabelle II näher, so fällt sofort 
das Verhalten der Differenzen auf: die numerischen Beträge zeigen 
von links nach rechts eine ausgesprochene Abnahme. Infolge 
dessen besitzen die Zahlen in senkrechter Richtung einen ausge- 
sprochenen Gang. Da ein solches Verhalten bei den numerischen 
Tafeln durchaus gewöhnlich ist, so entsteht die Frage nach dem 
allgemeinen Grunde dieser Erscheinung. 

Um diese zu beantworten, ordne ich den im vorigen § ge- 
fundenen Ausdruck für J F (x) nach steigenden- Potenzen des 
Tafelintervalls w und erhalte: 

AF{x) = l{n\ A x'"^ + (w - Dl 4i x'"^ + {n - 2)i A^ o;""^ -\ 

+ (2)i A^^, X + (l\ ^„. J CO 

+ [(n)2 Ao x""-"^ + (n-l)2 A ^"'^ H h (8)2 ^«.3 x+&)^ -^nj ce>« 

+ [^3 Ao X""-^ H hWs A.4^+(3)3 A-3] ''^ 

+ • 

4- (n)^ Ao o)" 
oder: 

AF{x) = F' {X) 'o^-^r^ (:^) |!+ JP-' {x)'^^-\ YF^'-^x)^^ 

+ J^^^)^ (1) 

worin allgemein 

F^ {X) = limes ^ÄJ^ 
Ax=^ ^^ 

die k^ Derivierte der Funktion F (x) bedeutet. Übergehend zu 
den höhern Differenzen, welche aus J F (x), A^ F (x) u. s. w. in 
genau gleicher Weise abgeleitet werden, wie A F {x) aus F (x)^ 
ergeben sich für dieselben folgende Ausdrücke: 

A^ F {X) = Af r' (X) ü? + 4'^ i?"" (X) cü^ -\ h ^f F'' ix) io'' 

A^ F ix) = Af r" (X) 0/ H \- Af F'' (x) -'' 



O) 



J" F (x) = ä'-;^ F" {X) oi'\ 
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in welchen die KoeflBzienten JL folgende einfache Beziehung mit 
einander eingehen: 

(r+ 1) ^41, = k • {Af + Af-'y (2). 

Hiernach sind die /c^" Differenzen in Bezug auf das TafeUntervall 
o) von der Ordnung k. Ist das Intervall w so klein gewählt, dass 
in der Reihe 

J'F(x) = Af F' (X) w' + 4^, i^+^ {x)io''^' + 

das Anfangsglied dem numerischen Betrage nach überwiegt, so 
sind die Grössenänderungen der Differenzen 

JF(x), A^F(x\ J^F(x), 

wesentlich durch die Ausdrücke 

F'{x)a), F''{x)a)^, F'"(x)w^, (3) 

gegeben. Hierdurch ist aber auch zugleich die Ursache aufge- 
deckt, welche bewirkt, dass bei genügend kleinem Intervall die 
Glieder des Differenzenschemas einer tabulierten Funktion den 
oben in Tabelle II beobachteten Gang zeigen. 

§ 7. Die soeben gefundenen Resultate haben vorderhand 
selbstverständlich nur für ganze rationale Funktionen Geltung; 
indessen lassen sie sich mit einigen Modifikationen ohne weiteres 
auf jede, durch eine für alle in Betracht kommenden Argumente 
convergente Reihe darstellbare Funktion übertragen. Diese Ein- 
schränkung ist notwendig ; denn setzt man z. B. für F (x) den 
reziproken Wert von 1 — x und tabuliert die Funktion 1:1 — x 
mit den ganzzahligen Vielfachen von w als Argument, so enthält 
die Tafel auch das Argument o; = 1 mit dem Funktionswerte oo. 
Diese Unstetigkeit breitet sich nun in dem Differenzenschema 
nach oben und unten aus, was zur notwendigen Folge hat, dass 

* Die Bestimmung der Koeffizienten A^ welche bei andern Fragen der Differenzen- 
rechnung, spezieil bei den Untersuchungen über numerische Integration, von grosser 
Bedeutung sind, geht folgendermässen vor sich : Man setzt zunächst k = 2 und be- 
rechnet mit Hülfe der aus (1) sich ergebenden A}^ A!i;\ Ai\ ^?> die 

Werte von A^^ Af\ Af^ u. s. w., wobei noch zu bemerken ist, dass alle Koeffizienten, 
bei welchen der untere Index kleiner ist als der obere = Null zu setzen sind. 
So findet man z. B. 

2 . ^»> == 2 [A? 4- A?] = 2 ^?>; also A?> = A? - 

3 . A? = 2 [A*' + Ai^] = 2 ' |. = 3; also ^f = 1 

u. s. w. , 
Sind alle A^ berechnet, so setzt man k = 3 und bestimmt in gleicher Weise alle 
A^ u. s. w. 



9 
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auf jeder Zeile des Schemas die Differenzen von einer gewissen 
Stelle ab unendlich werden. 

Anhangsweise mag noch auf folgende Tatsache hingewiesen 
werden: Wird das Tafelintervall w-raal kleiner gewählt, so treten 
an die Stelle der Reihe (3) die Glieder 

F'ix)co:n, F^ix)w^:n^, F" (x) cw» : «», (4) 

Man hat es also in der Gewalt, durch Verkleinerung des Inter- 
valls die Differenzen einer höhern Ordnung bis zur Unmerklich- 
keit herab zu drücken. Als Beispiel mögen nachfolgende zwei 
Täfelchen dienen, von denen das zweite gegenüber dem ersten 
halbiertes Intervall besitzt. 



Tabelle III. 



X 


Log X 


I 


TT 


TU 


X 


Log X 


I 


U 


33 


1,51851 


1297 






33 


1,51851 


653 




34 


1,53148 


1259 


38 


+ 2 


33,5 


1,52504 


644 


" 


35 


1,54407 

1 


1223 


— 36 


+ 3 


34 


1,53148 


634 


10 


36 


; 1,55630 


1190 


33 


+ 1 


34,5 


1,53782 


625 


9 


37 

1 


1,56820 


1158 


32 




35 


1,54407 


616 


9 


38 


1,57978 








35,5 


1,55023 







Die Vergleichung zeigt, dass die Differenzen des zweiten 
Täfelchens sich zu den entsprechenden des ersten — rund 
gerechnet — verhalten wie 1:2, 1:4. Wollte man das 
Intervall soweit verkleinern, dass die zweiten Differenzen die 
halbe Einheit der fünften Dezimale nicht übersteigen, so müsste 
der Verkleinerungsdivisor n aus der Gleichung 0,5 • n^ rund 
gleich 35 bestimmt werden, womit sich für n rund 8 ergeben 
würde. 

§ 8. Zurückkehrend zur Tabelle I bilden wir mit denselben 
Grössen a die neue Tabelle: 
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«5 


ß^ 






a* 


ß» 


n 


h 


«8 


ßs 


n 


h 


a» 


ß. 


n 




«1 

• « • • 









welche sich von der ursprünglichen dadurch unterscheidet, dass 
die a in umgekehrter Reihenfolge tabuliert sind. Sucht man 
auch hier die Differenzen der verschiedenen Ordnungen, so er- 
geben sich zwischen ihnen und den entsprechenden der Tabelle I 
folgende Beziehungen: 

^1 = «4 — % = — (% — 0^4) = — &4 
iSg = fl^3 — 0^4 = — («4 — (h) = — h 

ßz = 0.2 (h = («8 O2) = &2 

ebenso /?4 = -^ 61 

u. s. w. 

Tl = ßt — ßl = &4 — &3 == Cs 

U = ßz — ßi ^ h — h == C2 

ebenso 7-3 = ^1 

u. s. w. 

^1 = r2 — ri = — (C3 — C2) = — d2 

ebenso ög = — rfi 

u. s. w. 
Die vorgenommene Änderung lässt sich so auffassen, als ob die 
erste Spalte der ursprünglichen Tabelle um eine Zeile herum 
umgeklappt worden sei. Man kann. also sagen: 

Das Umklappen verändert in allen Spalten des Differenzen- 
Schemas die Beihenfolge der Glieder und in den ungeraden Spalten 
überdies die Vorzeichen, 

§ 9. Die einzelnen Glieder des Differenzenschemas entstehen 
aus den Grössen 0'*^'' Ordnung — Funktion s werten — nur durch 
Addition und Subtraktion. Infolge dessen müssen zwischen ihnen 
notwendig lineare Beziehungen bestehen, die in folgendem auf- 
gesucht werden mögen. Dem Gedankengang Bruns* folgend, 

* Loc. cit. 
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setze ich zu diesem Zwecke mit der willkürlichen Grösse t die 
Reihe ?=+« 

S = y^ap tf = ■■■ (ht-\- Uit' + Ua t^ -\- a^t* -\- a^ t^ -\ (5> 

an, multipliziere dieselbe successive mit 1 — i, (1 — 0^ (1 — ^)^ 

(1 — ty --- und erhalte: 

S(l — t) = «1 t -\- {üi — ai) t^ -{- (Os — ag) t^ + K — az) t^ 

4- (ag — aj ^« + 

S{l—t)^=: ai^ + (a2— 2ai)^»+(a8— 2a2+ai)jJ8 + (a4— 2as 

-\-a,)t^-{-(a^—2a^-\-(H)t^--{- 

Sf (1—^)3 = ait-\-{a2—3(h)t^4-(^9—Sa2-\-Sa^)t^-\-{a^—3(h 

-\- 3 «2 — «i) ^*+ (% — 3a4-|- 3as — a») ^^-|~ 

Sa—t)^= CHt+[(q)oa2-^iq),ck]t^ 

+ [(q)o (h— iq)i cii 4- {q)i «i] t^ 

+ [(g)o «4— (3)i % + (g)« 0^2— (3)8 Ol] ^4 

+ [(3)o %— (g)l «4 + (e)« 0^8 — (3)8 0^8 + (3)4 Ö^l] t^ 
(6) + : 

+ [(3)0 Ö^g+1 — (3)1 <^i + (3)2 ö^g-1 — 

+ (—1)« (3)g «1] ^-^^ 

+ [(3)0 «^+2— (3)1 ^g+l + (3)2 ög — 

+ (—1)^ (q\ a^] t^+^ 

+ 

Betrachten wir nun die in obigen Reihen vorkommenden ... Oi, 
ag, «8, «4, as ... als die Differenzen 0'"' Ordnung des die Tabelle I 
bildenden Differenzenschemas, so erkennen wir ohne weiteres^ 
dass die Differenzen g^^ Ordnung als Koeffizienten der Reihe er- 
halten werden, welche entsteht, wenn wir die Reihe S mit der 
binomischen Entwicklung (1 — 0« = (q)^ — (g)i t + {q)^ t^ — 
(3)3 t^ + multiplizieren. Es ist in der Tat für : 

q = 1, «2 — (^1 — hl as — a2 = 62, «4 — «8 = 08, «5 — cLa = h 

3 = 2, as — 2 «2 + «1 = ^2 — öl = ^1 
«4 — 2 «3 + a2 = &3 — 02 = ^2 
«5 — 2 ch + «3 = *4 — h = Cb 



g = 3, a4 — 3 «8 + 3 «2 — ai = C2 — Ci = di 

«5 3 «4 + 3 as «2 = Cs 6*2 = ^2 
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§ 10. Ausser der Tabelle I sei jetzt noch ein anderes, von 
jenem vollständig unabhängiges Diflferenzenschema angesetzt 
worden : 

Tabelle IV. 



^1 

B, A 

B, 



Beide Tabellen verbinden wir nun Glied für Glied vermittels der 
willkürlichen Faktoren X ujid /^, so dass das Grössensystem 
a^ = Xüp -\- [i Ap^ ßp = Xbp + fJL Bp^ rp = kcp + fx Cp 

u. s. w. 
entsteht. Werden nun die Grössen a, ß^j ... in derselben Weise 
wie die a, b, c ... und A, B, C ... tabellarisch zusammengestellt, 
so entsteht wieder ein DifFerenzenschema. Um dies zu beweisen, 
genügt es, die allerersten Glieder zu betrachten. Es ist 
«2 — «1 = (Aa2 + fJiAi) — (Xüi + fiAi) 
= X (ag — ai) + ^ Ua — Ai) 
= X bt + fi Bi, 
woraus a^ — aj = ßi folgt, d. h. ßi ist in der Tat die zu aj und «g 
gehörige Differenz. Was hier zunächst nur für zwei Komponenten 
nachgewiesen wurde, gilt, wie leicht einzusehen, auch für Ver- 
bindungen aus beliebig vielen Komponenten. Zur Erläuterung 
möge der folgende einfache Fall behandelt werden : Die Tabelle IV 
enthalte dieselben Zahlen wie die Tabelle I, jedoch um eine Zeile 
verschoben, so dass die Beziehungen 

Ap = ctp-^\, Bp == öp-f-i, Cp = Cp^i u. s. w. 
bestehen. Ferner sollen die beiden Multiplikatoren X und // den 
speziellen Wert y besitzen. Dann nehmen die aus der Verbin- 
dung der beiden Tabellen entspringenden Grössen die Gestalt 

«/> = Y (ö^p + CLp+i), ßp =-j (bp + bp^\), Tp =-j ißp + Cp-t-i) 

u. s. w. an. 

Zur Orientierung diene folgendes Beispiel: 
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Tabelle V. 



X 


F(^x) 


I 


II 


in 





5 
6 


2 






1 


7 


8 


12 






14 


14 


18 


12 


2 


21 


36 


24 


12 




40 


38 


30 


12 


3 


59 


56 


36 


12 




96 


74 


42 


12 


4 

1 
1 


133 
194 


98 

122 


48 




5 


255 









In dieser Tabelle bedeuten die fettgedruckten Zahlen die 
vorhin besprochenen Zwischenwerte; dieselben bilden, wie 
leicht zu sehen, ein vollständiges Differenzenschema, die so- 
genannte Zwischentäbellej welche weiterhin von grossem Nutzen 
sein wird. 

§ 11. Bei der Tabulierung von F unk tions werten, die aus 
andern als ganzen, rationalen Funktionen abgeleitet worden sind 
und somit nur mehr oder weniger genaue Näherungswerte dar- 
stellen, werden bei der Abbruchsstelle immer Abrundungsfehler 
begangen, die, wenn auch noch so klein, doch gross genug sind, 
um die Differenzen der höheren Ordnungen ganz bedeutend zu 
beeinflussen. Um die Wirkung dieser Fehler, die wir, um sie 
von denen der andern Kategorie zu unterscheiden, unvermeidliche 
nennen, übersehen zu können, mögen folgende zwei Tabellen, 
von denen die eine mit und die andere ohne Ahrundung erstellt 
worden ist, zum Vergleiche herangezogen werden. 
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Tabelle VI. 













.. - 


X 


i^(a;)=x* 


I 


n 


III 


IV 


1,5 


5,06 


1,49 








1,6 


6,55 


1,80 


0,31 


0,04 




1,7 

1 


8,35 


2,15 


0,35 


0,03 


O.Ol 


1,8 


10,50 


2,53 


0,38 


0,06 


+ 0,03 


1,9 


13,03 




0,44 




0,02 






2,97 


- 


0,04 




2,0 


16,00 


3,45 


0,48 






2,1 


19,45 











Tabelle VII. 















X 


F{x) — x^ 


I 


II 


m 


IV 


1,5 


5,0625 


1,4911 








1,6 


6,5536 


1,7985 


0,3074 


0,0396 




1,7 


8,3521 


2,1455 


0,3470 


0,0420 


0,0024 


1,8 


10,4976 


2,5345 


0,3890 


0,0444 


0,0024 


1,9 


13,0321 


2,9679 


0,4334 


0,0468 


0,0024 


2,0 


16,0000 


3,4481 


0,4802 






2,1 


19,4481 











Als den beiden hier angeführten Tabellen gemeinsames Merk- 
mal erkennt man die Abnahme der numerischen Beträge von links 
nach rechts (§ 6). Während aber die Differenzen der Tabelle VII 
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stetig kleiner werden und schliesslich in Null übergehen, zeigen 
schon die dritten Differenzen der Tabelle VI auffallende Unregel- 
mässigkeiten, welche besonders dadurch charakterisiert sind, dass 
grössere Differenzen mit kleinern abwechseln, eine Erscheinung, 
die dann in der vierten Differenzenreihe die alternierenden Vor- 
zeichen hervorruft und damit nicht nur das Nullwerden der Dif- 
ferenzen verhindert, sondern sogar ein ziemlich rasches Anwachsen 
derselben nach rechts hin bedingt. 

Einen noch bessern Einblick in das Wesen der Abrundungs- 
fehler gewährt die nachfolgende Tabelle VIII, die wir etwa mit 
dem Namen Fehlertabelle bezeichnen wollen. Um dieselbe zu er- 
halten, verbinden wir die Tabellen VI und VII mit Hilfe der in 
§ 10 entwickelten Beziehungen 

Up = kap + fiAp 
ßp = ^bp + fi Bp 
Yp =z X Cp + u Gp u. s. w., 

indem wir in denselben, um unnötigen Nullenbalast zu vermeiden, 
X = 100 und « = — 100 setzen, also dass (fi (x) = 100 {F(x) — F(x)) 
wird. Wir erhalten: 

Tabelle VIIL 



■ 1 

1 

X 


(fi (x) 


I 


II 


TT! 


IV 


V 


VI 


vn 


1,6 


0,25 


















0,11 














1,6 


0,36 


+0,15 


+0,26 


+0,04 










1,7 


-0,21 


+0,45 


+0,30 


-1,20 


1,24 


+ 4 






1,8 


+0,24 




0,90 




+2,76 




9 






t 


0,45 




+1,56 




5 




+17 


1,9 


0,21 


- ■ 


+0,66 




2,24 




+ 8 








+0,21 




0,68 




+ 3 




18 


2,0 


0,00 


+0,19 


0,02 


+0,08 


+0,76 


2 


5 




2,1 


+0,19 


+0,25 


+0,06 


1,16 


1,24 








2,2 


+0,44 

1 


0,85 


1,10 












2,3 


; -0,41 
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Die Beträge der 7. Differenzenreihe sind schon so stark an- 
gewachsen, dass selbst die Zehntel der 7. Differenzen der mit 
den Äbrundungsfehlern tabuherten Funktion F (x) beeinflusst 
werden. Es ist diese Tatsache ein deutlicher Fingerzeig für den 
Berechner mathematischer Tabellen, bei der Tabulierung von 
Funktionswerten möglichst viele genaue Dezimalen mitzunehmen 
und zwar besonders dann, wenn grosse Tafelintervalle häufige 
und weitgehende Interpolationen notwendig machen. 

§ 12. Übergehend zum Studium der zufälligen Fehler,* die 
sich bei der . Berechnung einzelner Zahlenwerte einschleichen und 
durch Unachtsamkeit entstehen, möge folgende Annahme gemacht 
werden. Man denke sich, dass eine Rechnung mit den Grössen 
a durchgeführt werden solle, dass man aber diese Grössen nicht 
genau kenne, sondern statt dessen nur die mit den Fehlern s 
behafteten Zahlen 

« = a + £. 

Dann gilt, wenn man die Differenzen der Grössen a, a und s, 
also «^,4-1 — Up^ ap^x — ap und s^+i — Sp u. s. w. bildet und 
mit einander verbindet, der Satz : Die Fehler der Differenzen 
sind zugleich die Differenzen der Fehler,** An der Hand dieses 
Satzes lässt sich die Wirkung eines einzelnen Fehlers bequem 
verfolgen. 

Wir nehmen zuerst an, dass in der ersten Spalte des 
Schemas nur ein Ghed den Fehler £ besitze, dass dagegen die 
übrigen Glieder der Spalte fehlerfrei berechnet und angesetzt 
worden seien. Dann enthält die Fehlertabelle in der ersten 
Spalte das eine von Null verschiedene Glied s, sonst aber nur 
Nullen, so dass das Diflferenzenschema der Fehler folgende Ge- 
stalt annimmt: 



* Vgl. Rice: The Theory and Practice of Interpolation 1899. 
** Es ist in der Tat Ctp^i — a« die fehlerhafte Differenz, 

a^_j_l — üp die wahre Differenz, 
somit («p+i — ap) — ((lp-\-\ — ö^p) der Fehler der Differenz. Nun ist aber 

(«/»+! «/») (ö^p4-l *CLp) = («/>+! Ö^p-fl) («/> — <^'/i) = i-V+\ — ^i)) 

gleich der Differenz der Fehler w. z. b. w. 
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Tabelle IX. 





< 


.. . 












Hx) 


I 


n 


TTI 


IV 


V 


VI 


VIT 

















+ s 
















+ e 




7s 












+ s 




6s 











+ 5 




5s 




+ 21s 







+ £ 




4e 




4- 15 s 






+s 




3s 




+ 10S 




-35 s 


e 




2s 




+ 6e- 




20 s 






— c 




+ 3£ 




10 s 




+ 35 s 







+ £ 




4e 




+ 15 6 











— £ 




+ 5s 




21s 












+ e 




6e 
















s 




+ 7 s 

















+ « 





Die Faktoren, mit denen s innerhalb einer Differenzenreihe 
multipliziert wird, sind nichts anderes als die gewöhnlichen Bi- 
nominal-KoefiBzienten mit alternierenden Vorzeichen. Der Beweis 
hierfür ergibt sich aus der in § 9 abgeleiteten Formel (6) für 

S {l—t){ 

Setzt man nämlich in derselben alle a bis auf eines, etwa a^y 
welches den Wert e haben mag, gleich Null, so geht die rechte 
Seite über in 

(g)o £ — (g)i s + (?)2 s — (g)3 ^ + + (— 1)^ {q)i s. (7) 

Aus dem Verhalten der Zahlen in Tabelle IX ist zu entnehmen, 
dass die V^irkung des Fehlers s an Ausbreitung und Grösse zu- 
nimmt, je weiter man mit der Bildung der Differenzenreilien 
fortschreitet und dass femer die stärkste Fehler Wirkung in jeder 
Spalte stets in derjenigen Zeile stattfindet, auf welcher der Fehler 
in der ersten Spalte begangen worden ist. Demnach können tvir 
aus den Werten der Differenzen verschiedener Ordnungen nicht nur 
einen Schluss auf das Vorhandensein eines Fehlers ziehen, sondern 
sogar seinen Ort, sowie bis zu einem geioissen Grade wenigstens 
auch seine Grösse bestimmen. 

Zur Erläuterung einige Beispiele! 
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Als erstes wähle ich die aus 

F(x) = 10:^3 — 101 a;2.— 109 a; + 1799 
berechnete 

Tabelle X. 



X 


Fix) 


I 


n 


m 


c 


TTT+c 


Fix) 





1 

+ 1799 


200 










+ 1799 


1 


+ 1599 




142 








+ 1599 


1 




342 




+ 60 




+ 60 




2 


+ 1257 




82 








+ 1257 


1 




424 




+ 55 


+ 5 


60 


. 


3 


+ 833 




27 








+ 833 






451 




+ 75 


15 


60 




4 i 


+ 382 




+ 48 








+ 387 






408 




+ 45 


+ 15 


60 




5 ' 


21 




+ 93 








21 


I 

1 




-310 




+ 65 


5 


+ 60 




6 


331 1 


152 


+ 158 








331 


7 


483 

1 












483 



Nach den in § 5 angestellten Untersuchungen über ganze 
rationale Funktionen sollten sämtliche Differenzen dritter Ord- 
nung den Wert 60 annehmen. Dem gegenüber zeigen aber die 
numerischen Beträge der obigen Tabelle deutiich hervortretende 
Abweichungen, welche sich nur dadurch erklären lassen, dass 
man annimmt, die Tabulierung sei eine fehlerhafte. 

Die Auffindung des Fehlers gestaltet sich in unserem Bei- 
spiel um so einfacher, als wir die genauen Beträge der dritten 
Differenzen schon zum voraus kennen. Der einzuschlagende Weg 
ist kurz folgender: 

Neben die erste abweichende Differenz notieren wir in der 
mit c* überschriebenen Spalte den Unterschied zwischen der rich- 
tigen Differenz 60 und der fehlerhaften 55, also 5, bilden damit 
die übrigen Glieder der binomischen Entwicklung (7) für q = B 
und e = 5 und tragen dieselben ebenfalls in c ein. Nun be- 



Abgekürzte Bezeichnung für Correction. 
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rechnen wir die Spalte III + c. Da in derselben sämtliche Be- 
träge gleich 60 werden, so dürfen wir mit Sicherheit schliessen, 
dass der begangene Fehler fünf Einheiten umfasst. 

Um nun noch den Ort zu ermitteln, wo der gefundene Fehler 
sich eingeschUchen hat, erinnern wir uns des im Anschluss an 
die Tabelle IX ausgesprochenen Satzes. Nach demselben kann 
nur 382 der fehlerhafte Funktionswert sein. 

Als zweites Beispiel diene folgende, der Schlömilch'schen 
Logarithmentafel entnommene 

Tabelle XL 



X 

40 


Log X 


I 


TT 


TTT 


c 


1 

III+c 

1 


1 
Log X 

L 

1,60206 


1 1,60206 










1 

1 


! 
i 


1072 












1 


41 


' 1,61278 

1 


1047 


25 












1,61278 


42 ; 

1 


1 1,62325 


1022 


25 


+ 1 




^ 


1 


, 1,62325 


43 


1,63347 

1 


998 


24 


+ 2 




+ 


1 

2 


1,63347 

r 


44 


1,64345 




22 




M 






1,64345 






976 




+ 5 


— 4 


+ 


1 


1 


45 


1,65321 


959 


17 


12 


+ 12 




1 




■ 1,65821 

i 


46 


1,66280 




29 










1,66376 






930 




M3 


-12 


1 


1 


1 


47 

1 


1,67210 




16 










1,67210 






914 




2 


f 4 


+ 


2; 




48 


1,68124 


896 


18 


1 






1 

11 


: 1,68124 

1 


49 


1,69020 


877 


19 








1 

! 


1,69020 


50 


, 1,69897 

1 












1 


1,69897 



Bildet man, wie dies in der Spalte III -|- c geschehen ist, die 
Summen aus den dritten Differenzen und den in oben beschrie- 
bener Weise gefundenen Correctionen, so erhält man die Diffe- 
renzen, wie sie sich ergeben würden, wenn der Funktionswert 
Log 46 richtig, also mit 1,66276 angesetzt worden wäre. 



A 
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§ 13. In den beiden soeben vorgeführten Beispielen wurde^ 
entsprechend der der Tabelle IX zu Grunde gelegten Annahme 
die Existenz von nur einem zufälligen Fehler vorausgesetzt. Von 
dieser Einschränkung wollen wir nun Umgang nehmen und zeigen, 
wie man auch jetzt noch im stände ist, sowohl die Grösse als 
auch den Ort der Fehler zu ermitteln. 

Als erstes Beispiel dieser Art diene die durch Tabulierung 
der Funktion F{x) = 2x^ -\ dx^ + 7 x + 1 entstandene 



Tabelle XIL 



i 
X 


Fix) 


I 


11 


ITT 


Ci 


III+c, 


Ci 


III+c, 


Fix) 





1 














1 


1 






14 












1 


1 


15 




22 










1 


15 


2 


51 


36 


34 


+12 




+12 




+12 


51 


3 


121 


70 


46 


+12 




+12 




12 


121 


i 


237 


116 


53 


+ 7 


-f- 5 


+12 




12 


237 


5 


406 


169 


80 


+27 


15 


+12 




12 

1 
1 


411 


6 


655 


249 


92 


+12 


+15 


+27 


15 


12 


655 






341 




—28 


5 


33 


+45 


12 




7 


996 




64 












9S1 


8 


1401 


405 


121 


+57 




+57 


45 


12 

i 


1401 






526 




3 




3 


+15 


12 




9 


1927 




118 












1927 


10 


2571 


644 


130 


+12 




+12 




+12 


2571 






774 












1 


11 


3345 














3345 

1 



Die erste von 2 • 81 = 12 abweichende Differenz ist 7. 
Man setzt somit rechts in die Spalte Ci die Verbesserung -{- 5 
und darunter der Reihe nach — (3)i- 5, -f (8)2- 5, — (8)3- 5. 
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Die grösste Pehlerwirkung findet nun von — 15 zu + 15 statt, 
was darauf schliessen lässt, dass 406 ein fehlerhafter Funktions- 
wert ist. Jetzt berechnet man die Spalte III -f c^ und erhält die 
verbesserten Differenzen. — Die Tatsache, dass die letzteren nur 
teilweise den Wert 12 annehmen, lässt erkennen, dass ausser 
dem fehlerhaften Funktionswert 406 noch ein anderer ebenfalls 
fehlerhaft sein muss. Eine ähnliche Untersuchung wie die vorige 
lässt auch diesen auffinden. Es ist der Funktionswert 996. Die 
Übereinstimmung der Differenzen mJR-{- Ci + c^ beweist zudem, 
dass die Verbesserungen + 5 für 406 und — 15 für 996 wirkUch 
diejenigen sind, welche uns die richtigen Werte 411 und 981 liefern. 

Als zweites und zugleich letztes Beispiel denken wir uns die 
Sonnendeklinationen für den Monat April 1906 von zwei zu zwei 
Tagen berechnet und tabuliert. (Siehe Tabelle XIII.) 

Dieselben Überlegungen, die uns bei der Betrachtung der 
Tabelle XII zur Entdeckung der fehlerhaften Funktionswerte 
fahrten, lassen solche auch hier wieder auffinden. Allerdings 
muss bemerkt werden, dass die Ermittlung hier bedeutend mehr 
Arbeit und Sorgfalt erfordert, als dort. Die Ursache hegt we- 
niger in der grössern Fehlerzahl, als vielmehr darin, dass es 
dem Rechner wohl nicht geUngen wird, schon das erste Mal die 
richtigen Verbesserungen zu finden. Die Neigung der dritten 
Differenzen, den Wert — 2,7 " anzunehmen, wird ihn ganz sicher 
veranlassen, als erste Verbesserung in Ci nicht + 5)5", sondern 
eher-f- 5,4" zu wählen, um so mehr, als diese Annahme in der 
Spalte 1114- Ci zu den scheinbar noch bessern Werten — 2,7", 
— 2,7", — 2,4", — 2,6" führt, als es die in der Tabelle XIH 
notierten sind. Erst am Schlüsse seiner Rechnung angelangt, 
wird der Rechner bemerken, dass jene gewählte erste Verbes- 
serung mit 5,4" zu klein angenommen wurde. Diese Tatsache 
hängt damit zusammen, dass die erste Verbesserung die zweite 
und auch die dritte Verbesserung beeinflusst. Diese Beein- 
flussung existiert natürlich auch bei der Tabelle XII; allein 
dort sind nach § 5 die dritten Differenzen konstant, während 
dies bei obigem Beispiel nicht der Fall ist, ja wegen der unver- 
meidlichen Abrundungsfehler überhaupt nicht der Fall sein kann. 
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Deklination 
der O 


o 


5 3 58,6 

5 49 50,1 

6 35 17,0 

7 20 16.7 

8 4 46,5 

8 48 44,1 

9 82 6,7 
10 14 51,9 

10 56 57,p 

11 38 19,1 

12 18 55,6 

12 58 43,7 

13 37 40,6 

14 15 43,6 




•-< + 


2.7" 
2,6 
2.7 
2,3 
2,8 
2,4 
2,7 
2,9 
2,6 
2,8 
2,8 
-2,7 




;S* 


vO 00 00 vO 

+ + 1 




6^ 

+ ;? 
t^ + 


2,7" 

- 2,6 

- 2,7 

- 2,3 

- 2,8 
2,4 

- 8,7 
+ 15,1 
-20,6 
+ 3.2 

- 2,8 

- 2,7 




^ 


- 3,6" 
+ 10,8 
-10,8 
+ 3,6 


• 

M 
M 


«4 
+ 

1— ( 


- 2,7" 
2,6 

- 2.7 

- 2.3 
+ 0.8 
-13,2 
+ 2,1 
+ 11,5 
-20,6 
+ 3,2 

- 2,8 

- 2,7 


'S 




+ 5,5" 
-16,5 
+ 16,5 
- 5,5 


in 




5: 

r^-^000000C<l^iOvO<N00t^ 
<N OO" od 00 vo od Of »^ O od (N (N 

1 + T + 1 + + 1 + 1 1 




i 


— 21,9" 
24,6 
32,7 
18.9 
37,7 
31,4 
44,6 
42,5 
31,0 
51,6 
48,4 
51,2 

-53,9 




1— » 

1 
1 


er: 

•«1 


45 51,5 
45 26,9 
44 54,2 
44 35,3 
43 57,6 
43 26,2 
42 41,6 
41 59,1 
41 28,1 
40 3b,5 
39 48,1 
38 56,9 
38 3,0 




Deklination 
der O 


O 


5 3 58,6 

5 49 50,1 

6 35 17,0 

7 20 11,2 

8 4 46,5 

8 48 44,1 

9 32 10,3 
10 14 51,9 

10 56 51,0 

11 38 19,1 

12 18 55,6 

12 58 43,7 

13 37 40,6 

14 15 43,6 




April 
1906 


^H 


00lO^-c^-JC2£E:2c3S3SS5a 




IL Die Formeln und Methoden der Interpolation. 



§ 14. Es sei y = F (x) eine für die Argumente Xo» ^i--.^« 
eines Intervalls von x' bis rz;" numerisch bekannte Funktion. 
Dann besteht die Aufgabe der Interpolationsrechnung darin, den 
Wert von F [x) für irgend ein anderes Argument desselben Inter- 
valls mit Hilfe einer der numerischen Berechnung leicht zugäng- 
hchen Funktion näherungsweise zu bestimmen. 

Wir behandeln in folgendem nur den Fall, dass y = F(x) 
durch eine, gewissen Anforderungen genügende, ganze rationale 
Funktion rj = 0{x) ersetzt werde. 

Für diese wählen wir den Ausdruck 

(1) 7J = ÄoX^' +Ä^ X'"-^+ + A 

und knüpfen daran die Forderung, dass in den n -{- 1 Stellen 
Xo, Xi, Xn die Identitäten 

(2) ^0 = 2/0, ?i=2/i, Vn = yn 

stattfinden. Um diesen Bedingungen zu genügen, müssen die 
Ao, Ai, die Gleichungen 

A rr? + ^1 x^-^ + -\- An = yo, 

Ao x^ + A^ xf-^ + + An = «/i, 

Ao xl + Ai xl~^ + + An = yn 

befriedigen. Da diese Gleichungen Unear in den Unbekannten A 
sind, so würde die Berechnung derselben gar keine Schwierig- 
keiten bieten. Nun kommt es uns aber nicht sowohl auf die 
Bestimmung der einzelnen Werte der A, als auf ihre durch (1) 
gegebene Funktion an ; und diese können wir direkt berechnen, 
indem wir aus dem gleichzeitigen Bestehen der Gleichungen 
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rj — An — An-l X — 

yo — An — A-i Xo — 



— A^x"" = 0, 

— Ao x^ = 0, 



yn — An — An-l OOn — — Aq X^ =0 

weitere Schlüsse ziehen.* Zunächst können wir diese n -{- 2 
Gleichungen als ein vollständiges System von homogenen Gleich- 
ungen der n-{-.2 Werte 

auffassen. Nun aber muss für ein solches immer die Determi- 
nante ihrer Koeffizienten null sein. Man hat also 



V 


1 


X 


X' ... 


x" 


Vo 


1 


\ 


xl... 


/Y»tl 


Vi 


1 


X, 


X* • - - 


x\ 


Vn 


1 


\ 


xl... 


K 



= 



oder nach den Elementen der ersten Spalte entwickelt 



(3) 



X X* . . »X^ 

Ix^-'-x^ 


V 


Ix ...af 

JL X4 - - .Xa 

lx^...x?^ 


2/0+ 


Ix ...rc" 

\x^...x^ 

J X2...X2 


«/!-• 


..+ 


lx...x^ 

n n 


jlX^ . . 'X^ 
n n 


lx...xl 

n n 



Ix 



Ix. 



0' 



_ 1^1 



/Y«n 



IXn- 



//.« 

^n-1 



2/n=0^ 



Beachtet man ferner den Satz, dass eine Determinante ihr Vor- 
zeichen ändert, wenn man zwei Parallelreihen miteinander ver- 
tauscht, so geht (3) über in 



1Xq...x^ 




Ix ...iC" 




l%--^o 




1 Xq. . .a^ 


X Xt • • .x^ 




J. X^ • • .Jb.. 




J. iC . . .X 




X x^ . . .^^ 


1 Xn' • 'Xn 


'i— 


1 ^2" • •*^2 


2/0 + 


IX^'-'X^ 


2/1 + --H- 


1 Xn. - -i^ 


Ix.'.xl 

n n 


n tt 


1 iZ/ ... x^ 

n n 


loj ...o;** 



2/«. 



Dividiert man diese Gleichung durch den Koeffizienten von tj und 
rechnet die einzelnen Determinanten aus,** so erhält man nach 
einigen Umformungen den Ausdruck 

i trr n 

[X — Xq) [X — X^) . . .{X — Xi — \} \X ^""Xi-L-ij ... \X Xn) 



~ Z^ (x^ — , 

1 = 



^) [Xl Xl) — (^i — Xi—i) \Xi — Xi-^x) — \Xi Xn) 



'Vi. 



* Vgl. Netto: Elementare Algebra, Akademische Vorlesungen 1904. 
** Siehe Anmerkung * auf der nächsten Seite. 
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welcher durch Einführung der Funktion 

(p (x) = (x — Xo) (x — Xt) {X—X2)'.-{X-^0Ct) 

übergeht in den folgenden einfachem : 



(X—Xn) — TT (X 
1 = 



-Xi) 



(4) 



Zj <P i^i) • (^ — ^i) 
i= 



Die Herleitung von (4) zeigt uns, dass es eben nur diese 
eine Funktion gibt, welche den in (2) aufgestellten Bedingungen 
entspricht. Die Formel heisst die Lagrange'sche Interpolation»- 
Formel. Um die Richtigkeit derselben einzusehen, genügt es zu 
zeigen, dass z. B. für o; = a;o, 7 = Vo wird. 

Setzen wir somit in (4) x = Xq, so erkennen wir sofort, 
dass alle Summanden mit Ausnahme des einen, der mit j/o multi- 
pliziert ist, verschwinden; und dieser erhält den Koeffizienten 1.** 

Zahlenbeispiel. Gegeben: log sin 20<* 31' = 9,5446630, 

log sin 20« 34' = 9,5456745, 
log sin 20« 36' = 9,5463472, 
log sin 20» 37' = 9,5466832, 
log sin 20» 40' = 9,5476893. 

Man berechne log sin 20« 38'. 



* Die Ausrechnung der Determinanten von der vorgelegten Form geht wie 
folgt vor sich: Es ist 



1 Xo ^ 3^ 

i. Xi CCi Xfi 

X Qu% jC% CC% 

i- Xz Xz Xz 



= (Xi— Xo) (X,— Xo) {Xz—X^ 



** Es ist zunächst 
nach Hopital 
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J- Xi — Xo Xi (Xi — Xo) Xi (Xi — Xo) 
1 Xj Xo Xj (Xj — Xo) X^ (Xj Xo) 
1 X8— Xo X, (X«— Xo) xJ (Xg— Xo) 

1 

1 Xj Xi Xj (,Xj — Xi) 
1 Xj— Xi X, (Xj— Xi) 
= {Xi Xq) {X2 Xq) {X^ Xq) 

(^8— ^1) (^3 — ^1) 



— (Xi—OCo) (X,— Xo) (X,— Xo) 



(Xt—Xo) (Xt—Xo) (^t—^a) 

{X^—Xt) (Xt—Xi) 



L Xx Xi 

1 Xtoi 
1 Xtoi 

1 Xt 
1 Xg 



limes 



^' (Xo) (X— Xq) 
(fix) 



= -S-, also tmbesümmt. Nun ist aber 



^' (Xo) (x— Xo) 



_ /(^) _ 1 



X -rxö f (Xo) 
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Die Lagrange' sehe Formel lautet för fünf gegebene Grössen : 



>? = 



+ 
+ 
+ 
+ 



\CC — 00\ 






{Xi — Xq 
[X — CCq 



(x ■ 



-Xq 
-Xq 



(X^ — OOq 
[X OCq 



\X^ Xq 



{X — Xt) ix — Xs) {X — SO4) 



(X — Xt) (X —00s) (^ — ^a) 



(Xf — Xi) (Xt—OCs) (5Ci — X4) 
(X — Xi) (X — Xz) (x —^Xi) 



(x^ — Xi) (X^ — Xs) (X^-^Xi) 
\X — X^j (X — X^f \X — X^f 



(X^ — OOi) (X^ — X2) {00$ — X4J 
(X — Xi) (x — X2) (x — x^) 



Vo 



Vi 



y% 



2/3 



2/4. 



^etzt man hierin statt a\), Xi, «j, Xz, x^, x die Winkel 20» 21', 
20« 34', 20« 36', 20« 37', 20« 40', 20« 88' und fflr t/,, yuyi,ys,y4. 
die ihnen entsprechenden Funktionswerte, so erhält man: 



^ = — 



+ 



-^ • 9,5446680 + ^ ■ 9,5456745 — J • 9,5468472 
405 27 5 
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9,5466832 + ^V ' 9,5476898 

Ol 



= — 0,1885365(5) -f 2,4748045(0) — 13,3648860(8) 

+ 19,8005281(2) + 0,8251089(5) 
== 9,5470189. 

§ 15. Die Interpolationsformel von Lagrange zeichnet sich 
wohl durch ihre Gleichartigkeit in Bezug auf die Argumente 
aus, aber sie ist in dieser Form doch wenig geeignet zur Inter- 
polation selbst und zwar vor allem deshalb, weil man bei Hin- 
zunahme neuer gegebener Grössen nicht nur neue GUeder hinzu- 
fügen, sondern die frühern auch ändern muss. Sie gibt z. B. 
bei einer gegebenen Grösse 

bei zwei 



V = 



/y» /y« 



-z — z^* yu 



X\ Xq 



bei drei 



rj = 



[ßj X\ ) {X X2 } 
(Xq OOi) (Xq X2) 



^0 + 



[X »^0/ \P^ '^2/ 

\Xi — Xq) [Xi- — X2) 



Vi 



+ 



\X Xq) \X 0(j\) 

kX^ — Xq) [X^ — Xi) 



y^ u. s. w. 
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Die Lagrange'sche Formel bedarf somit einer Umformung 
und zwar in dei* Weise, dass man hei Hinzunahme neuer gegebener 
Grössen nur neue Glieder hinzuzufügen braucht, ohne die vorher- 
gehenden ändern zu müssen. 

Wie dieser Forderung am rationellsten entsprochen werden 
kann, möge in folgendem gezeigt werden. 

§ 16. Wir setzen nach Gauss* in (4): 
-^0 = 2/0, 

/>» /y» /y /y 

^1= - — —' 2/0 + -z — —' Vi 1 



OuQ CCi 



«Cj Ji/Q 



^ {X —Xi) {X —X2) , {X —Xo) (X —X2) . (X —Xo) (X ~x^) 
•^2= Tir— TTTi: ITT 'y^-T-FZ ^\ /^ ITX 'Vl-rTZ ZrrTZ, ITT- 2/2> 



(^0 — ^1) (^0 — ^2) \^1 — ^^) (^1 — '^2) (^2 — ^) (^2 — ^1) 

-er \P^ ^1) (^ «^ä) (^ ^S) I (^ ^q) (^ X2) \X 00%) 



+ 



(Xo — ^^1) (^0 ^2) (^0 ^^) 
(^2 — ^0) (^2 — ^\) (^2 — ^^3) 



(^1 — ^0) (^1 ^^2) («^l — ^^) 

^' ~^ te-^o) (^3-^) (i^Cs-^r^) ' ^" 



Hieraus folgt durch successive Subtraktion: 

-^0 = 2/0, 

'2/0.2/1 



^1 — Jlq — (^ — X^ 



.Jx.2~ — -^1 — [X — Xq) [X — Xi) 



+ 

OüQ Xi Xi Xq 

2/0 



(Xq Xi) (Xq — X2) 



+ 
+ 



2/1 






(^1 Xq) {Xi X2) 

Vi 



+ 



(^ ^1) (^ X2) \Xq x^) 

Vi 



(aJi— aso) {oCi—Xi) {x^—Xi) ' (Xi—Xo) («2— «i) (Xi—Xs) 



+ 



(^3 ^0) te ^^1) (^ — ^2) , 



* Gauss^ Werke III. Vgl. auch Markoff, Differenzenrechnung (deutsch voa 
Friesendorff und Prümm) 1896. 
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und somit durch Addition 

? = 2/0 -\- ix — OOo) [oco, ^i] -{- (x — Xo) (x — x^) [odq^ x^^ x^ 

<5) + (^ — ^o) {x—Xi) (X — X2) [xo^' Xu X2, Xs] H 

Hierin ist: 



(6) 



2/0 



2/1 



Xq. Xi — • 



\Xo^XijO02J = 
I ^, Xij ^2) Xq\ 



2/0 



Vi 



Vi 



\X^—Xi) (OCq—X^) \X\—Xq) {Xi—X^) \X2~Xq) \X2'~Xi) 



(oOo-Xt) (Xo-Xi) (Xo-oJb) (:3:5i-^) (^1-^2) (^-^) 
iX2-oco) (x^-Xi) (x^-x^y (OOz-OOo) (Xb-Xi) (Xb-X^)' 



Betrachten wir die Beziehungen (6) etwas näher, so erkennen 
wir, dass sie sich nicht ändern, wenn die Buchstaben Xo, ooi^ x^, 
a^, ... gleichzeitig wie die Buchstaben 2/07 Pu ^2? 2/3? ••- unter sich 
vertauscht werden. Gestützt auf diese Tatsache findet man aucl^ 

2/1—2/0 



(7) 



[Xo, Xi] — 



[Xq^ X\^ X%\ 



[^0} X\^ ^2, 0C%\ 



X\ ""^Xq 

[OOij OC2] — [Xq^ Xi] 
OC2 Xq 

[OOij ^, OOs] — [Xq^ X\^ x%\ 
X% Xq 



In der Form (5) ist die Interpolationsformel bereits von 
Newton gegeben worden. Betreffs der Anordnung erhalten wir 
das folgende Schema: 

Tabelle XIV. 



^ 
% 


Vi 




[xq, Xi, a^] 


L^^ ^f ^> ^s] 




« 




^ 


Vi 


fe^sl 


[hf Xi, Xsl 


[Xj, X2f Xq, X4J 


lXo> Xi, Xg, Xg, X4] 


[X^j Xj, Xg, X3, 


X4, Xß] 


H 

^A 


Va 




[xa, Xs, X4] 


1 

[X2, Xs, X4, Xß] 


[Xj, a^, Xg, X4, X5J 






2^5 


Vs 















/Ol 
VS 2 



0/, 
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§ 17. Aus den Untersuchungen des vorigen Paragraph^» 
geht hervor, dass die Argunaente ocq, ^i, ^2, a^, ... in jeder belie- 
bigen Reihenfolge genommen werden können. Wählen wir also 
z. B. die Substitution 

2 3 4 

4 15 

so geht die Formel (5) über in 

i ? = 2/3 + (^-^3) [^, 0(k] + (OO-X^) (X-X2) [002, ^3, ^4] 
-j- (X-Xs) {X-OO2) (OC-Xi) [OCijX^^OCsjXi] 

(8) I 4" (^-^3) (aJ-a;2) (x-Xi) (x-Xi) [ooi, 0C2^ x^^ ^4, x^'\ 

-f- (X-iCa) (X-X2) {X-X^) {X-OOi) {X-Xf,) \Xq^ Xi^QC^Xs, X^^ X^'\ 

+ 

Bemerkenswert ist noch die zu (8) symmetrische Formel,* nämlich 

9 = 2/2 + {X-Xi) [Oh, Xz] -\- {X-X2) (X-Xz) [0(h,0(h, x^l 
+ {OO-X2) (X-OOs) (X-OCi) [Xi,OC2,003,Xi] 

(9) ! + {X-0C2) (oo-o(k) (x-Xi) (x^Xi) [Xo, ooi,ah,Xs, x^] 

-f- (X-OO2) (OO-Ok) {Xr-Xi) (OO-Xi) (X^Xo) [Xo, Xi, X2j Xs^ Xi, iZJö j 

+ • • 

Dieselbe wird aus (5) mit Hilfe der Substitution 

/O 1 2 3 4 5 
\2 3 1 4 5 

in genau gleicher Weise abgeleitet, wie (8) aus (5). 

Zahlenbeispiel. Man berechne die scheinbare Deklination des 
Mars für den 13. April 1906 unter der Voraussetzung, dass sie 
für den 9., 10., 12., 14. und 15. April bekannt ist. 

Nach dem Nautical Almanae für 1906 wird 

filr den 9. AprU ^ = 17« 5' 38,3", 

10. „ d = 170 18' 20,2", 

12. „ d = 170 43' 14,5", 

14. „ 5 -- 180 7' 28,8", 

15. „ d == 180 19' 20,8". 

Bildet man mit diesen Werten die der Tabelle XIV entsprechende 
Zahlentabelle in der Weise, dass man für Xq, Xi, x^, o^s, X4^ der 



) 



w 



» 



M 



n 



w 



w 



» 



yy 



* Die Symmetrieaxe ist eine Horizontale der Tabelle XIV und liegt in d«r 
Mitte zwischen den beiden Werten y^ und 1/3. 



—IV- .1«_? 



/ 
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Reihe nach die Zahlen 9, 10, 12, 14, 15 und fllr y«, Vu Vn 2/8. 2/4 
die ihnen entsprechenden DeJdipationen setzt, so ergibt sich die 

Tabelle XV. 



Xo— 9 
»1 = 10 


«/e = 17« 5' 38,3" 
Vi — 11 18 20,2 


12' 41,90" 
12 27,15 


— 4,92" 


0,01 " 


Xi — \2 
Xt — U 


y^ — \l 43 14,5 
^8 — 18 7 28,8 


12 7,16 

11 52,00 


5,00 
— 5,05 


0,01 


Xi 15 


yi = J8 19 20,8 









Nun wird nach (6) 

, = 170 5' 38,3" + 4 . 12' 41,90" — 4 • 3 - 4,92" — 4.3.1- 0,01 " 
= 170 5' 38,3^' + 50' 47,60" - 59,04" - 0,12" = 17» 55' 26,7(4)"; 

nach (8) 

Tj = 18« 7' 28,8" — 1- 12' 7,15" + 1.1. 5,05"— 1.1.2.0,01" 
= 18« 7' 28,8" — 12' 7,15" + 5,05"— 0,02" = 17« 55' 26,6(8)" 

und endlich nach (9) 

, = 170 43' 14,5" + 1.12' 7,15" + 1 • 1 . 5,00" +1-1.3. 0,01" 
= 170 43' 14,5" + 12' 7,15" + 5,00" + 0,03" = 17« 55' 26,6(3)". 

Der Nautical Almanac gibt ^ = 17» 55' 26,7". 

§ 18. Die Formeln (5), (8) und (9) wollen wir nun, entsprechend 
der uns gestellten Aufgabe, auf den speziellen Fall anwenden, dass 
die Argumente o;©, ^1? ^2? •-. die Glieder einer arithmetischen 
Reihe erster Ordnung sind, also dass 

ist. Dann wird, wenn man Xq =^ a setzt, 
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und analog 



Xi = a -\- io^ 0C2 = a -\- 2wj ' 



X—\ = a — iO^ X—2 ^^ d — 2 fe>. 



Um ' däbßi mit der ziemlich allgemein eingeführten Encke* sehen 
Bezeichnungsweise für die diesen Argumenten entsprechenden 
Funktionswerte 
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••'1 y-2, 2/-1, 2/0, 2/i, «/2, ... 
Übereinstimmung zu schaffen, setzen wir allgemein für den p^^ 
Funktionswert das Symbol 

{a -{- p (0^ 0) 
und finden so für p = ... — 2, — 1, 0, 1, 2, ... 

. .., 2/_2 = (a — 2 CO, 0), y-i = (a — w, 0), yo = (a, 0), 
2/1 = (ö^ + ^«>, 0), ^2 = (a + 2 fe>, 0), ... 
Subtrahiert man jede dieser Grössen von der folgenden, so 
ergeben sich die Differenzen erster Ordnung, deren allgemeine 
Form durch die Beziehung 



(a +p + 1 .0, 0) — (a + p .e>, 0) = (a +p + I CO, D* 

festgelegt sein mag. In gleicher Weise erhält man auch die 
Differenzen zweiter, dritter u. s. w. Ordnung aus 

(a -fp + -j CO, 1) — {a-{-p — -i- CO, 1) = ia-{-p co, 2), 
(a+^ + 1 CO, 2) — (a + p CO, 2) = (a +"7+? ^^' 3)' 



Nach diesen Festsetzungen geht die Tabelle XIV über in die 
folgende : 



* Die Schreibweise ci~\-p-^ — co soll daran erinnern, dass die durch (a -\-p -|- — ^ o 1 ) 
angedeutete erste Differenz aus den beiden Funktionswerten erhalten wurde, deren 
Argumente a -\- p -\~ — (q zum arithmetischen Mittel haben. 
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Setzt man noch ^ 5=5 a -j- to und rj = / (a -|- tw\ so geht mit 
Rücksicht auf obige Tabelle die Formel (5V über in die Newton' sehe 
Formel für äquidistonte TniervaUe 

/(« + ^..) = (a,0) + ^(a + f.,l) + ^";^^^^-f (a + ..,2) + . 



<u 



= (a,0) + Ma + -i-w,l) + 



+ 



<(< 



t{t — \) 
1-2 



(a + «, 2) 
■2) 



1 •2-3 



(a+|w,3) + 



oder kürzer 

r 

/■ (a -f- ^ "') 
(10) 



(a, 0) + (Ol (a + Y «>, 1) + (Oa (a + «, 2) 
+ (^)8 (a + I '", 3) + (t\ (a + 2 w, 4) 
+ it\ (a + T w, 5) + (Oe (a + 3 «, 6) 



+ (^7 (a + |<«, 7)+---. 
In gleicher Weise erhält man statt (8)* 

- -5- w, 1) + 



f{a-{-t CO) 



(a, 0) +— (a 



(a, 0) H- ^ (a — -i- w, 1) + 



1 • 2 • «ä 
^ « + 1) 



(a, 2) + 



+ 



(^ + 1) M^ — 1) 
1 • 2- 3 



1 • 2 
1 



(a, 2) 



(a — 1.«,3) + 



oder besser 

f{a-\-t w) 

(11) 



und statt (9) 

f {a-\-t CO) 



(12) 



(a, 0) + «), (a — 4 «>, 1) + (« + D« («, 2) 
+ (^ + 1)3 (a — Y '0, 3) + {t + 2)4 (a, 4) 
+ {t + 2)5 (a — l «, 5) + « + 3)» (a, 6) 
+ *« + 3)7(a — -i-«, 7)+--- 

(a, 0) + (Ol (a + 1 CO, 1) + «)2 (a, 2) 

+ (f + Da (a + 1 w, 3) + « + 1)4 (a, 4) 

+ (i + 2)5 (a + -i- <o, 5) + « + 2)« (OS, 6) 

+ (^ + 3)7 (a + -i-<y, 7)+.--. 



* Bei direkter Ableitung von (11) aus (8) würde die rechte Seite lauten: 
(a + 3 iü, 0) + (« - 3)i (a +-f <o, 1) + (< - 2)^ (a + 3 ^o, 2) +;..- 
Die Umwandlung in (11) geschieht dadurch, dass man die Phase t um drei Einheiten 
vergrössert und statt a überall a — 3 ^o setzt. 
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Die Formeln (10), (11) und (12) sind, wie sich gezeigt hat, blosse 
UmgestaltuBgen der Formel von Lagrange. Die Interpolation nach 
diesen Formeln kommt also, wenn man jedesmal bis zur n**"* Dif- 
ferenz geht, darauf hinaus, dass man an die Stelle der vorgelegten 
Funktion y =z F {x) eine gewisse ganze rationale Funktion n^ 
Qracjes tj ^^ (x) setzt, für welche die in (2) angegebenen Bezior 
hungen bestehen. Die rechnerische Benutzung der obigen Formeln 
besteht im wesentlichen darin, dass man die Änderung ermittelt, 
welche die vorgelegte Funktion erfahrt, wenn man von f (a) auf 
das zu interpolierende f {a-\-to)) übergeht. Da es nun vorteilhaft 
ist, lange Rechnungen zu vermeiden, so wird man die Grösse to) 
möglichst klein zu halten suchen. Denagemäss wählt man für a 
stets eines der beiden Tafelargumente, die den Wert x einschliessen 
und berechnet die Phase t aus 

X =^ a-^-tü) oder t = {x — a) : w^ 

so dass t stets als positiver oder negativer Bruch erscheint. 
Dieser Bruch lässt sich sogar zwischen den Grenzen und + 0,5 
halten, wenn von den beiden x einschhessenden Argumenten das- 
jenige gewählt wird, welches a am nächsten liegt.* 

§ 19. Bevor wir mit der Untersuchung der drei Formehi 
(10), (11) und (12) fortfahren, wird es zweckmässig sein, zum 
bessern Verständnis der weiterhin zu behandelnden Fragen die 
in jenen zur Verwendung kommenden Differenzen im Zusammen- 
hang mit den übrigen einer nähern Betrachtung zu unterziehen. 
Zu diesem Zweck bedienen wir uns des nachfolgenden Differenzen- 
Schemas : 



Vergl. die Beispiele des III. Abschnittes; 
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Ein Blick auf diese Tabelle lehrt, dass die einzelnen Glieder des 
Dififörenzenschemas durch Symbole von der Form (a -\- p w, q) 
bezeichnet sind. Yon den beiden Bestandteilen a -\- p w und q 
des Symbols ist der zweite, nämlich die Zahl q, innerhalb einer 
Spalte konstant und stimmt mit der Ordnungsnummer der an- 
gesetzten Diflferenz tiberein. Der erste Bestandteil soll Zeilen- 
nummer heissen. Er ist, wie die Tabelle erkennen lässt, längs 
einer Zeile konstant. Was speziell die Zahl p anbetrifft, so ist 
sie ein gerades oder ungerades Vielfaches von -|, je nachdem q 
gerade oder ungerade ist. Daraus folgt, dass der Ausdruck 2p — q 
stets eine gerade, positive oder negative Zahl ist. So ist z. B. in 

(a + 1 <o, 3) p = -j, q = S und 2 p — g = 5 — 3 = 2. 

Das Verhalten des Ausdrucks 2 p — q ist für die obige Tabelle 
charakteristisch; denn es ist notwendige und zugleich hinreichende 
Bedingung daftir, dass irgend eine Diflferenz (a -\- p w^ q) dem 
Diflferenzenschema der obigen Tabelle angehört. Wenden wii* 
dieses Kriterium auf die im vorigen Paragraphen abgeleiteten 
Formeln an, so erkennen wir, dass in der Tat alle dort in 
Rechnung gezogenen Diflferenzen jene Bedingung erfüllen. Zur 
bessern Übersicht sind sie in der Tabelle durch fettere Schrift 
kenntlich gemacht. 

§ 20. Zurückkehrend zu den in § 18 abgeleiteten Inter- 
poiationsformeln, wollen wir nun zeigen, in welcher Weise die- 
selben zur Herleitung einiger anderer, ebenfalls im Gebrauche 
stehender Formeln dienen können. 

Wir denken uns zunächst einmal das Diflferenzenschema der 
Tabelle XVII um die durch (a, 0), (a, 2), ... gehende Horizontale 
umgeklappt. Dann ändern nach § 8 die Diflferenzen ungerader 
Ordnung ihre Vorzeichen und ausserdem vertauscht sich 

{a-\-p w, q) mit (a — p w^ g), 

so dass z. B. (10) übergeht in 

f(a — tw) = (a, 0) - (Ol (a — T ^^ D + (^)2 i^L — co, 2) 

— (^^3 iß. — I io, 3) + (t\ (a — 2 io, 4) 

— («5 (ö^ — I ^^, 5) + ©6 (a — 3 CO, 6) 



m 
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Ersetzt man hierin t durch — t und berttcksichtigt die Beziehung 

(— t)p = i—iy{t-\-p — Di 
so findet man 

ria-\-t w) = (a, 0) + (Ol (a—^w, l) + «-|-l), (a — w,2) 

4- (f + 2)8 (o — f«j 3) + (/ + 3)4 (a — 2 «, 4) 
+ « + 4)5 (a-|.o,5) + (i + 5)6(a — 3«,6) 

+ « + 6)T(a — 1«>,7) + ---. 

Ersetzen wir ferner in der Formel (11) a durch a -\- lo und < 
durch t — 1, so verwandelt sie sich in 

f(a + t<o)= (a + CO, 0) -{-{t — Dl (a + j .«>, 1) + (0* (a + w, 2) 

+ (Os (« + -!■ '". 3) + « + 1)4 (o + «>, 4) 
+ « + Ds (a + -i- «>, &) + (< + 2)« (a +.<.;, 6) 
+ « + 2), (a + |w, 7)+--.. 

Bildet man jetzt hieraus und aus (12) das arithmetische Mittel 
und setzt zur Abkürzung 

(g, g) + (a + oj, q) ,„ , i . „. 

SO ensteht die Formel 

(/•(« + «.<>) = (a + l«,0)+(02(a+Y«, 2) + (< + 1)4 (a+|«>, 4) 

+ (^ + 2)6(a + 7<«,6) + --- 

+ «—4-) [(ö+Y"'. 1)+tW« («+T^''. 3) 
+ -r(« + l)4(a + T«>,5) 

+^(^ + 2)6(a + l«,7)+-.-]. 

Die hierin vorkommenden Grössen 

(a + 4- «, 0), (a + i- <o, 2), (a + 4 w, 4), (a + ^ w, 6), • • - 

sind, wie oben angedeutet wurde, arithmetische Mittel aus zwei 
aufeinander folgenden Differenzen derselben Ordnung der Tabelle 
XVn. Sie bilden somit nach § 10 selbst wieder ein Differenzen- 
schema, dessen Gheder wie dort in unzweideutiger Weise durch 

* Es ist 

, ,, -t(-t-l){-t-2)... {-t-p + l) , ,^, <«+ l)« + 2)...(<+p-l) 
<-')"= 1.2-3 p == ("^> 1-2.3 ... p 

= (-1/ « I-2J-1V. 
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Zeilen- und Ordnungsnummer vermittelst des Symbols (a-\-pio,q) 
bezeichnet werden können, nur dass jetzt die Differenz 2 p — q 
stets eine ungerade^ positive oder negative Zahl ist. 

§ 21. Was die Anwendung der einen oder andern der 
Formeln (10) bis (14) betrifft, so lässt sich darüber etwa fol- 
gendes festlegen: 

Die Benutzung der Formel (10) ist dann geboten, wenn man 
die Werte {a — lo, 0), (a — 2 oj, 0), ... nicht zur Hand hat, d. h. 
bei Interpolationen in der Nähe des Anfangs einer mathematischen 
Tabelle. Aus demselben Grunde empfiehlt sich die Formel (13) 
zu Interpolationen in der Nähe des Schlusses. Kennt man aber 
alle in der Nähe einer Zeile befindlichen Differenzenwerte, was 
in der Mitte einer Tabelle stets der Fall ist, so zieht man vor, 
die Interpolation mit Hilfe der Formeln (11), (12) oder (14) aus- 
zuführen. Dabei steht dem Rechner die Wahl zwischen diesen 
drei Formeln vollständig frei ; denn jede derselben arbeitet völlig 
sicher, sobald höhere als die in Rechnung gezogenen Differenzen 
auf das Resultat keinen Einfluss mehr auszuüben vermögen. 
Trotzdem wird sich aber einem Rechner, der mit etwelcher 
Überlegung arbeitet, die Frage aufdrängen, welche Formel von 
den ihm zur Verfügung stehenden zu bevorzugen sei. Die nahe- 
üegende Antwort lautet, dass im allgemeinen die Formel mit 
den kleineren Koeffizienten die bequemere sein wird. Es wird 
sich also jetzt darum handeln, eine Vorstellung von dem Ver- 
halten der Koeffizienten zu gewinnen. Zu diesem Zwecke sind 
nachstehend die Koeffizienten zweiter bis siebenter Ordnung ta- 
buliert worden, und zwar für die Zehntel der Phase zwischen 
den Grenzen -(- 1 bis — 1, unter Fortlassung der Phase t = 0,0, 
da für diese die betrachteten Koeffizienten sämtlich verschwinden. 
Die Überschrift gibt die Ordnungsnummer der Koefläzienten an, 
darunter folgen spaltenweise die Phase und die in Einheiten der 
vierten Dezimale angesetzten Koeffizienten für die einzelnen For- 
meln (11), (12) und (14). Da es hier nur auf die absoluten Be- 
träge ankommt, so sind die Vorzeichen fortgelassen. Dafür sind 
den Formelnummern die Zeichen -f- oder — angehängt, welche 
angeben, dass die Phase positiv oder negativ genommen worden ist.* 



Vgl. Bruns. loc. cit. 
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2. 


Ordnung. 


. S 




1 

t 


(11) + 


(11)- 


(12) + 


(12)- 


(14)+ 


(14) 


0,1 


550 


450 


450 


550 


450 


550 


0,2 


1200 


800 


800 


1200 


800 


1200 


0,3 


1950 


1050 


1050 


1950 


1050 


1950 


0,4 


2800 


1200 


1200 


2800 


1200 


2800 


0,5 


3750 


1250 


1250 


3750 


1250 


3750 


0,6 


4800 


1200 


1200 


4800 


1200 


4800 


.0,7 


5950 


1050 


1050 


5950 


1050 


5950 


0,8 


7200 


800 


800 


7200 


800 


7200 


0,9 


8550 


450 


450 


8550 


450 


8550 


1,0 


10000 








10000 





10000 







3- 


Ordnung. 






t 


(11) + 


(11)- 


(12) + 


(12)- 


(14) + 


(14) 


0,1 


165 


165 


165 


165 


60 


110 


0,2 


320 


320 


320 


320 


80 


280 


0,3 


455 


455 


455 


455 


70 


520 


0,4 


560 


560 


560 


560 


40 


840 


0,5 


625 


625 


625 


625 





1250 


0,6 


640 


640 


640 


640 


40 


1760 


0,7 


595 


595 


595 


595 


70 


2380 


0,8 


480 


480 


480 


480 


80 


3120 


0,9 


285 


285 


285 


285 


60 


3990 


1,0 

















5000 
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4. 


Ordnung. 






t 


(11) + 


(11) 


(12) + 


(12)- 


(14) + 


(14) 


0,1 


87 


78 


78 


87 


78 


87 


0,2 


176 


144 


144 


176 


144 


176 


0,8 


262 


193 


198 


262 


193 


262 


0,4 


886 


224 


224 


336 


224 


336 


0,5 


391 


234 


234 


391 


234 


391 


0,6 


416 


224 


224 


416 


224 


416 


0,7 


402 


193 


193 


402 


193 


402 


0,8 


336 


144 


144 


336 


144 


336 


0,9 


207 


78 


78 


207 


78 


207 


1,0 

























S- 


Ordnung. 






t 


(11) + 


(11) 


(12)+ 


(12) 


(U)+ 


(14) 


0,1 


33 


33 


33 


33 


6 


10 


0,2 


63 


63 


63 


63 


9 


25 


0,3 


89 


89 


89 


89 


8 


42 


0,4 


108 


108 


108 


108 


4 


60 


0,5 


117 


117 


117 


117 





78 


0,6 


116 


116 


116 


116 


4 


92 


0,7 


104 


104 


104 


104 


8 


96 


0,8 


81 


81 


81 


81 


9 


87 


0,9 


45 


45 


45 


45 


6 


58 


1,0 














■ 






* - - " J 

d 4 ■*•*■* 
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6. Ordnung. 



t 


(11)+ 


(11) 


(12)+ 


(12)- 


(U) + 


(14)- 


0,1 


17 


16 


16 


17 


16 


17 


0,2 


34 


30 


30 


34 


30 


34 


0,3 


49 


40 


40 


49 


40 


49 


0,4 


61 


47 


47 


61 


47 


61 


0,5 


68 


49 


49 


68 


49 


68 


0,6 


70 


47 


47 


70 


47 


70 


0,7 


64 


40 


40 


64 


40 


64 


0,8 


51 


30 


30 


51 


30 


51 


0,9 


30 


16 


16 


30 


16 


30 


1,0 

_ .... 

























7. 


Ordnung. 






t 


(11) + 


(11) 


(12) + 


(12) 


(14) + 


(14) 


0,1 


7- 


7 


7 


7 


1 


1 


0,2 


14 


14 


14 


14 




3 


0,3 


20 


20 


20 


20 




6 


0,4 


23 


23 


23 


23 




8 


0,5 


24 


24 


24 


24 





10 


0,6 


24 


24 


24 


24 




11 


0,7 


21 


21 


21 


21 




11 


0,8 


16 


16 


16 


16 




9 


0,9 


9 


9 


9 


9 




6 


1,0 
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An der Hand der vorstehenden Zahlen lassen sich jetzt 
leicht bestimmte Aussagen über die Zweckmässigkeit der in 
Betracht kommenden Formeln machen. 

Vergleicht man zunächst innerhalb derselben Formel und 
Ordnung die Koeffizienten, die zu entgegengesetzten Phasenwerten 
gehören, so erkennt man sofort, dass sich in (11) die negative 
Phase, in (12) und (14) dagegen die positive als die günstigste 
erweist. Beschränkt man sich jetzt auf positive Phasenwerte 
bei (12) und (14) und auf negative bei (11) und hält für jede 
Ordnung die Zahlen unter (11)_, (12)-|- und (14)-|- gegeneinander, 
so ergibt sich, dass die Koeffizienten der Formel (14) in allen 
Ordnungen, besonders aber in den höhern, die kleinsten Werte 
annehmen. Daraus schliessen wir, dass für Interpolationen in der 
Mitte einer mathematischen Tabelle diese Formel die bequemste ist. 

§ 22. Die Berechnung der Koeffizienten, namentlich bei den 
höhern Differenzen, bildet gewöhnlich den lästigsten Teil der 
Interpolationsarbeit. Aus diesem Grunde sind die Werte der 
Koeffizienten vielfach tabuliert worden. Solche Tafeln, besonders 
für die Formel (14) finden sich in den meisten Tafelsamralungen, 
die für Zwecke der Astronomie oder Geodäsie zusammengestellt 
sind. Hat man keine Hilfstafeln zur Hand, so ist es zweck- 
mässig, statt (14) die Formeln (11) und (12) zu verwenden und 
zwar deshalb, weil diese leicht in andere, der Berechnung besser 
zugängliche umgeformt werden können. Betrachtet man nämlich 
die Koeffizienten dieser beiden Formeln, so erkennt man sofort, 
dass jeder von ihnen Teiler aller folgenden Koeffizienten ist. 
Gestützt auf diese Tatsache bilden wir, mit (11) beginnend, 
zunächst die Grössen 

Ä = t, B = (t-\-l):2, C=(t—l):S, D = (^ + 2) : 4, 
E=(t—2):5,F={t-\-3):6, G = (t — 3):7, 

u. s. w. 

Mit Hilfe derselben geht dann (11) über in 
f{a -f tiü)= {a,0)-\-Ä(a — ^(o, 1) + ^5 (a, 2) 
+ ABC (a — I io, 3) + ÄBCDia, 4) 
+ ÄBCDE (a — ^ w, 5) + ABGDEF (a, 6) 
+ ABCDEFG [ia — ^ w, 7) -\ 
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= (a, 0)-\-A 



{a — Y (0, 1) -\- B 



(110 



(ö, 2) + C{(a — 1 io, 3) + 1> { (o, 4> 

+ E^(a — ^w,5) + F{{a,6) + G{{a—iio,7)-\-...}}M 

oder kürzer 

/•(a + < «) = (a, 0) + ^- 1, I = (a — Y«, 1) +5- H, 

n = (a, 2)+ c-in, ni -= (a — l «>, 3) + D • rv^ 

IV = (a, 4) + j6? . V, V = (a — yw, 5) + F-VI, 

VI = (a, 6) + G • Vn, VII = {a — jw, 7) +---. 

Ebenso erhält man statt der Formel (12), wenn 

A = t, B = (t—1):2, C = « + 1) : 3, D = (t — 2) : i, 

E = {t -\- 2) : 5, F=(jt — 3):6, G = (t -\- 3) : 7, 

u. s. w. 
gesetzt wird, 

f(a-\-t(o) = {a,0)-^Ä-I, I = (o + |-«», l) + 5-n, 

n = (o, 2)+ C-III, in = (a + -i-w, 3) + i)-IV, 

IV = (a, 4) + i7- V, V = (a + l «>, 5) + i?'' VI, 

VI = (a, 6) + G' vn, VII = (a + l«, 7) +---. 

In gleicher Weise lassen sich auch die Formeln (10) und (13) 
behandeln. Setzt man nämhch: 

A = t, B = it—l):2, C=it—2}:S, D = (< — 3) : 4, 

E=(t — i):5, F={t — S):6, G = {t — 6):7, 

u. s. w., 
so geht (10) über in 

f(a-^tco) = ia,0)+Ä'l, I = (a + i«>,l)4-B-II, 



(120 



(lOO 



n = (o + w,2)-|- C-lII, ni = (a + |«>,3) + Z»-IV, 
IV = (a + 2«>,4)+i:- V, V = (a + -|w,5) + J?'- VI, 



VI = (a + 3 w, 6) + G . vn, VII = (o + |«>, 7)+---. 

Bildet man endlich auch noch die Grössen 

A = t, B = it -{- l) : 2, G = (t -\- 2) : S, Z» = (^ + 3) : 4, 
E = {t -\- 4) : 5, F = (t -\- 5) : 6, ff = « + 6) : 7, 



u. s. w. 



und substituiert dieselben in (13), so entsteht die Formel 
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(130 



n = (a — w,2)+ C-Ill, III = (a — |w,3) + i) -rv, 



< 



IV = {a — 2 CO, 4).+ ^. Y, V = {a — Y^o, 5) + F' VI, 

VI = (a — 3ro,6) + (?- VE, VII = (a — 1^0,7)+.... 

Die Berechtigung zu den in (HO, (12'), (lOO und (13') gemachten 
Ansätzen erhellt sofort, wenn man. die Hilfsgrössen I, II, m,..-.. 
wieder eliminiert. 

§ 23. Bei den bis jetzt abgeleiteten Interpolationsformeln 
für aequidistante Intervalle, über deren praktische Verwendbarkeit 
wir uns noch zu verbreiten haben, besteht ohne Ausnahme der 
Nachteil, dass bei der Ausführung mehrerer Interpolationen 
innerhalb eines und desselben Intervalls die konstant bleibenden 
Differenzen immer wieder als Faktoren der veränderiichen Koef- 
fizienten auftreten. 

Das einzige Mittel, diesen Nachteil zu eliminieren, besteht 
in der Absonderung dieser konstanten Differenzen. Zu dem Ende 
entwickeln wir die Formeln. (10), (11)", (12), (13) und (14) nach 
Potenzen der Phase t 

Wir beginnen mit der Formel (10) und erhalten, indem wir 
die Zähler der darin .vorkommenden Binomialkoefflzienten aus- 
multiplizieren: 
i = t, ■ '^'■•: 'v 

t{t — l) {t — • 2) = fl '^ 3i»; -\r2t, . 

< (# — 1) « — 2) (« — S) = t* — 6<» + 11«« — 6i, 
t{t—l){t — 2) {t — 8) (i — 4) = «5 _ 10^ 4- 35i» — 50<* + 24«, 
t {t — 1) {f — 2)(t — 3) {t — 4) (« — 5) = «« — 15<« + 85^ — 225<8 

+ 274«» — 120«, 
t(t—D(t — 2) (« — 3) .(« ~ 4) {t — 5) (t — 6) = f — 21«« ' 

+ 175«» — 735«* 4- 1624«» — 1764«» + 720«, 

«... , ...;. 

Die vollständige Entwicklung von (10) ergibt also : 
/(« + « CO) = (o, 0) + [«] (a + \co, 1) + [|' — ^] (a + «>, 2) 

+ L6 -2 +3J^'' + ^^'''^^ +124-4 + ^-4] ("^ + ^"''^^ 



54 



, r <» t* It* 5f , t-,^ ,5 ;,, 

+ LTäö-T2+24-T2+TJ^«+2*''^> 



+ l72Ö-i8+li4 - 16 + -36Ö— eJ <« + ^ «*' 6) 

+ 1.5040 240^^144 48^^90 20 + tJ ^" ^^ 2 '«'. '> 



+ 

oder nach Potenzen von t geordnet: 

f {a -{- 1 w) ^ {a, 0) 



-\-t 



(a-\-\w, 1) — i(a4-«»,2)4-|(a4-|«»,3)— ^(a+2.«,4) 



+<* 



^1 (a + «., 2) - ^ (a+ I «,, 3) + ii (a + 2 «;, 4) 



^2(a+T«»,5) + 3gQ 



20 



fl , . 3 ... l . . ^ .- . 7 



+*' 



^(a+l«, 3) - -i- (a + 2 <«», 4) + H7 (a + |«,5) 



6 



24 



29 



_JL(a_|_3«,,6) + ^(a + |..,7)- 



16 



90 



+t' 



'^(a + 2«»,4)-^(a + |«>,5) + ^(a4-3«», 6) 
-;^ (a + !<«», 7)+... 



(»'> 



'^*\'m ^" + ^ '"' ^^ - 2^ (« + T «*. 7) + ...} 

+i5^(« + ^'''^)-- 
+ 
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Für die Formel (11) findet man: 
t = t, 

(f-\-l)t(t — l) = t'~ t, 
(t -j- 2) « + 1) < (« — 1) = if* + 2t» — <ä — 2t, 
{t + 2) « -j- 1) < (< — 1) {t — 2) = t<^ — bt» + U, 
(t 4- 3) (f -|- 2) {t-\-l)t(t — 1) (t — 2) = «• + 3<» — 5t* — 15«» 

+ U^ 4- 12^ 
{t + 3) (« 4- 2) {t-\-V)t(t — 1) {t — 2)(t — Z) = «' — 14«» 

+ 49«» — 36«, 



also 



f{ß-\.tco) = (a, 0) + [«] (a — i «, 1) 4- [|+|] («, 2) 

+ [e - d (« - 2 '"' ^> + [24+12-24-T2J ^''' *^ 



"^ L72O "■" 240 144 48 "•" 180 '^ 6OJ ^ ' ' 
~ L5040 360 ' 



5040 360 ' 720 140 



'''-Tiö](«-T-'7) 



+ • 

oder nach Potenzen von « entwickelt: 



+* 



f{ß-{-tw) = (a, 0) 
(a — I «, 1) 4- 1 (a, 2) — ha — \io, 3) — ^ (a, 4) 



6 



12 



l+^^*~^"''^^+^^'*'^^~liÖ^'*~^"''''^~*" 



+««{ i- (a, 2) - ij- (a, 4) 4- ^^ (a, 6) - ■ 



2 

ri 



24 



180 



4-«» 



(a-4«>, 3)+^(a, 4) 



g V" 2 -' -" I 12 '"' *' 24 



(a — I «, 5) 



48 



720 



(11 '0 



•- ■ Z ■VJ^' 



56 



+ ^*{ ^ («, 4) - l]j («, 6) 
+ --•• 



■1 



360 



(a — 4«, 7) 



1 



m 



Hieraus leitet sich die (12) entsprechende Formel wie folgt ab: 
Man denkt sich das Differenzenschema der Tabelle XVn um die 
Zeile (a, 0) umgeklappt (§§ 8 und 20) und t durch — t ersetzt. 
Auf diese Weise erhält man : 

f(a-\-t(o)=^ {a, 0) 
(a + 4 <ü, 1) _ 1 (a, 2) — i- (a+ I <o, 3) + ^ (a, 4) 



-{- — (a + ^co,5)--^ (a, 6) - j^ {a + ho,7) + — 
+ <»[ ^ (a, 2) - ^ (a, 4) + j|^ (a, 6) - -• 



+ *' 



1 11 1 

ia + ^o), 3) — j2 (a, 4) — 2^ (« + T*". 5) + ^ («, 
7.-1 



6 
+ 



6) 



720 



(a + I .«>, 7) — 



+ ^*(^(a,4)-jij(a, 6)+-.} 



+ <* 



+.^« 



f 1 



+ t' 



120 

i 

1 

720 

1 
5040 



(a + 1 «>, 5) — ^jjT- (a, 6) 



240 



^ (a+l «, 7) + 



360 



1 



(a, 6) 



(a + 1 



(B'O 



I '0, 7)---j 



+ -! • 

In gleicher Weise ergibt sich auch die (13) entsprechende Formel, 
nämlich : 
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+ t 



+ t» 



f (a -\- t lo) = . (a, .0) 
(a — ^«,l)+|(a-«,2)+|(a — |w,3)+j{a— 2w,4) 

1 ^7 7 

+ ^(a— 3w,6)+^(a — 1«, 7)+-.- 



860 : 



20 



+ t 



3{ 



^ (a — 1 .0, 3) + j (« — 2 «>, 4) + ^ (a — I «>, 5) 
+ 4(« — 3<'>,6)+|^(a — I«, 7) + --- 



16 



90 



■^(,a — 2w,4:)-^^{a — jco,5)+~{a — Sw, 6) 



+ <' 



24 



12 



144 



+ 48(« — 1«>, 7) + 



+ i» 



11^ 



+'im^' 



3 «>, 6) + 



240 



(a 



|«,7) + --} 



+ ■ 



^ (a — 1«, 7)+-.-j 



5040 



m 



Übergehend zur Umformung der Formel (14), findet man: 

2^—1 ==2^ — 1, 
t {t i- 1) = P — t, 

(t +.1) t (t — 1) (t — 2)=^t* — 2t» — t^-h 2t, 

(2t — 1) (t + 1) t (t — 1) (f. — 2 = 2^0 — 5i* + ö«ä _ ,2t, , . 

(i +.2) (^ + 1) ^ (^ — 1) (t — 2) (t — 3) = ^« — 3<« — 5^* + 15«» 

, +it^—12t, 

(2t —l)(t + 2) (t + l)t{t — 1) (t — 2) (t —3) -- 2P — 7*8 — It" 

+ 35«* — 7t» — 28«2 -f- 12«, 
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Die vollständige Entwicklung von (14) ei^bt also: 

f(a+tw) = ia + ^<o,0)-^ [t-^ («+!«>, 1)H- [|-|] (a+i«',2) 

+ Ll2Ö~i8 + 48-12öJ ^"^ + ^ '"' ^^ 

+ L72Ö~2iÖ~li4 "T"48'+"18Ö~6ÖJ ^" + « '"' ''^ 



"*" U040 1440 1440 '^ 28! 



«s 



288 1440 860 



^+8y(«+^''''^> 



+ 



oder nach Potenzen von t geordnet: 



+ t 



+ <* 



(a+i«,l)-i-(a+i-<«,2)+^(a+|w,8)+^(a+>,4) 

-liö^''+^"''^^~Ä^'^^"''^^ + 8iö^"+^*'''^^ + - 
i(o+j«,2)-j(a:+ir«>,8)-^(a+|-«,4)+^(a+|«,6) 



+ «» 



+ < 



*^ 



^(aH-|-«>,3)— ^(a + |.«,4) + ^(a + i«;, 6) 

1 • 1 1 

2j(a + -i-«,4) — jg(a + y«,5) — ^(a + l«, 6) 






(14'^ 



f 
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Für die Anwendungen ist es zweckmässig, die beiden von f 
freien Glieder cüeser Formel zu einem einzigen Glied zu ver- 
einigen: nach den Untersuchungen in § 10 und mit Rücksicht 
auf die Tabelle XVH ist nämlich 

(a + i ^e>, 0) = i- (a, 0) + ^{a + w, 0) 

= ^ (a, 0) + ^ (a, 0) + I {a + ^iü, 1) 

= {a,0) + ^ia + \w, 1) 

oder 

(a + 1 ..>, 0) — ^ {a + \<o, 1) = (a, 0). 

Hat man also eine grosse Anzahl von Interpolationen inner- 
halb eines und desselben Intervalls auszufahren, so benutzt man 
eine der soeben entwickelten filnf Formeln und berechnet ein für 
allemal die Koeffizienten von t, t^, t^ u. s. w. 





IIL Anwendung der Interpolationsformeln auf die 
Bestimmung der zu gegebenen Argumenten 
gehörigen Funktionswerte und umgekehrt 



§ 24. Die Verwendung der in den §§ 18, 20, 22 und 23 
•des vorigen Abschnittes entwickelten Interpolationsformeln fllr 
sequidistante Intervalle erstreckt sich sowohl auf die Benutzung, 
als auch auf die Herstellung der mathematischen Tabellen. Die 
beiden Anwendungsgßbiete ergänzen sich logischerweise, indem 
-die Herstellung der Tabellen notwendig der Benutzung voraus- 
gehen muss. Wenn in folgendem ein anderer Weg, als der aus 
dieser Tatsache sich ergebende, eingeschlagen wird, so hat dies 
seinen Grund einmal darin, dass ein Rechner weit häufiger in 
die Notwendigkeit versetzt wird, ehie fertige Tafel zu benutzen, 
als selbst eine solche herzustellen. Dazu kommt, dass an den 
Interpolationsformeln, bevor dieselben bei der Herstellung der 
Tabellen verwendet werden können, Umformungen vorgenommen 
werden müssen, die aus pädagogischen Gründen am besten erst 
-dann behandelt werden, wenn die praktische Verwendbarkeit der 
Ausgangsformehl durch eine Reihe von Rechnungsbeispielen dar- 
getan ist. 

Jedermann, der genötigt ist, bei seinen Rechnungen mathe- 
matische Tabellen zu gebrauchen, hat wohl schon die Erfahrung 
gemacht, dass die weitaus grösste Zahl der in Frage kommenden 
Funktionswerte nicht direkt der Tabelle entnommen werden 
kann, sondern erst durch Interpolation aus den darin tabulierten 
gefunden werden muss. 

Betrachten wir den Aufbau der hiezu notwendigen Formeln, 
so lässt sich sofort erkennen, dass der zur Interpolation erfor- 
derhche Zeitaufwand lim so grösser ist, je mehr Differenzen in 
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Rechnung gezogen werden müssen.- Er reduziert sich somit auf 
ein Minimum, wenn alle höhern als die ersten Differenzen ver- 
schwinden. In diesem Falle gehen nämlich sämtliche Formeln 
über in die eine: 

f{a +tco) = (a, 0) + ^ (a + I lo, l)/ 

Die Interpolation besteht dann ganz einfach darin, dass man 
Proportionalteile** der ersten Differenz zu dem Ausgangsgliede 
hinzuaddiert, eine Operation, die schon Ptolomäus bei seiner 
'Sehnentafel in Aussicht nahm und durch Beifügung einer Diffe- 
renzen-Kolonne zu erleichtern suchte. 

In § 7 wurde darauf hingewiesen, dass es theoretisch möglich 
ist, durch Verkleinerung des Tafelintervalls die Differenzen einer 
höhern Ordnung bis zur UnmerkUchkeit herabzudrücken. Wenn 
eine solche Verkleinerung trotzdem nicht in allen uns zur Ver- 
fügung stehenden Tabellen in oben gewünschter Weise durch- 
geführt werden konnte, so hat dies seinen Grund darin, dass 
mit derselben der räumüche Umfang der Tabelle, und damit 
selbstverständlich auch der Herstellungspreis sehr rasch anwächst. 

Nach diesen einleitenden Bemerkungen gehen wir nun über 
zu den Beispielen. 

Beispiel /. Man bestimme log 32^54 unter der Voraussetzung^ 
dass die siebenstelligen Logarithmen der Zahlen SO, 31, 32, 33,. 
34, 35 und 36 bekannt sind. 

Aus der siebenstelligen Logarithmentafel von Vega Seite 2 
entnimmt man: 



* Das Verschwinden der höhern als der ersten Differenzen tritt ein, wenn 
diese letztem sämtlich einander gleich werden, also wenn : 

. . .(a — -f iOj 1) zrz (a — 4- Wj 1) == (a -|- -5- w, 1) = (a -|- -f f?>, 1) . . . 

** Daher die Überschrift P. P. (Partes proportionales) über der Spalte, in 
welcher sich in den meisten Logarithmentafeln die kleinen Täfelchen vorfinden. 
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Tabelle XVIII/ 



X 


Lag X 


I 


II 


III 


IV 


30 


1,477 1213 














142404 








31 


1,491 3617 




4521 










187888 




+ 277 




32 


1,505 1500 




4244 




22 






183689 


4116(6) 


+ 255 


— 23(6) 


38 


1,518 5139 




8989 




25 






129650 




+ 230 




34 


1,531 4789 


125891 


— 8759 


+ 218 


17 


35 


1,544 0680 


122845 


8546 






36 


1,556 8025 











Da es sich hier um eine Interpolation in der Mitte der Tabelle 
handelt, so ist es zweckmassig, eine der drei Formeln (11), (12), 
<14) zu benutzen. Um dem Lernenden Gelegenheit zu geben, 
■die Resultate, wie sie die drei von einander verschiedenen For- 
meln zu liefern im stände sind, mit einander zu vergleichen, 
sollen hier alle zur Verwendung gelangen. 

1. Berechnung nach (11), Man erhält 

f(a\tco) = (a, 0) + {t\ (a — \ w, 1) + {t + 1), {a, 2) 

+ {t^ 1)3 (a — 4 CO, 3) + it + 2), (a, 4). 
Hierin ist, da f«> = 1, 



t 



0,254 



z\i setzen. Nach obiger Tabelle wird ferner 

(a, 0) = + 1,5051500 

(a — I oj, 1) = + 0,0137833 
(a, 2) = — 0,0004244 

(a — I CO, 3) = + 0,0000277 
(a, 4) = — 0,0000022. 



Die Differenzen sind in Einheiten der siebenten Dezimale ausgedrückt. 
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Der oben angesetzte Wert von t ergibt: 



t -\- 1 = 1,254 
t — 1 = — 0,746 
< + 2 = 2,254 



log (Ol = 9,4048387 

log (t + l)j = 9,2021012 
log {t + 1), = 8,5977187« 
log (t + 2)« =- 8,3486126« 

Es wird also 

log 82,254 = 1,5051500 + 0,254 • 0,0137883 

— 0,1593 • 0,0004244 

— 0,0396 • 0,0000277 
+ 0,02 • 0,0000022 

== 1,5051500 + 0,0035022 

— 0,0000675 — 0,0000011 
oder endlich 

log 32,254 = 1,5085836. 

2. Berechnung nach (12). Es wird 
f(a + tw) = (a, 0) + (t\ (a + I w, 1) + ((k (o, 2) 

+ (^ + 1)3 (a + y (o, S) + (,t + Ih (a, 4). 

Hier ist: 

{a, 0) = + 1,5051500 

(a + \ w, 1) = -{- 0,0133639 

(a, 2) = — 0,0004244 

(a + I c«, 3) -- + 0,0000255 

(a, 4) = — 0,0000022 

log (Ol = 9,4048337 
log {f)i = 8,9765425« 

log (t + Ds = 8,5977187„ 
log (t + D* = 8,2377029 

Es wird somit 

log 32,254 = 1,5051500 + 0,254 • 0,0133639 
+ 0,0947 • 0,0004244 

— 0,04 • 0,0000255 

— 0,02 • 0,0000022 
= 1,5051500 + 0,0083944 
+ 0,0000402 — 0,0000010 oder 

log 32,254 = 1,5085836 



t — 1 = — 0,746 
* + 1 =^ + 1,254 
t — 2 = — 1,746 
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3. Berechnung nach (14). Es wird 

f{a + t ai) ^ {a, 0) + (t), (a + Y ">, 1) + (O2 ia-\- ^ o,, 2) 

+ (^— |) I (Oa (a + i ö*, 3) + (i + 1)4 (a + y CO, 4). 



Hierin ist: 



(a, 0) = + 1,5051500 

(a + i «>, 1) = + 0,0183639 

(a + I «, 2) = — 0,0004116(5) 

(a + 1 w, 3) = + 0,0000255 



log (t). 


— 9,4048337 


t — 1 


— 0,746 


log (Oa 

log (t-'^)^ {th 


— 8,9765425« 

— 7,8903563 


' 1 

t + 1 


— — 0,246 

— 1,254 


log (^ + 1)4 


— 8,2377029 






Man erhält also 






log 32,254 


— 1,5051500 - 


f 0,254 • 


0,0133639 



+ 0,09474 • 0,0004116(5) 
4- 0,008 • 0,0000255 
— 0,02 . 0,0000023(5) 
= 1,5051500 + 0,0033944 
+ 0,0000390 + 0,0000002 oder 
log 32,254 = 1,5085836. 

Die Vergleichung (]er nach den drei verschiedenen Formehi 
(11), (12) und (14) erzielten Resultate erzeigt völlige Übereinstim- 
mung. In vorstehendem Beispiel wurde log 32 als Ausgangswert 
gewählt, d. h. die Interpolation fand nach vorwärts statt. Es ist 
selbstverständlich, dass man zu demselben Resultate gelangen 
muss, wenn log 33 als Ausgangswert angenommen und nach 
rückwärts interpoliert wird : 

1. Berechnung nach (11). Da in diesem Falle die Phase 
negativ zu nehmen ist, so lautet die zu verwendende Formel 
wie folgt: 

f{a — taj) = (a, 0) — (Ol (a — I ö>, 1) + {t\ (a, 2) 

— {t + 1)3 (a — 4 ö>, 3) + (^ + 1), (a, 4). 
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Hiefin ist nun: 



t = 0,746* 
(a, 0) = + 1,5185139 
(a — y ö>, 1) = -j- 0,0133639 
(a, 2) = — 0,0003989 
(a — i w, 3) = + 0,0000255 
(a, 4) ' = — 0,0000025 
log (Ol = 9,8727388 t — \ = — 0,254 
log {th = 8,9765425« t -\- 1 = 1,746 

log it -\- l)s = 8,7414654» t — 2 = — 1,254 
log (t + D* = 8,2377029 
Es wird somit 

log 32,254 -- 1,5185139 — 0,746 • 0,0133639 
+ 0,0947 • 0,0003989 
4- 0,055 • 0,0000255 

— 0,02 • 0,0000025 

= 1,5185139 — 0,0099694 
+ 0,0000377 + 0,0000014 oder 
log 32,254 = 1,5085836. 

2. Berechnung nach (12). Es wird 
f{a — tco) = (a, 0) — (Ol (a -{- ^ ü>, 1) -\- (t -\- 1)^ (a, 2) 

— it-\- 1)3 ia + j <o, 3) -\- (t -\- 2)i (a, 4). 



Hierin ist: 



(a, 0) = + 1,5185139 

(a + |- w, 1) = -j- 0,0129650 

(a, 2) = — 0,0003989 

(a + I a>, 3) =- + 0,0000230 

(a, 4) = — 0,0000017 



log (Ol -= 9,8727388 

log {t + Da = 9,8137530 

log (t + Ds = 8,7414654„ 

log it + 2)4 = 8,5781059« 



i + 1 = 1,746 
t — 1 = — 0,254 
^ -f 2 = 2,746 



Es wird somit 

log 32,254 = 1,5185139 — 0,746 • 0,0129650 
— 0,65126 • 0,0003989 

* Das D^ative Vorzeichen ist schon in der Formel berücksichtigt. 
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-f 0,055 • 0,0000230 
+ 0,03«) • 0,0000017 
= 1,5185139 — 0,0096719 
— 0,0002597 + 0,0000013 oder 
log 32,254 = 1,5085836. 

3. Berechnung nach (14). Man erhält: 

f{a — ta))=^ (a, 0) — (Ol (a + i ß>, 1) + (« + 1), (a-\-\ a>, 2) 

— (^ + i) I a + Da (a + I- «, 3) + « + 2)4 (a + |ö>, 4). 

Hierin ist: 

(a, 0) = + 1,5185139 

(a + ^ ö>, 1) = -f 0,0129650 
(a + y ö>, 2) = — 0,0003874 
(a -\- \ü), 3) = + 0,0000230 
(o + I ü>, 4) — — 0,0000021 



log {t\ = 9,8727388 

log {t + Ds = 9,8137530 

log (^ + ^) ^ (^ + 1)2 = 9,4321497 



t -\-l = 1,746 
< + i- = 1,246 
t -\-2 = 2,746 



log (« + 2)* = 8,5781059,. 

Es wird also 

log 32,254 = 1,5185139 — 0,746 • 0,0129650 

— 0,65126 • 0,0003874 

— 0,270 • 0,0000230 
+ 0,03.7, • 0,0000021 

= 1,5185139 — 0,0096719 

— 0,0002523 — 0,0000062 
+ 0,0000001 oder 

log 32,254 = 1,5085836. 

Die Interpolation nach rückwärts ist immer dann zu empfehlen, wenn 

2t > l. 

Beispiel 2. Es sei die Sonnenkoordinaie X, bezogen auf das mitt- 
lere Äquinoktium des Jahresanfanges (1905,0), für 1905 Januar 
für 13'' 22"' 18,5' mittlere Berliner Ortszeit zu bestimmen.* 



* Man zählt in der Astronomie die mittlere Zeit von Mittag zu Mittag von 
bis 24 Stunden. Der Mittag des 0. Janaar ist also in astronomischer Zeitrechnung 
als 0,0 zu bezeichnen. 
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Aus dem Berliner Jahrbuch entnimmt man S. 22: 

Tabelle XIX. 



Januar 


X 


I 


Tl 


TTT 


0,5. 


0,1689530 


+ 86113 






1,0. 


0,1775643 


+ 85977 


136 


6 


1,5. 


0,1861620 


+ 85835 


142 


7 


2,0. 


0,1947455 


+ 85686 


149 


8 


2,5. 


0,2083141 


+ 85529 


157 


6 


3,0. 


0,2118670 


+ 85366 


163 




3,5. 


0,2204036 









Wie im vorigen Beispiel, so sind die Differenzen auch hier 
in Einheiten der siebenten Dezimale ausgedrückt. Der oben 
gegebene Zeitpunkt fällt nach Mitternacht und ist somit nach 
astronomischer Zeitrechnung Januar 0,5571586 (in Bruchteilen 
des Tages).* Da wir in voriiegendem Beispiel am Anfang der 
Tafel zu interpoheren haben, so ist es angezeigt, die Formel (10) 
ZM verwenden. Dieselbe reduziert sich auf 

f (a -\-.taj) = (a, 0) + (Ol ia + |- oj, 1) + (t), {a + ö>, 2) 

+ (^3 (ä + 1 ö>, 3). 

Hierin ist 

<o = 0,5 und t CO r^ 0,0571586, also i^ = 0,0571586^ : 0,5*^ = 0,1148172. 

Da 2 ^ <C 1 ist, so empfiehlt sich die Interpolation nach vorwärts. 

Demgemäss setzen wir: r 

(a, 0) = -f 0,1689530 

(a + I w, 1) == -j- 0,0086113 
(a -\- w 2) = — 0,0000136 



(a + f «, 3) 



0,0000006 



* Zur Verwandlung von Stunden, Minuten und Sekunden in Bruchteile des 
Tages ist am Schlüsse des laufenden Abschnittes die allen praktischen Bedürfnissen 
genügende Tafel beigegeben. 



68 



' 



log (t)x ^ 9,0581116 t — 1 ^ — 0,8856828 
log \t)» -= 8,7043598« < — 2 = — 1,8856828 
log {% = 8,5027071 

Hieraus folgt: 

X = 0,1689530 + 0,11432 • 0,0086113 

+ 0,051 • 0,0000136 

= 0,1689530 + 0,0009844 

+ 0,0000007 oder also 
X = 0,1699881, 

wobei das dritte Glied völlig unmerklich ist. 

Boispiel 3. Man berechne sin 45^ 40' 12" aus: 

sin 400 ^ 0,6427876 
sin 41" := 0,6560590 
sin 420 = 0,6691306 
sin 480 = 0,6819984 
sin 440 = 0,6946584 
sin 450 = 0,7071068 
sin 460 = 0,7193898 

Mit Hilfe dieser Werte bilden wir die 

Tabelle XX. 



X 


sin X 


I 


II 


ITT 


r 

40» 


0,6427876 


132714 






41» 


0,6560590 


130716 


1998 


40 


42» 


0,6691306 


128678 


2038 


—40 


43« 


0,6819984 


126600 


2078 


38 


44» 


0,6946584 


124484 


2116 


38 


45» 


0,7071068 


122330 


2154 




46« 


0,7193398 




• 
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Der zu interpolierende Funktionswert findet sich am Schlüsse 
der vorstehenden Tabelle. Es ist somit angezeigt, die Fonnel 
<13) zu benutzen. Dieselbe lautet: 

f {a -\- t eo) = {a, 0) -{- it)i {a — ^ w, i) + «-|-l), {a — w,2) 

+ (« + 2)3 (a - I a> 3). 
Hierin ist: 

ö> = 10; t(o^ 0,67», also t == 0,67. 

Da hier 2 i = 1,34 > 1 ist, so ist es zweckmassig, den Funktions- 
wert 0,7193398 als Ausgangswert zu wählen und rückwärts zu 
interpolieren. In diesem Falle geht die obige Formel Ober in 

f{a — t(o) = (a, 0) — {t)x (a — i w, 1) + {% (a — <o, 2) 

— (% (a — I «*, 3). 
Hierin wird nun: 

t = 0,33 
(a, 0) = -|- 0,7193398 

(a — I ü>, 1) -= 4- 0,0122330 
(a — co,2) =-: — 0,0002154 

(a — I cu, 3) = — 0,0000038 



t—l= — 0,67 
t — 2 = — 1,67 



Es ist ferner: 

log (t)i ^ 9,5185139 
log (Oj = 9,0435587» 
log (f)3 = 8,7891539 
also 

sin 45» 40' 12" = 0,7193398 — 0,33 • 0,0122330 

+ 0,11055 . 0,0002154 
4- 0,06 • 0,0000038 
= 0,7193398 — 0,0040369 
+ 0,0000237 -j- 0,0000002 oder 
sin 45» 40' 12" = 0,7153268. 

Wählt man den Funktionswert sin 45» als Ausgangswert und 
setzt dementsprechend t = 0,67, so liegt der gesuchte Funktions- 
wert vollständig ausserhalb der in Rechnung zu ziehenden Dif- 
ferenzen und die Bestimmung desselben erfolgt somit durch eine 
sogenannte Extrapolation. 

Beispie/ 4. Man berechne log cos 21" 16' 28,4". Die sieben- 
stellige Logarithmentafel von Vega ei'gibt auf Seite 417: 
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Tabelle XXL 



X 


Log cos X 


I 


21» 16' 0" 


9,9693704 


82 


21» 16' 10" 


9,9693622 






1 


82 


21« 16' 20" 


i 9,9693540 


82 


21» 16' 30" 


; 9,9693458 

i 


82 


21» 16' 40" 


, 9,9693376 

f 
\ 





Da hier schon die ersten Differenzen konstant sind, so redu- 
zieren sich die uns zur Verfügung stehenden Interpolationsformeln 
auf die ersten zwei Glieder, welche för alle gleich lauten. Es ist 

f(a-]-tw) = (a, 0) + (Ol (« + Y ö>, 1). 

Hierin ist 

o) = 10", t w = 3,4" und somit t = 3,4" : 10" = 0,34. 

Es wird also: 

log cos 210 16' 23,4" = 9,9693540 — 0,34 • 0,0000082 

=: 9,9693540 — 0,0000028 
= 9,9693512. 

Wenn, wie in diesem Beispiel, alle hohem als die ersten 
Differenzen verschwinden, so lässt sich mit Leichtigkeit auch 
die umgekehrte Aufgabe lösen, d. h. zu einem gegebenen Funk- 
tionswert das Argument bestimmen. 

Beispiel 5. Es sä log cos x = 9,9693598. Man bestimme 
den zugehörigen Winkel, 

Vergleicht man den hier gegebenen Funktionswert mit deu 
in der Tabelle XXI tabelüerten, so erkennt man sofort, dass 

9,9693540 < 9,9693598 < 9,9693622, 

also auch 

210 16' 20" > o; > 210 i6' 10". 

Da nun aber x = 2lo 16' 10" -f t w ist, so wird 10 > ^ ö> > 0, 
also endlich 1 > ^ > 0. 
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Unsere Aufgabe besteht also in der Bestimmung der Grösse t 
aus der Gleichung 

9,9693598 = 9,9693622 — t • 0,0000082. 
Hieraus ergibt sich 

^ == 24 : 82 = 0,29, also 

X = 210 16' 10'' -|- 2,9" = 210 16' 12,9". 
Nicht so einfach gestaltet sich die Bestimmung des Winkels, wenn 
ausser den ersten Differenzen auch noch höhere in Rechnung 
gezogen werden müssen. Das hiebei einzuschlagende Verfahren 
möge an zwei Beispielen erläutert werden. 

Beispie/ 6. Man bestimme die Zahl, deren Logarithmus 

2,4971871 beträgt. 

Es ist: 

Tabelle XXll. 



X 

1 
1 


Log X 

1 


I 


II 


1 

312 ; 


2,4941546 


18897 




313 

1 


2,4955443 


13853 


44 


814 , 

1 


2,4969296 

1 ' 


13810 


43 


315 


! 2,4983106 


13765 


45 


316 


2,4996871 







Hieraus ergibt sich die Ungleichung: 

2,4969296 < 2,4971371 < 2,4983106, 
also auch 

SU < X < 315. 

Nun ist aber x =^ 314 -{- t co, somit 

< ^ < 1. 

t bestimmt sich nach (11) aus der Gleichung: 

2,4971371 = 2,4969296 + {t\ • 0,0013853 

— (t -\- 1)2 • 0,0000043 

= 2,4969296 + t • 0,0013853 

— t^ . 0,0000022 — t . 0,0000022 
oder 22 t^ — 13831 t + 2075 = 0. 



n 
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Eine Wurzel dieser Gleichung liegt, wie oben gezeigt wurde, 

zwischen und 1. Den genauem Wert derselben bestimmen wir 

durch Annäherung mit Hilfe der Regula falsi. Durch Erraten 

findet man 

0,1 < ^ < 0,2. 

Bezeichnet man diese beiden t einschliessenden Werte mit U und U 
und ihre Funktionswerte bezw. mit / (<i) und f (U\ also dass 

t^ = 0,1 ; /* (^i) = + 692,12 

t^ = 0,2; f iU) = — 690,12 
und setzt diese Werte in der Formel 



ein, so wird 



t = k —f (ti) . 

^ = 0,1 + 692,12 
= 0,1 + 0,06 



r(U) — fiti) 
1382,24 



= 0,15. 
Es ist somit 314,15 die gesuchte Zahl. 

Wegen der Kleinheit der zweiten Differenzen lässt sich das 
t mit grosser Annäherung auch aus der Gleichung 

2,4971371 = 2,4969296 + t • 0,0013853 oder 
13853 t = 2075 zu t = 2075 : 13853 = 0,1498 bestimmen. 

Beispiel 7. Es sei gegeben sin x = 0,660919017. Wie gross 
ist der zugehörige Winkel? 

Zur Lösung dieser Aufgabe benutzen wir folgende Tafel: 

Tabelle XXIIL 



X 


sin X 


I 


TT 


Tu 


390 


1 

0,629320391 


13467219 






40« 


0,642787610 


13271419 


195800 


4042 


41« 


0,656059029 




~ 199842 








13071577 


301833 


— 8981 


42» 


0,669130606 


12867754 


203823 




43« 


0,681998360 
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Es ist 0,656059029 < 0,660919017 < 0,669130606, 
also auch 41» < x < 42» oder, da w = 1« ist, 

< ^.< 1. 
Nun erhält man nach (14) die Gleichung: 

0,660919017 = 0,656059029 + t • 0,013071577 

— (f)i • 0,000201832 

— (i — h-i- (f)i • 0,000003981 

oder 

0,660919017 = 0,656059029 + t • 0,013071577 

— t^ • 0,000100916 
+ t • 0,000100916 

— t» • 0,000000664 
+ t^ • 0,000000995 

— t • 0,000000332 

oder geordnet: 664 t»-\- 99921 t^ — 13172161 t + 4859988 = 0. 
Wir finden als erste Annäherung 

0,3 < ^ < 0,4, also 
ty = 0,3; f{ti) ^ + 917350,518 
ti = 0,4 ; /• (ti) = — 392846,544. 

Substituiert man diese Werte in der im vorigen Beispiel ange- 
führten Formel, so erhält man 

^ = 0,3 + 917350,518 • ^gj^^ 

oder 

^ = 0,3 + 0,1 . 0,7 = 0,37, somit 
o; = 410 + 0,370 = 41^370 = 410 22' 12". 

Beispie/ 8. Dem „Berliner Astronomischen Jdhrbuch^^ für 1905 
entnehmen wir die Sonnendeklinationen vom 1. bis 5. August für 
den mittleren Berliner Mittag, nämlich: 

189 8' 43^2" 
ir> 53' 36,1" 
17^ 38' 11,5" 
IT 22' 29,7" • 
i?^ 6' 31,0" 

und vergleichen dieselben mit den im „Nautical Almanac^^ für den 
mittleren Mittag in Greenwich notierten: 
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iS« 8' 9,7" 
17° 58' 2,0" 
17" 35' 86,7" 
170 21' 54,8" 
17" 5' 55,0". 

Es soll nun auf Grund der sich ergebenden Verschiedenheiten 
der Längenunterschied Greenwich-Berlin berechnet werden. 

Die Lösung dieser Aufgabe besteht ganz einfach darin, dass 
z. B. irgend ein Deklinationswert des N. A. in das B. J. hinein- 
interpoliert wird.* Mit Hilfe der Angaben des B. J. bilden wir 
folgendes Täfelchen: 

Tabelle XXIV. 



Datum 


Becl. app. d. 


I 


II 


III 


Aug. 1. 


180 8' 43,2" 


15' 7,1" 






„ 2. 


17 53 36,1 




17,5" 








15 24,6 


-17,4 


+ 0,3" 


„ 3. 


17 38 11,5 


15 41,8 


17,2 


+ 0,3 


« 4. 


17 22 29,7 


15 58,7 


16,9 




„ 5. 


17 6 31,0 









Der N. A. gibt nun z. B. für den mittleren Mittag des 2. August 
in Greenwich den Wert 17« 53' 2". Ein Blick auf obige Tabelle 
zeigt, dass 

17» 38' 11,5" < 17« 53' 2" < 17« 53' 36,1. 

Es wird also nach (14) 

170 53' 2" = 17» 53' 36,1" — (^)r • 15' 24,6" — («)» • 17,4" 



oder 



= 17« 53' 36,1" — t -15' 24,6" 
— t* • 8,7" + t • 8,7" 



87 t^ + 9159 t — 341 



0. 



* Man würde selbstverständlich zum gleichen Resultate gelangen , wenn man 
z. B. einen Deklinationswert des B. J. in den N. A. hineininterpoUeren würde. — 
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Hieraus findet man durch successive Annäherung (vgl. Beispiel 6> 

t = 0,0372^ = 53"* 35« 

als Längenunterschied Greenwich-Berlin. Nach den oben benutzten 
astronomischen Jahrbüchern beträgt der genaue Wert 53"* 34,84*. 

Beispiels. Man berechne die scheinbare Bektaszension des Monden 
für 1905 Oktober 15 für 13^ 7^ 12' mittlere Solot hurner Ortszeit. 

1. Berechnung nach dem Berliner Jahrbuch, Im Moment, 
da Solothum 13** 7"* 12* mittlere Zeit hat, zeigt die Sternwarten- 
uhr in Berlin schon 13^ 7"* 12* + 23"* 26,13**, oder 

13»* 30"* 38,13*. 

Um also diese Bektaszension des Mondeis für den genannten 
Zeitpunkt zu finden, interpoliert man im Berliner Jahrbuch den 
jenem entsprechenden Funktionswert. Hiezu benutzen wir die 

Tabelle XXV. 



Datum 
1906 


AR. app. d. C 


I 


II 


ni 


IV 


Okt. 14. 


2'' 3" 2,97» 


22» 


50,75» 








„ 14,5. 


2 25 53,72 






13,90» 




• 






23 


4,65 




2,70» 




„ 15. 


2 48 58,37 


23 


21,25 


16,60 


2,32 


0,38» 


» 15,5. 


3 12 19,62 






18,92 




0,50 






23 


40,17 


19,83 


1,82 


0,50 


„ 16. 


3 35 59,79 


24 


0,91 


20,74 


1,32 


0,50 


„ 16,5. 


4 0,70 


24 


22,97 


22,06 






„ 17. 


4 24 23,67 













Der oben angegebene Zeitpunkt liegt 1** 30"* 38,13* nach Oktober 
15,5. Nun ist (in Teilen des Tages) 13^ 30"* 38,13* = 0,0629413^ 
= t cü, also 

t = 0,0629413^ : 0,5^ == 0,1258826. 



* Berlin (Sternwarte) liegt nämlich 23" 26,13" östlich von Solothum ; daher 
der Unterschied in den Ortszeiten. 



76 



} 



Al3 Interpolationsformel benutzen wir die Formel (14) und erhalten 

AR. = 3" 12"" 19,62« + 0,12588 • 23» 40,17« 
— 0,055 • 19,83« -j- 0,007 • 1,82« 
= 3" 12"" 19,62« 4- 2° 58,77« — 1,08« + 0,01« 
= 3" 15"' 17,32«. 

2. Berechnung nach dem Nauticai Almanac. Da Qreenwich 
30" 8,71« westlich von Solothurn liegt, so zeigt um IS* 7" 12« 
die dortige Uhr erst IS«« 7" 12« — 30» 8,71« ==12'' 37» 3,29«. 
Zur Ausführung der Interpolationen bedienen wir uns der 

Tabelle XXVI. 



Datum 

Okt. 15. 

1906 


AR. app. d. ^ 


I 


n 


9" 
10 


3" 8» 12,28« 
3 10 9,60 


Im 


57,32« 
57,45 


0,13« 


11 


3 12 7,05 




57,58 


0,13 


12 
13 


3 14 4,63 
3 16 2,34 




57,71 


0,13 
0,135 

0,14 


14 


3 18 0,19 




57,85 
57,99 


0,14 


15 


3 19 58,18 









Der Zeitpunkt, für welchen wir interpolieren müssen, liegt 
37» 3,29« nach 12''. Es ist demnach 

^ w = 37» 3,29« = 2223,29«, 

also, da w = V = 3600« ist, 

t ^ 2223,29« : 3600« = 0,61758. 

Man findet somit: 

AR, = 3" 14» 4,63« + 0,61758 • 1» 57,71« 
— 0,118 • 0,135 

= 3" 14» 4,63« + 1» 12,69« — 0,01' 
= 3" 15» 17,31«. 
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Das hier behandelte Beispiel stellt eine der wichtigsten Auf- 
gaben der Tafelinterpolation dar, zeigt es doch, wie man den 
Ephemerideninhalt des Berliner Jahrbuch oder des Nautical AI- 
manac für einen beliebigen Erdort und für einen beliebigen Zeit- 
punkt umzurechnen im stände ist, sobald von ersterem die 
Länge in Bezug auf Berlin oder Greenwich bekannt ist. 

§ 25. Die vorgeführten Beispiele haben wohl zur Genüge 
bewiesen, dass die zur Interpolation herangezogenen Formeln 
sehr genaue Resultate zu liefern im stände sind. Es soll nun 
an zwei weitern Beispielen gezeigt werden, in welcher Weise 
die Formeln (lOO bis (130, die in § 22 entwickelt wurden, ver- 
wendet werden können. Wir beschränken uns dabei auf die 
beiden Formeln (11') und (120, indem wir die eine zum Kontrol- 
lieren des mit Hilfe der andern gefundenen Resultates benutzen. 

Beispiel 10. Man berechne für den 18, Februar 1906 für 7* 
i5"* 22^^ mittlere Zürcher Ortszeit die Deklination des Mars, 

Der Nautical Almanac gibt für den mittleren Mittag in 
Greenwich : 

Tabelle XXVII. 



Datum 
1906 


Dekl. app. d. ^ 


I 


n 


TTT 


Febr. 16. 


3» 8' 52,8" 










„ 17. 


3 27 15,6 


18' 
18 


22,8" 
19,9 


■ 2,9" 


0,2" 


„ 18. 


8 45 35,5 






-3,1 








18 


16,8 




0,2 


„ 19 


4 3 52,3 






3,3 








18 


13,5 




-0,1 


„ 20. 


4 22 5,8 


18 


10,1 


-3,4 


0,1 


„ 21. 


4 40 15,9 


18 


6,6 


3,5 




„ 22. 
' 


4 58 22,5 











Nun beträgt die Längendifferenz Greenwich -Zürich 34"" 12,29* 
(östlich). Der in vorliegende Tabelle hineinzuinterpolierende Zeit- 
punkt liegt also 
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7* 15- 22« — 34- 12,29' = 6" 41- 9,71* 

nach Februar 18,0. Es ist somit: 

< = 6^ 41- 9,71« = 0,2785846^. 

Die Rechnung liefert nun, wenn man mit der Formel (HO beginnt, 
nachstehende Zahlen: 



A 






+ 0,2785846 


B 






+ 0,639 


C 






0,24 


m 






0,2" 


c- ni 






0,0" 


(a, 2) 






3.1" 


n 






3,1" 


BU 






2,0" 


{a — ^w, 1) 




+ 


18' 19,9" 


I 




+ 


18' 17,9" 


A-I 






5' 5,8" 


(a, 0) 


+ 


3« 


45' 35,5" • 


fia-[-tw) 


+ 


3» 


50' 41,3" 



Anderseits erhält man nach (12'): 



J 


+ 0,2785846 


B 


0,361 


C 


+ 0,43 


in 


0,2" 


c-ni 


0,1" 


(a, 2) 


3,1" 


n 


3,2" 


B-n 


+ 1,1" 


(a + y ö>, 1) 


+ 18' 16,8" 


I 


+ 18' 17,9" 


Ä-I 


+ 5' 5,8" 


(a, 0) 


+ 3» 45' 35,5" 


r(a-\-toi) 


+ 3» 50' 41,3" 
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Beispiel 11. Es soll die Rektaszension des Mondes für den 
16. Oktober 1905 Vormittags 10^ 16"^ 12,5' M. E. Z. berechnet werden. 

1. Berechnung nach dem Berliner Jahrbuch. Der Zeitpunkt, 
für welchen in das B. J. hineininterpoliert werden muss, liegt 
nach astronomischer Zeitrechnung 22»» 16"» 12,5« — 6"» 25,16» 
= 22»» 9"» 47,34» nach Oktober 15,0 oder 10»» 9"» 47,34" nach 
Oktober 15,5 mittlere BerUner Ortszeit.* Es ist somit 

t'w = 10»» 9"» 47,34" = 0,4234646^ 

oder, da das Tafelintervall (co) 0,5«» beträgt, 

t = 0,4234646^ : 0,5^ = 0,8469292. 

Da nun 2 ^ > 1 ist, so ist es zweckmässig, nach rückwärts zu 
interpolieren, d. h. die dem 16,0. Oktober entsprechende Rekta- 
szension als Ausgangswert zu wählen und 

t = — 0,1530798 

zu setzen. Unter Benutzung der dem B. J. für 1905 entnommenen 
Tabelle XXV erhalten wir, wiederum mit (11') beginnend, nach- 
folgende Zahlen: 



« 

A 
B 
C 
D 
IV 




0,1530708 
— 0,4235 

0,384 
+ 0,46 

0,50' 


i)-IV 




0,23' 


{a — ^w, 3) 




— 1,82' 


ITT 

c- in 




+ 1,59' 
0,61' 


ia, 2) 

n 




-7 


- 20,74' 

- 20,13' 


^•n 




— 


- 8,52' 


{a — ^co, 1) 




23"' 40,17' 


I 
A -1 

(a, 0) 
f(a — t (o) 


4- 3" 
3" 


23'° 48.69« 

3-" 38,69' 

35-" 59,79' 

32'° 21,10' 



* Der mittlere Meridian der Zone, für welche die M. E. Z. (mitteleuropäische 
Zeit) gilt, geht durch Stai-gard in Pommern, d. i. ein Ort, dessen geographische Länge 
in Bezug auf Berlin 6'" 25,16* (östlich) beträgt. Demzufolge geht die M. E. Z. der- 
jenigen von Berlin 6"* 25,16" voraus. 
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Andererseits erhält man nach (120: 



Ä 


0,1530708 


B 


0,5765 


C 


0,282 


D 


0,54 


IV 


0,50' 


D-IV 


+ 0,27« 


{a-\-^w, 3) 


+ 1,32« • 


in 


+ 1,59« 
— 0,45' 


c-in 


(a, 2) 


+ 20,74' 
— 21,19' 


11 


B-11 


12,22' 


(a + 1 «, 1) 


-f- 24" 0,91' 


I 


+ 23"" 48,69« 


Ä-1 


3"' 38,69« 


(a, 0) 


3" 35"" 59,79« 
3" 32°' 21,10« 


f {a — t (o) 



2. Berechnung nach dem Nautical Älmanac. Die M. E. Z. 
geht derjenigen von Greenwich, welche dem N. A. zu grande 
gelegt ist, um 1*» voraus. Der Zeitpunkt, für welchen wir inter- 
polieren müssen, liegt somit 22»^ 16"* 12,5* — 1^ = 21»» 16"* 12,5«^ 
nach Oktober 15,0. Mit Rücksicht darauf, dass die Mond-Ephe- 
meriden des N. A. das Interwall co =^ l^ = 3600" besitzen, wird 

t = 972,5« : 3600« = 0,270139. 
Zur Ausführung der Interpolation bedienen wir uns des nach- 
stehenden Täfelchens: 

Tabelle XXVllL 



Datum 
Okt. 16. 1906 


AR. app. d. C 


I 


II 


19" 

20 

21 

m 

22 
23 
24 


3" 27- 61,48« 
3 29 50,16 
3 31 48,98 
3 33 47,94 
. 3 35 47,04 
3 37 46,28 


l" 58,68' 
1 58,82 
1 58,96 
1 59,10 
1 59,24 


0,14' 
0,14 
0,14 
0,14 
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Unter Benutzung der Formel (HO erhalten wir: 



Ä 


+ 0,270139 


B 


.+ 0,635 


II 


+ 0,14» 


jB-n 


+ 0,09» 


(a — 1 w, 1) 


+ 1» 58,82' 


I 


+ !■" 58,9P 


A-1 


+ 32,12" 


(ö, 0) 


+ 3" 31- 48,98" 


fia + tco) 


+ 3" 32» 21,10» 



Andererseits wird nach (12'): 



A 


+ 0,270139 


B 


0,365 


II 


+ 0,14» 


B-II 


0,05» 


(a + l CO, 1) 


+ 1™ 58,96» 


I 


+ 1"" 58,91» 


A-l 


+ 32,12» 


(a, 0) . 


+ 3" 31°- 48,98' 


/■ (a + < ü>) 


+ 3'' 32"' 21,10» 



§ 26. Anschliessend an die soeben ausgeführten Interpola- 
tionen mag noch auf eine einfache, leicht zu behaltende Regel 
auftnerksam gemacht werden, die von den Rechnern bei der 
Benutzung der Formel (12')+ mit Vorteil angewendet wird. 
Schreiben wir dieselbe schematisch in der Form: 



(«, 0) + 1 



[(« 



f {a^t w) 
, D] 



I 



II=(a,2) 



«+1 



[(« — y CO, 3)] 



III 



[(a+a>,0)] 



I-(a+|«,l)4- 



t-1 



2 



11 



[(« + ö>, 2)] 



in=(a+i«>,3)+-- 



und fassen die hierin vorkommenden Grössen I, II, III, ... als 
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verbesserte Differenzen und die Produkte t • I, — ^ — II, -~— • III, . . 

als die angebrachten Verbesserungen auf, so ist ohne Mühe erkenn- 
bar, dass alle Differenzen durch die angebrachten Verbesserungen 
den ihjien an der andern Seite des horizontalen Striches gegenüber- 
stehenden Differenzen genähert werden* Die letztern sind im Schema 
in eckige Klammern eingeschlossen. Die grosse Bedeutung dieser 
Regel besteht darin, dass man bei deren Befolgung von der Be- 
rücksichtigung der Vorzeichen vollständig absehen kann. Zur 
Erläuterung diene das Beispiel 10. 

Die bei der Interpolation nach (12') zur Verwendung kom- 
menden Differenzen sind nach Tabelle XXVII: 

(a, 0) = + 30 45' 35,5", 

(a + 4- a>, 1) = + 18' 16,8", 

K 2) ^ = — 3,1", 

(a + I c^, 3) = — 0,2". 

Ordnet man dieselben, ohne weitere Rücksicht auf die Vorzeichen 
zu nehmen, in derselben Weise, wie sie sich in der Tabelle vor- 
finden, und zieht zwischen den sie enthaltenden Zeilen eine 
Horizontale, so erhält man: 



3« 45' 35,5" 


[18' 19,9 "J 


3,1" 


[0,2"] 


"4» 3' 52,3"] 


18' 16,8" 


[3,3"] 


0,2" 



* Um die Richtigkeit dieser Regel zu beweisen, denken wir uns die betrachtete 
Formel mit der Differenz (a -|- -y ö>, 3) abgebrochen, also dass III = {a -\- -^ w^ 3) 
wird. Ausserdem wollen wir zunächst voraussetzen, es sei III >> 0. Dann wird wegen 

0<^^<1 auch < ^:i . III< (a f -i- ce>, 3) oder (a, 2)< (a, 2) + ^^ . III 
<I (a -|- ce>, 2), welche Ungleichung besagt, dass die Differenz (a, 2) durch die Ver- 



besserung 



t+i 



III in der Tat der ihr an der anderen Seite des Striches gegen- 



überstehenden Differenz (a -^- 0)y 2) genähert wird. Ist III •< 0, so ist wegen 
> ^^ . III > (a + 4-ce>, 3) auch (a, 2) > (a, 2) f ^^ • III >(a -f o>, 2). 
Ähnliche Überlegungen, wie die soeben angestellten, führen zu den weitern Un- 
gleichungen (a + -|- a>, 1) ^ (a + ^ o>, 1) -| ^ — • II ^ (a — -f a>, 1) 

und (a, 0) ^ (a, 0) -|- ^ • I ^ (a -\- ö>, 0), wodurch die Richtigkeit der Regel 
offenbar dargetan ist. 
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Die Verbesserung der zweiten Differenz durch die dritte wird nun: 

0,43 • 0,2" = 0,1" 

und ist so anzubringen, dass 3,1" dem 3,3" genähert wird. Folglich 
wird die verbesserte zweite Differenz 3,2". Aus dem gleichen 
Grunde erhält man für die erste: 

18' 16,8" + 0,361 • 3,2" = 18' 17,9". 

Der gesuchte Funktionswert wird jetzt: 

30 45' 35,5" + 0,2785846 • 18' 17,9" = 3» 50' 41,3", 

weil die Verbesserung eine Annäherung an 4^ 3' 52,3" bewirken 
soll. Die hier benutzte sogenannte Strichregel verdient um so 
grössere Beachtung, als ihr Anwendungsgebiet sich nicht bloss auf 
die Formel (12')-|- beschränkt, sondern, wie gezeigt werden soll, 
sich auch auf die Formel (ll')_ ausdehnen lässt* Dieselbe lautet: 

/ (a — ^ o>) = (a, 0) — A . I, I = (a — I o>, 1) — 5 . n, 

n -- (a, 2) — C • in, III = (a — I o>, 3) — D . IV, 

IV = (a, 4) — jE . V, V = (a — I c^, 5) — i^ . VI, 

Vn, VE =(a — |c^, 7) , 



VI = (a, 6) — ö 
worin wie bei (12')_j. 

^ == ^, jB = (^ — 1) : 2, C 

E = {t -]- 2) \ b, F 



(^ + 1) : 3, D = (^ — 2) : 4, 
(^ — 3) : 6, (? = (^ + 3) : 7. 



Um einzusehen, dass die Regel in der Tat auch bei Benutzung 
dieser Fonnel angewendet werden kann, genügt es, dieselbe wie 
<12')-|- schematisch zu schreiben und die dort angestellten Be- 
trachtungen zu wiederholen. Auf das Beispiel 11 angewendet, 
erhält man nach Tabelle XXV (in Beispiel 9) 



[3" 12"' 19,62'] 


23" 40,17- 


[18,92-] 


1,82' 


[0,50'] 


3" 35"' 59,79« 


[24" 0,91'] 


20,74" 


[1,32'] 


0,50- 



Die Verbesserung der dritten Differenz durch die vierte wird 

0,46 • 0,50^ = 0,23« 
und ist so anzubringen, dass 1,82' dem 1,32' genähert wird. 



* In § 21 wurde darauf hiDgewiesen, dass die beiden Formeln (11)- und (12)-j- 
mit zu den empfehlenswertesten gehören. Das gleiche gilt selbstverständlich aucli 
noch für die aus ihnen abgeleiteten Formeln (11')— und (12')-|-. 
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Folglich wird die verbesserte dritte Differenz 

1,82« — 0,23« = 1,59«. 

Hiermit wird die zweite: 

20,74« — 0,384 - 1,59« = 20,13« 

und somit die erste: 

23"' 40,17« + 0,4235 • 20,13« --= 23°^ 48,69«. 

Multipliziert man diese endlich noch mit 0,1530708 und sub- 
trahiert das Produkt von 3^ 35"" 59,79«, so erhält man 

3h 32°^ 21,10« 

als den gesuchten Funktionswert. 

§ 27. Eine zweite Regel, die ebenfalls bei den Formeln (12')+ 
und (11')- empfohlen wird, gründet sich auf die Erscheinung, dass 
die in ihnen auftretenden Correctionsfaktoren A^ B, C, ... beim 

weiteren Fortschreiten gegen den Grenzwert IL ^ konvergieren.* 

Li 

Die Verbesserung der cf'^ Differenz ist also für ein grösseres q 
beinahe gleich der Hälfte der nächsthöheren Differenz. Da nun 
nach der Strichregel die Verbesserung immer so anzubringen ist^ 
dass die verbesserungsbedürftige Differenz der ihr an der andern 
Seite des Striches gegenüberstehenden Differenz genähert wird, 
so ist einleuchtend, dass für den Fall, dass die zur Verwendung 
kommende Formel nur bis zur g^-" Differenz benutzt werden soll^ 
die (g -|- 1)^^ Differenz dadurch angenäherte Berücksichtigung 
findet, dass man als q^^ Differenz statt des deyn Striche zunächst 
liegenden Wertes das Mittel aus den beiden dem Striche anliegenden 
Werten benutzt. 

Wollte man also die im Beispiel 10 zur Verwendung kom- 
mende Formel (120+ mit der Differenz (a, 2) abbrechen, so würde 
die an (a -j- 4- <^j 1) =: 18' 16,8" anzubringende Verbesserung 
0,361 • (3,1" -f 3,3") : 2 oder 0,361 • 3,2" sein. 



* Um dies zu beweisen, betrachten wir die Faktoren der Differenzen gerader 
und ungerader Ordnung gesondert. Der allgemeine Ausdruck der ersteren ist 

^ ~ T t 1 ,. vt n 1 



—■) also linies 1 1 

2' ^^^ «:.-oo Vq 2j 



q q 2 q — oo\~q 2J 2 

Für die letztern findet man ebenso allgemein 






2 

— rr^ + TT, JllSO 
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Anderseits ergäbe sich . für den Fall, dass man iin Beispiel 
11 die Formel (11')- mit (a —4 <^/3) abbrechen würde, als ver- 
besserte dritte Differenz: 

(1,82^ + 1,320 : 2 -- 1,57^ (statt 1,59«). 

Der Unterschied zwischen dem genauen Wert und dem nach 
•der Regel berechneten ist so klein, dass schon die verbesserte 
zweite Differenz mit dem oben gefundenen Werte 20,13' überein- 
stimmt. 

• 

§ 28. Nachstehend geben wir noch einige Beispiele, an 
denen die Verwendung der dritten Forme] gruppe (§ 23) erläutert 
werden soll. 

Beispiel 12. Bei der Berechnung der Sonnenfinsternisse nach 
•der Methode von Hansen* ist es unter anderem notwendig, dass 
man zuerst für einige in der Nähe der wahren Konjunktion (Neu- 
mond) liegende und gleich weit, etwa 0,05 Tag, von einander 
-abstehende Zeiten aus den Mondtafeln die Längen und Breiten 
-des Mondes berechnet oder aus einer Ephemeride durch Inter- 
polation herleitet. Zur Erläuterung des hiebei einzuschlagenden 
Yerfahrens möge nächfolgendes Zahlenbeispiel dienen: 

Am 19. August 1906 findet in den nordlichen Polargegenden 
■eine Sonnenfinsternis statt. Der Nautical Almanac ergibt hiefür 
als mittlere Zeit für die wahre Konjunktion in Rektaszension 

Aug. 19., 12h 33in 47^5. 

oder nach Verwandlung der Stunden, Minuten und Sekunden in 

Teile des Tages 

Aug. 19,5234665. 

Es sollen nun vermittelst Interpolation die Mondlängen für die 
Epochen 



Aug. 


19,42 


Yi 


19,47 


7J 


19,52 


n 


19,67 



» 



. 19,62 



berechnet werden. 



* Vergl. Theorie der Sonnenfinsternisse und verwandten Erscheinungen von 
P. A. Hansen. Abhandlungen der K. Sächsischen Gesellschaft der Wissenschaften IV. 
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Nach dem Nautical Älmanac ist: 

Tabelle XXIX. 



August 
1906 


X a 


I 


TT 


TTT 


IV 


17,5. 


119« 


83' 


55,0" 


6» 


20' 


2,0" 








18. 


125 


58 


57,0 








4' 9,7" 














6 


24 


11,7 




+ 5,4" 




18,5. 

1 

1 


182 


18 


8,7 


6 


28 


26,8 


4 15,1 


+ 1,0 


4,4" 


19. 


188 


46 


35,5 








4 16,1 




4,4 










6 


32 


42,9 




- 3,4 




19,5. 


145 


19 


18,4 








4 12,7 




—3,5 




, 






6 


36 


55,6 


4 9,3<n 


6,9 


8,7 


20. 


151 


56 


14,0 








4 5.8 




—3,9 










6 


41 


1,4 




—10,8 




20,5. 


158 


37 


15,4 








3 55,0 




-4,0 




1 
1 






6 


44 


56,4 




—14,8 




21. 


165 




22 


11,8 


6 


48 


36,6 


3 40,2 






21,5. 


172 

1 

1 


10 


48,4 















Wir benutzen die zu August 19,5 gehörige Länge als Aus- 
gangswert und finden: 



Epoche 


Aug. 19,42 


Äitg. 19,47 


Ätig. 19,52 


Aug. 19,67 


Aug. 19,62 


t w 


0,08 


0,03 


+ 0,02 


+ 0,07 


+ 0,12 


t 


0,16 


— 0,06 


+ 0,04 


+ 0,14 


+ 0,24 



1. Berechnung nach (11"). Zur bessern Übersicht ordnen 
wir die in dieser Formel vorkommenden Klammergrössen wie 
folgt an: In einer 1. Spalte stehen die in Rechnung zu ziehenden 
unveränderten Differenzen der Tabelle XXIX. Hierauf folgen in 
den mit t^ t% t^, t^ überschriebenen Spalten die Zahlenwerte, wie 
sie sich ergeben, wenn man die gegebenen Differenzen mit den 
ihnen vorgesetzten KoeflBzienten multipliziert. 

Es wird also: 
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Differenzen 


t 


t^ 


P 


«* 


+ 6» 32' 42,9" 


+ 6» 32' 42,9" 








+ 4' 12,7" 


+ 2' 6,3«" 


+ 2'6,3»>" 






3,4" 


+ 0,5(n" 




— 0,5m" 




— 3,5" 


+ 0,2«" 


+ 0,1«" 


0,2«" 


— 0,1«" 




+ 6« 34' 50,1..," 


+ 2' 6,5" 


0,8(8) 


0,1»," 



Hierin enthält die letzte Zeile die ausgerechneten KoeflBzienten 
von ^, ^2^ t^ und ^*. Setzt man nun für t der Reihe nach die 
oben gefundenen Grössen, so ergeben sich für die gesuchten 
Mondlängen folgende Werte: 

Aug. 19,42 X(l = 1440 16' 11,2" 

19,47 X^ = 144« 55' 37,4" 

19,52 X^ = 145» 35' 6,2" 

19,57 ^ J =: 146« 14' 37,5" 

19,62 i « = 146« 54' 11,3". 

Die Glieder in t^ und t^ sind völlig unmerklich. 






2. Berechnung nach (12"). Wir erhalten: 



Differenzen 


t 


«ä 


<» 


^* 


+ 6« 36' 55,6" 


-}-6<> 36' 55,6" 








+ 4' 12,7" 


2' 6,3«" 


+ 2' 6,3«" 






6,9" 


+ 1,1»." 




— 1,1«" 




3,5" 


0,2«" 


-f 0,1«" 


-f 0,2«" 


— 0,1«" 




+ 60 34' 50,1(1," 


+ 2' 6,5" 


0,8« 


0,1«" 



Vergleichen wir die hier gefundenen Koeffizienten von t bis t^ 
mit den entsprechenden der Formel (11"), so ergibt sich völlige 
Übereinstimmung. Daraus folgt aber notwendig, dass auch die 
Resultate mit einander übereinstimmen müssen. 



3. Berechnung nach (14"). Die Interpolation nach dieser 
Formel liefert ebenfalls dieselben Klammerausdrücke und somit 
auch die gleichen Resultate. Es ist in der Tat: 
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Differenzen 


t 


t^ 


t» 


«* 


— 6» 36' 55,6" 


+6» 36' 55,6" 








+ 4' 9,2., " 


2' 4,6»," 


+ 2' 4,6*" 


• 




6,9" 


0,5»," 


+ l,?«»" 


. 1,1»>" 




3,7" 


0,2»," 


+ 0,1»," 


4- 0,3«," 


0,1»," 


. 


+ 6» 34' 50,1a," 


4-2' 6,5" 


0,8»," 


0,1»," 



Wenn, wie im vorliegenden Beispiel, die gesuchten Funktions- 
werte äquidistante Intervalle besitzen, so ist es angezeigt, dass 
man die gefundenen Grössen durch Bildung von DifFerenzenreihen 
prüft. Nach den in den §§ 6, 12 und 13 gemachten Ausführungen 
müssen nämlich dieselben einen genügend regelmässigen Gang 
zeigen. Auf unser Beispiel angewendet, findet man: 

1440 16' 11,2;' 



144 55 37,4 

145 35 6,2 

146 14 37,5 



39' 26,2'' 
39 28,8 
39 81,3 



2,6" 

2,5 

2,5 



39 33,8 
146 54 11,3 

Der Gang der Differenzen ist derart regelmässig, dass die Möglich- 
keit einer fehlerhaften Berechnung der Funktions werte völlig aus- 
geschlossen ist. 

Beispie/ 13. Anschliessend an das soeben durchgerechnete Bei- 
spiel sollen die fünf ausgewählten Epochen für wahre Zeit um- 
gerechnet werden,* 

Zur Lösung dieser Aufgabe bedienen wir uns der nach- 
stehenden, A&m. Nautical Älmanac für 1906 entnommenen Tafel: 



* Hansen bezieht die Elemente der Sonnenfinsternis ohne Ausnahme auf die 
wahre Zeit des 1. Meridians (Grreenwich). 
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Tabelle XXX. 



August 
1906 


1 Zeitgl. im 
miUL Mittag 

j - 


I 


n 


17. 
18. 
19. 
20. 
21. 
22. 


4". 4,60". 
3 52,12 
3 39,14 
3 25,67 
3 11,71 

2 57,28 

1 


12,48" 

12,98 

13,47 

13,96 

14,43 


0,50" 
0,49 
0,49 
0,47 



1^ beträgt, findet man für 



Mit Rücksicht darauf, dass o) - 
t nacheinander die Werte: 

0,42, 0,47, 0,52, 0,57 und 0,62 

1. Berechnung nach (11^'). Wie im vorhergehenden Beispiel, 
so bilden wir auch hier das Schema: 



Differenzen 


t 


P 


— 12,98' 


12,98" 




0,49" 


0,24«" 


— 0,24»." 




— 13,22»." 


— 0,24«." 



Setzt man hierin für t der Reihe nach die oben angegebenen 
Grössen, so findet man für die Zeitgleichiing nachstehende Werte : 

Aug. 19,42 Ztgl. = 3- 33,55' 

= 3 32,87 

= 3 32,19 

= 3 31,52 

= 3 30,85. 

2. Berechnung nach (12'0. Man erhält: 



» 


19,47 


»> 


w 


19,52 


» 


w 


19,57 


» 


w 


19,62 


» 



Differenzen 


t 


t^ 


— 13,47« 


13.47« 




0,49« 


+ 0,24(5,« 


— 0,24(5)« 




— 13,22(1^,« 


— 0,24(»« 
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Die Übereinstimmung der Klammerausdrücke mit denjenigen 
der Formel (11") war zu erwarten. Sie beweist die Richtigkeit 
der berechneten Schemas. Um überzeugt zu sein, dass auch die 
gefundenen Werte für. dija Zeitgleichung wirklich die gewünschten 
sind, bilden wir die Differenzenreihen. Wir finden: 

3"» 33,55'' 

~ 0,68^ 
3 32,87 

— 0,68 
3 32,19 

— 0,67 
3 31,52 

— 0,67 
3 30,85 

Die Möglichkeit einer fehlerhaften Berechnung der Funktions- 
werte (Zeitgl.) ist infolge des regelmässigen Ganges der ersten 
Differenzen völlig ausgeschlossen. Die Umwandlung der fünf 
Epochen in wahre Zeit geht nun in der Weise vor sich, dass 
man von jeder die ihr entsprechende Zeitgleichung subtrahiert. 
Man erhält demnach: 

Aug. 19,42 — 3"* 33,55* ^ Aug. 19,10»» 1"^ 14,45- 

Aug. 19,47 — 3"^ 32,87« ■= Aug. 19,11»» 13°» 15,13- 

Aug. 19,52 — 3"» 32,19* -= Aug. 19,12»» 25»»» 15,81" 

Aug. 19,57 — 3»»» 31,52* = Aug. 19,13»» 37»»» 16,48* 

Aug. 19,62 — 3»»» 30,85* = Aug. 19,14»» 49»»» 17,15*. 

Beispiel 14. Über den Planeten Tlietis liegen folgende Beobach- 
tungen vor: 

1903 MittL Zeit AR. 

Petersburg Juni 12,14^ 15"^ 24' 16^ 46"" 55,55' 

Marseüle Juni 12,10^ i« IS' 16^ 47'^ 1,29' 

Algier Juni 12,8^ 48"' 4' 16^ 47"^ 8,77'. 

Es sollen dieselben mit den vorausberechneten Ephemeridemverten 
verglichen werden. 

Als solche benutzen wir die von P. Neugebauer berech- 
neten OppositionS'Ephemeriden, Dieselben ergeben für die Um- 
gebung des 12. Juni nachstehende Werte (vergl. Berliner Jahr- 
buch 1905): 





Tabelle im 


Yerwi 


istioi 

1 


Dezimalen 


Hin. 


; 1 


0,0416667 


1 


2 


0,0833338 


2 '■ 


3 


0,1250000 


3 


4 


0,1666667 


4 


5 


0,2083333 


5 


6 


0,25(K>000 


6 


■ 7 


0,2916667 


7 


8 


0.3333333 


8 


9 


0,3750000 


9 


10 


0,4166667 


10 


; 11 


0,4583333 


11 


12 

t 


0,5000000 


12 


i ^^ 


0,5416667 


13, 


1 14 


0,5833333 


14; 


15 


0,6250000 


15 1 


16 


0,6666667 


16 1 


17 


0,7083333 


17 


18 


0,7500000 


18' 


19 


0,7916667 


19 


; 20 


0,8333333 


20 


21 


0,8750000 


21 


22 


0,9166667 


22 


23 


0,9583333 


23 


24 

1 
1 

1 


1,0000000 


24 

25 
26 
27 

28 
29 


1 
1 
1 




30 
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Tabelle XXXT. 






m M. Z. 


AR. 


I 


TT 


m 


Juni 9. 


16'« 49°' 


27,09« 


— 55,99« 






„ 10. 


16 48 


31,10 


— 55,63 


+ 0,36» 


+ 0,08 


„ 11. 


16 47 


35,47 


55,19 


+ 0,44 


+ 0,08 


„ 12. 


16 46 


40,28 


54,67 


+ 0,52 


+ 0,09 


„ 13. 


16 45 


45,61 


— 54,06 


+ 0,61 


+ 0,07 


„ 14. 


16 44 


51,55 


— 53,38 


-f 0,68 


+ 0,08 


„ 15. 


16 43 


58,17 


52,62 


+ 0,76 




„ 16. 


16 43 


5,55 









Die hier gegebenen Werte für AB. gelten für den mittleren 
Mittag in Berlin. Um eine Vergleichung zwischen BeobacMung 
und Rechnung vornehmen zu können, wird es notwendig sein, 
vorerst die drei Beobachtungszeiten 

Petersburg: Juni 12,14^ 15"* 24" 
Marseille: Juni 12,10»» 1"» 13« 
Algier: Juni 12,8** 48™ 4« 

auf Berliner Zeit zu reduzieren. Nun ist: 

Längenunterschied Berlin-Petersburg = — V" 7" 36,5' 
„ Berlin-Marseille =- + 0»» 32"* 0,27* 

„ Berlin-Algier = -f 0*" 41"* 26,3". 

Die reduzierten Zeiten sind somit: 



Berliner 



Juni 12,13»* 7"* 47,5«' 

■^^^"^^4. { J^m 12,10»* 33"* 13,37«' 

mittl. Zeit T • i c^r^y. r.^™ o/^ o. 

Juni 12,9»* 29"* 30,3' 
oder in Teilen des Tages: 

Juni 12,5470776 
Juni 12,4397369 
Juni 12,3954895. 
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Eine weitere Vorbereitung zur Lösung der gestellten Aufgabe 
besteht in dem Interpolieren der AB, der Tabelle XXXI für die 
soeben gefundenen Zeiten. 

1. Berechnung nach (11"). Um die in dieser Formel vor- 
kommenden Klanunerausdrücke zu finden, bilden wir das Schema: 



Differenzen 


t 


t^ 


t^ 


55,19« 


55,19» 






— 0,52« 


+ 0,26» 


+ 0,26» 




+ 0,08» 


0,01» 




+ 0,01» 




54,94» 


4- 0,26» 


+ 0,01» 



Unter Berücksichtigung, dass der Phase t nacheinander die Werte 

0,5470776 
0,4397369 
0,3954895 

beigelegt werden müssen, erhält man: 



Juni 


12,5470776: 


AR. = 


~ 16" 46"" 10,30» 


Juni 


12,4397369 : 


AB. - 


- 16" 46"' 16,17» 


Juni 


12,3954895: 


AB. = 


= 16" 46°' 18,59». 



Durch Vergleichung der so gefundenen Werte mit dem beobach- 
teten ergibt sich endlich 

für Juni 12,13»» 7"> 47,5«: J AR. = -{- 45,25^ 
für Juni 12,10»» 33°» 13,27«: J AR. = + *542», 
für Juni 12,9»» 29"» 30,3«: J AR. = -\- 45,18», 

wo mit J AB, die Differenz Beobachtung — Bechnung bezeichnet 
wird. 

2. Berechnung nach (12 '0. Das hier in Betracht kommende 
Schema lautet: 



Differenzen 


t 


t^ 


t» 


54,67» 


— 54,67» 






+ 0,52» 


— 0,26» 


-\- 0,26» 




+ 0,09» 


0,01» 




+ 0,91' 




— 54,94» 


+ 0,26» 


-|- 0,01« 
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Die hierdurch gefundenen Klammergrössen zeigen vollständige 
Übereinstimmung mit denjenigen der Formel (11"). Daraus ergibt 
sich aber mit Notwendigkeit, dass auch die Produkte mit t, t^ 
und t^ von den daselbst gefundenen nicht verschieden sein können^ 
oder es müssten sich denn bei der Bildung derselben Rechnungs- 
fehler eingeschlichen haben. 

§ 29. Bei der Besprechung der in Abschnitt II entwickelten 
Interpolationsformeln waren wir genötigt, dieselben, je nach dem 
Orte, an welchem eine Interpolation vorgenommen werden musste, 
in drei Gruppen zu teilen (vergl. § 21). Es erübrigt nun noch^ 
auf einen Fall aufmerksam zu machen, bei dem, entgegen der 
daselbst gegebenen Regeln, die Formeln der zweiten Gruppe, also 
(11), (12) und (14) und ihre Umformungen nicht verwendet werden 
können. Wir denken dabei an die Möglichkeit, dass der zu inter- 
polierende Funktionswert f {a -\~ t oj) nur wenige Intervalle 
von einer irregulären Stelle der tabulierten Funktion entfernt ist. 
Um einzusehen, dass in diesem Falle in der Tat die genannten 
Formeln ihren Dienst versagen, genügt es, darauf hinzuweisen^ 
dass die Wirkung der Irregularität an Ausbreitung zunimmt, je 
weiter man mit der Bildung der Differenzenreihen fortschreitet 
(vergl. die Ausbreitung des Fehlers e in der Tabelle IX, S. 26). 
Sie teilt sich somit den Differenzen der Zeile a oder den anlie- 
genden Zwischenzeilen um so eher mit, je kleiner die Entfernung 
des zu interpolierenden Funktionswertes von der irregulären 
Stelle ist. Um die Interpolation dennoch vornehmen zu können^ 
berücksichtigen wir, dass sowohl die Glieder 

(a, 0), (a + I a>, 1), (a -f- oj, 2), (a + | o;, 3), --., 
als auch 

(a. 0), (a — 4- ^'^, 1), {(^ — <^^, 2), (a — ^ co 3), -•• 
der Wirkung jener In-egularität nicht unterworfen sind, also dass 
die Formeln (10) und (13) mit ihren Umformungen anwendbar 
bleiben. 





IV. Herstellung der mathematischen Tabellen. 



§ 30. Wie aus der Einleitung (§ 2) ersichtlich ist, stellt eine 
mathematische Tafel im allgemeinen den numerischen Verlauf 
«iner Funktion dar. Sie bildet demgemäss eine Zusammenstellung 
von Zahlengrössen — Funktionswerten — , die teils direkt etwa 
aus einem die zu interpolierende Funktion darstellenden allge- 
meinen Ausdruck, teils indirekt durch Interpolation* berechnet 
zu werden pflegen. Was dabei den Zeitaufwand anbetrifft, den 
jede der beiden Berechnungsarten erfordert, so zeigt die Erfahrung, 
dass selbst weitgehende Interpolationen einer ausschliesslich 
direkten Berechnung vorzuziehen sind.** Dementsprechend dient 
•diese letztere Berechnungsart hauptsächlich zur Herstellung von 
Tabellen mit sogenannten groben Intervallen, während das Inter- 
polationsverfahren zur Verfeinerung derselben verwendet wird. 

§ 31. Um den Übergang von dem ursprünglichen groben 
Intervall zu dem feinen Intervall zu bewerkstelligen, bedient 
man sich gerne der sogenannten Interpolation in die Mitte, Die 
hiezu notwendigen Formeln werden aus den im II. Abschnitt 
entwickelten Formeln (10), (11), (12), (13) und (14) dadurch abge- 



* Nach Wolf muss die Ehre, zuerst eine richtige Interpolationsmethode gelehrt 
und bei der Herstellung von Tafeln verwendet zu haben, für Henry Briggs in 
Anspruch genomtuen werden. 

** In „Bruns, wissenschaftliches Rechnen" findet sich hierüber folgende bemerkens- 
werte Äusserung : Wenn z. B. bei der mit dem Intervall h ausgeführten Tabulierung 
von F (x) für den einzelnen Funktionswert durchschnittlich der Zeitaufwand X erfor- 
derlich ist, so verlangt die Tafelstrecke H r= N h den Aufwand N X, Legt man 
dagegen die Tafel zunächst mit dem „groben" Intervall H an und bringt das Inter- 
vall nachher durch N — 1 Interpolationen auf das „feine" Intervall h, so ergibt 
sich, wenn die einzelne Interpolation dur4:;h8chnittlich den Aufwand Y erfordert, für 
die Strecke H der Aufwand X -j- (N—l) F. Der Unterschied ist also (N- 1) (X— Y) 
und zwar zu gunsten der Interpolation, sobald X grösser als Y ist. Dieser Fall ist 
aber die Regel. Die Differenz X — Y ist häufig so bedeutend, dass es vorteilhaft 
ist, für N ziemlich grosse Werte, z. B. 200 oder 300, zu wählen, obgleich bei 
wachsendem N die Anzahl der in der Interpolationsformel mitzunehmenden Glieder 
wächst, so dass Y zunimmt und X — Y abnimmt. 
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leitet, dass man in denselben för t den "Wert -^ substituiert. 
Dadurch erhält man der Reihe nach: 



f{a-\-^w) = (a, 0) + ^ (a + I «>, 1) 



<10"') 



2 

5 7 «1 

~" 128 ^^ + ^ "'' *^ + 256 ^* + ^ "'' ^^ 

21 33 

-1024^"+ ^"'^^ + 2048 



{a-\-laj, 7) 






,(a, 0) + 2 (a - 2 



^"', 1) 



<11'") 



+ I («, 2) - j^ (a - i- «>, 3) 



+ 



128 



1024 



(a, 6) 



256 
5 



2048 



(a — I o», 7) 



<12'") 



/•(« +1 ß.) = (a, 0) + 2 (a + i <", 1) 

- I (a, 2) - ^ (a + 1 «>, 3) 

O Q 



128 



256 



— T7^ («. 6) — 2Ö48 (« + T «'. 7) 



1024 



+ 



<18"') 



r{a-\-\cu) 



^"', 1) 



= («. 0) + 2 (a - 2 

+ I (a - a>, 2) + ^ (a — I o), 3) 



. 35 , 
+ 128 (« 
j 21 . 



2 <", 4) + jS^ (a 



+ 



1024 



3 o>, 6) + 



256 
39 



2048 



— f ö», 5) 

(a — 1 ö>, 7) 
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/^^\* 



(14'") 



na-[-\aS) 



(a-\-^co, 0) 
3 



1 (n _1_ 1 



(a + I ä>, 2) 
5 



+ 19« (^ + T "'' ^) - T7^ ('^ + T *"' ß) 



128 
35 



1024 



+ 



63 



32768 



262144 



(a-\-\ oj, 10) 



Was die Anwendung der einen oder andern dieser fQnf 
Formeln betrifft, so gelten auch hier selbstverständlich die ih 
§ 21 gegebenen Regeln. Es wurde dort unter anderm hervor- 
gehoben, dass sich für die Interpolationen in der Mitte einer 
Tafel die Formel (14)+ ganz besonders eignet. Das gleiche gilt 
somit auch von (14'") und zwar um so mehr, als die Glieder 
mit den Differenzen ungerader Ordnung herausfallen. Diese Tat- 
sache macht es begreiflich, warum bei Interpolationen in die 
Mitte diese Formel mit Vorliebe verwendet wird. Auch wir 
werden uns in den nachfolgenden Betrachtungen derselben, wo 
immer es angeht, bedienen. 

Beispiel 1. Es soll zwischen log 90 und log 100 der nach- 
stehenden Tabelle der Wert für log 95 interpoliert werden. Es ist: 

Tabelle XXXII." 



X 


log X 


I 


11 


III 


IV 

6707 


V 


VI 
- 2376 


80 


1 ,9030900 




- 68395 


■ 




511525 




f 14445 




+ 2564 




90 


1,9542425 




— 53950 




4143 


1117 




1,9771212 


457575 


48799 


+ 10302 


3420 


1-1447 


— 848 


100 

1 


2,0000000 


413927 


43648 


+ 7606 


26% 


+ 869 


• 578 


110 

i 


2,0413927 




36042 




1827 




334 



* Mit Rücksicht auf die grosse Bedeutung, welche bei der Zweiteilung der 
Intervalle dieser Formel zukommt, wurde sie, entgegen der bisherigen Übung, bis 
zur zehnten Differenz (inklusive) entwickelt. Für den Fall, dass auch diese Glieder- 
zahl noch nicht genügen sollte, so ist es ein leichtes, mit Hilfe der zur Formel (16) 
gehörigen Anmerkung (*) weitere Glieder zu bestimmen. 

** Die Zwischen werte (vergl. § 10), deren Berechnung die Benutzung der Formel 
(14) und ihrer Umformungen notwendig macht, sind auch hier wie in den Beispielen 
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Ih (14'") ist nun: 

(a + I cü, 0) = + 1,9771212 

(a + l cü, 2) = — 0,0048799 

(a -|- -j w, 4) = — 0,0008420 

(a + j w, 6) = — 0,0000848, 
also : 

log 95 = 1,9771212 + ^ • 0,0048799 

O 

— ~ ■ 0,0008420 + Y^T • 0,0000848* 

. = 1,9771212 + 0,0006100 
— 0,0000080 + 0,0000004 
= 1,9777236. 

§ 32. Die Tatsache, dass für t =^ -^ alle Koeffizienten der 

Differenzen ungerader Ordnung der Formel (14) verschwinden, 
gestattet eine ähnliche Umformung, wie wir sie in § 22 mit den 
Formeln (10), (11), (12) und (13) vorgenommen haben. Betrachtet 
man nämlich die Koeffizienten 

(O2, {t + 1)4, {t + 2)6, {t + 3)8, {t + 4)io, '-, 

so erkennt man sofort, dass auch hier wie dort jeder von ihnen 
in allen folgenden enthalten ist. Bezeichnen wir die Teiler der 
Reihe nach mit A, B, C, D, E, ...^ so wird 



des dritten Abechnittes durch Fettdruck ausgezeichnet. Die vorzunehmende Ab> 
rundung geschieht auf „grad^, d. h. so, dass z. B. statt 847,5 in der Spalte VI 848 
geschrieben wird. Im fernem sind die Differenzen alle in Einheiten der siebenten 
Dezimale angesetzt worden. 

* Bei der numerischen Berechnung der hierin vorkommenden Produkte ist es 
empfehlenswert, die Koeffizienten der vierten und sechsten Differenz in der Form 

. ^ ^128 ~ 32 128' 1024 ~~ 256 + 1024 • 

zu schreiben. Ailch für die Koeffizienten der achten und zehnten Differenz: 
der Formel (14"^) kann eine ähnliche Umformung von Nutzen sein. Für 

35 V. j.^.. 1 . . 1 . 35 1,3 1,1 1 

32768 ^^"^ °"'° '«»?'«»«''«>»« =^768 = 1^ + 32768 = 1024 + 8192- 32768' 
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/•(« -f i ö>) = (a + 1- ü>, 0) +^ (a + i ö>, 2) + AB (a + I- CO, 4) 

+ ABC (a 4- l CO, 6) + ABGD ia-\-^w, 8) 
+ ÄBCDE (a + ^ ö>, 10) +• ■• 

_(a+|. CO, 0) + A{(a -\-^w, 2) + B{(a + 1 «,, 4) + 

+ C{(a -\- ^ o,, 6) + D{{a -\- ^ w, 8) -{- 

+ i;{(a + je, 10) +.-.}}}}}. 



Zur Bestimmung der A, B, C, D, E, 
Ziehungen : 

t(t — 1)_^ 



dienen nachstehende Be- 



dS) 



C = 



1-2 



f « + 2) it - 3) 



2 



jB 



f « + 1) « — 2) 1 



5-6 



t=^ 



D = 



3-4 

r (f + 3) (< — 4) 1 



7-8 



*=T 



<=i 



E 



f « + 4) (* — 5) ] 



9- 10 

u. s. w., 






welche ausgerechnet ergeben: 



A = 



8 



B 



16 



, = 



5 t. 7 „ 9 

24'^ = -32*^=-4Ö'^-'-^- 



Setzt man diese Werte in der oben gefundenen Formel ein, so 
geht dieselbe über in: 



(18) 



/•(a + |-co)=:(a + |a», 0) — |{(a + |-fl>, 2)-^{(a + |a>,4) 

-^{(« + T-, 10) --..}}}}). 



oder: 



13 5 
* Die in dieser Formel yorkommenden Koeffizienten -ät rsi nJ* "* ianen 

folgendes einfache Bildungsgesets erkennen; Bedeutet allgemein q die Ordnungs- 

,-1 



nnmmer der in Frage kommenden Differenz, so wird der zugehörige Koeffizient 



4? 
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<B) 



/'(a + |c«) = (a4-|ö>,0)— g .1, 



I=(a+ia»,2)-^.n, 



U=-(a + ic.,4)-A.iu,m=(a + |a>,6)-^.IV, 



Unter Benutzung dieser Formel erhält man zur Berechnung von 
log 95 folgende Zahlenwerte: 



• m 


0,0000848 


-in 

24 


— 0,0000177 


(a + i w, 4) 


— 0,0003420 


II 


— 0,0008248 


^ .TT 
16 


0,0000608 


(a + 1 oj, 2) 


— 0,0048799 


I 


— 0,0048191 


r- 


— 0,0006024 


(a + lo*, 0) 


+ 1,9771212 


f(a-\-\a>) 


+ 1,9777236 



§ 33. Bei der Verwendung der Formel (17) kann man sich 
wieder die Kücksicht auf die Vorzeichen ersparen. Fasst man 
»ämlich wiederum die Grössen I, II, III, ... als verbesserte Dif- 



13 5 

ferenzen und die Produkte ^ • I, :r^ • 11, ^rj • III, 

8 lo 24 



als die an- 



gebrachten Verbesserungen auf, so wird im allgemeinen die an 
(d -\- \ o), (^ angebrachte Verbesserung stets so toirken, dass sie 

(a-\-\ CD, q) von der halben Summe av^ (a — o^ q) und (a-\-' 2a)^ g) 



3 5 7 9 

* Bei der nomerischen Berechntmg der Produkte TB'^h ki ' ^I» hö * IV, -77: • V 

lo J4 od 40 

ist es von Natsen, die konstanten Koeffizienten wie folgt zu schreiben: 

A — i->.JL A— 1__L -Z-_i__JL A — i-i. 
16 "~ 4 16' 24~'4 24* 32 ~ 4 32' 40—4 . 40' 



•V- 

• '■ • 
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entfernte Auf unser Beispiel angewendet erhält man als Ver- 
besserung von (a + ^ ö>, 4): 

' ^ . 0,0000848 =: 0,0000177. 

24 

Dieselbe ist, ohne weitere Rücksicht auf das Vorzeichen zu 

nehmen, an (a -f- i ^j 4) so anzubringen, dass 0,0008420 von 

0,0001827 + 0,0006707 ..c^^can *.i> ^ a ^ -a 

-? 1 — i = 0,0004267 entfernt wird. Es wird somit: 

n ^ 0,0008420 — 0,0000177 = 0,0003248. 
In ähnlicher Weise erhält man fttr I: 

0,0048799 — ^ . 0,0008243 -- 0,0048191, 

Ib 

weil die Verbesserung eine Entfernung von 0,0052218 bewirken 
soll. Die an (a -f- y w, 0) anzubringende Verbesserung wird jetzt: 

\ . 0,0048191 -- 0,0006024, 

o 

somit, da 1,9771212 > 1,9722414 ist, 

/•(a + yc^) = 1,9771212 + 0,0006024 = 1,9777236, 

§ 34. Die Anwendung der oben ausgesprochenen Regel wird 
noch wesentlich erleichtert, wenn man nicht erst, wie es in Tabelle 
xyXTT geschehen, die Zwischenwerte ia-\- ~ co^ 0), (a -f- 1- co, 2), 
(a -f- "l" a>, 4), ... berechnet, sondern die DiflTerenzen, aus welchen 
jene gebildet werden, selbst in Rechnung bringt. In diesem Falle 
geht die Formel (16) über in: 



* Der Beweis für die Richtigkeit dieser Regel ergibt sich auf Grund folgender 

Betrachtung: Bedeuten (a — 0;, 9), (a, q\ (a -\- ö>, q) und (a -\- 2 cOf q) vier 

auf einander folgende Differenzen der qten Ordnung (vgl. Tabelle XVII), so besteht; 

zwischen {a -\- ^ cOy q) und (a -\- -y oj, q ] 2) die Beziehung : (a -\- ^ Wf q -\- 2) 

(ß 4- 2 COy q) 4- (a — cOy q) , , i x , . , , 1 . Iv 

BS -5^— L ^-^V~^ ~ (* ^" T ^> ?) o^®^ besser: (a -}- -i- a>, ^ -f- 2) 

f (^4~T^> 9) =^-^ n ' Nun hat andererseits die Verbesse- 

_i 1 
rung stets die Form: (a } -f ^»9) a / 1 9\ (*^ H"T ^j 9 I 2). Hieraus und 

aus der vorhergehenden Beziehung ergibt sich nun ohne weiteres die angeführte RegeL 



• •• 



« 



• •• 
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(18) 






9^ 
40 



{«Xl 



!))l)} 



worin die «o, «2, t«4, - - durch nachstehende Ausdrücke bestimmt sind : 

«0 = (a, 0) + (a 4- ce>, 0), 
CZ2 = (a, 2) + (ö^ + ^j 2), 
«4 = (a, 4) 4- (ö^ + ^) 4), 



Nach obigem wird jetzt die an a, oder (a, 2) + (a -f- <^> 9) 
angebrachte Verbesserung stets so wirken, dass sie die Summe 
(a, q) -\- {a -\- co, q) von der Summe aus der nächstfolgenden und 
nächstvorhergehenden Differenz entfernt. 

Die Formel (18) ist so genau und zugleich so bequem, dass 
man bei der Berechnung einer Tabelle immer gut tun wird, für 
Intervalle, die um eine Potenz von 2 von einander entfernt sind, 
die strengen Werte zu berechnen und dann die dazwischen 
liegenden mit Hilfe der Formel (18) zu interpolieren. 

Beispie/ 2, Es sei: 

Tabelle XXXIII. 



X 


log sin x 


I 


TT 


TM 


IV 


82» 


9,7242097 


233520 


18877 


+ 1927 




84« 


9,7475617 


216570 


16950 


+ 1613 


314 


86« 


9,7692187 


201233 


15337 


+ 1359 


254 


38« 


9,7893420 


187255 


— 13978 


+ 1157 


— 202 


40« 


9,8080675 


174434 


~ 12821 


+ 991 


— 166 


42« 


9,8255109 




— 11830 







Es soll das Intervall derselben durch wiederholte Anwendung 
der Formel (18) auf o) = 30' gebracht werden. 



1 
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Wir beschränken uns auf die Tafelstrecke von 34 bis iO^ 
und finden, indem wir die höhern als die vierten Differenzen 
vernachlässigen* und gemäss der oben gegebenen Regel auf 
die Vorzeichen keine Rücksicht nehmen, zur Berechnung von 
log sin 350 folgende Zahlenwerte: 

a, = 0,0000568 
oa = 0,0032387 
«0 = 19,5167804. 

Die Verbesserung von «2 durch «4 wird 

— ' a^ = 0,0000106 
Ib 

und ist so anzubringen, dass 0,0032287 von der Summe 0,0018877 
+ 0,0013978 = 0,0032855 entfernt wird. Folglich wird der ver- 
besserte Wert von «2 

0,0032287 — 0,0000106 -- 0,0032181. 

Hieraus ergibt sich für die an «o anzubringende Verbesserung 

^ . 0,0032181 = 0,0004023. 

o 

Da nun 9,7242097 + 9,7893420 = 19,5135517 < Oo ist, so muss 
0,0004023 zu «0 addiert werden. Es wird somit der verbesserte 
Wert von «o 

19,5167804 + 0,0004023 -- 19,5171827. 

Nimmt man davon die Hälfte, so ergibt sich 

log sin 350 =- 9,7585913. 

In gleicher Weise erhält man: 

log sin 370 = 9,7794630, 
log sin 390 -= 9,7988718. 

Um die weitere Teilung vornehmen zu können, vereinigen wir 
die hier gefimdenen Werte mit denjenigen der Tabelle XXXin 
und bilden nachstehendes Diflferenzenschema : 



* Die nächsten in Bechnung zu ziehenden Differenzen wären diejenigen der 
sechsten Ordnung. Die absoluten Beträge derselben sind 8 und 16, also so klein^ 
dass ihr Einfluss auf das Endresultat völlig unmerklich ist. 



i 
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Tabelle XXXIT. 






X 


Log sin x 


I 


n 


III 


IV 


34« 


9,7475617 


110296 








35« 


9,7585913 


106274 


— 4022 


+ 191 


• 


36« 


9,7692187 


102443 


3831 


+ 178 


13 


37« 


9,7794630 


98790 


— 3653 


+ 161 


— 17 


38« 


9,7893420 


96298 


— 3492 


+ 151 


10 


39« 


9,7988718 


91957 


— 3341 






40« 


9,8080675 











Die Kleinheit der vierten Differenzen enthebt uns der Notwendig- 
keit, dieselben niitberücksichtigen zu müssen. Die höchsten noch 
in Betracht fallenden Differenzen sind demnach diejenigen der 
zweiten Ordnung. Man hat somit zur Berechnung von log sin 36^ 30' 

«2 = 0,0007484 und 
«0 = 19,5486817. 

Hieraus ergibt sich für die an «o anzubringende Verbesserung 

1 . «2 = 0,0000936. 

o 

Dieselbe ist an «o so anzubringen, dass 19,5486817 von 19,5479333 
entfernt wird. Demnach wird der verbesserte Wert von «o 

19,5486817 + 0,0000936 = 19,5487753, 
Es wird somit 

log sin 360 30' .^ 9,7743876. 
In gleicher Weise findet man auch 

log sin 370 30' = 9,7844471. 

Die Vergleichung der mit Hilfe der Formel (18) gefundenen 
Funktionswerte mit den entsprechenden einer 7-stelligen Loga- 
rithmentafel beweist, dass sich die Teilung des Intervalls in der 
oben besprochenen Weise mit grosser Genauigkeit vollführen lässt. 
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§ 35. Die vorstehenden Betrachtungen über die Zweiteilung 
des Intervalls können noch anders gewendet werden. Berück' 
ßichtigt man nämlich die sich aus der Tabelle XVII ergebenden 
Beziehungen 

(a + I CO, 0) =- 1^ (ö + I «, 1) + (a, 0) 
und 

{a -\- ^ CO, 0) = ia -\- CO, 0) —l (a -\- \ iü, 1)* 
und setzt zur Abküi'zung 



2 



2 (ö + 1 w, 1) = «1, 
g (a -f y CO, 2) = oj, 
(a + 4- CO, 4) = «4, 



128 

u. s. w., 

so erhält man statt (14'") die beiden Formeln 

fia-\-j CD) — (a, 0) = «i — «2 + «4 — 

und " (19) 

(a + ö>, 0) — f (a -\-\ oj) = «1 + «2 — «4 + , 

welche offenbar nichts anderes darstellen, als die beiden neuen 
Differenzen erster Ordnung, in welche die ursprüngliche Differenz 
nach voUführter Zweiteilung zerfällt. Da nun mit den Differenzen 
erster Ordnung zugleich auch die der höheren Ordnungen gegeben 
sind, so kann man, ohne erst die Funktionswerte selbst zu bilden, 
die nächste Halbierung des Intervalls durch blosses Spalten der 
Grössen (19) bewerkstelligen u. s. w., bis man bei der letzten 
Halbierung angelangt ist. Ist diese Arbeit erledigt, so ergeben 
sich die Funktionswerte durch blosses Aufsummieren der berech- 
neten Differenzen erster Ordnung ; man hat nur darauf zu achten, 
dass bei der Spaltung jedesmal die Summe der beiden neuen 
Differenzen genau gleich dem Werte der alten ist. 

Beispie/ 3, Es sei: 



* Bezüglich der Herleitong derselben vergleiche die in § 23 im AnachlnsB an 
die Formel (14") gegebene Anleitung. 
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Tabelle XXXT 


w 

• 




Augtist 
1906 


MonddManz 
Spica 


I 


n 


in 


1,0. 
1,5. 
2,0. 
2,5. 
3,0. 


65« 16' 56" 
72 28 41 
79 40 6 
86 50 19 
93 58 30 


7» 11' 45" 
7 11 25 
7 10 13 
7 8 11 


+ 0' 33" 

— 20 

— 1 12 

— 2 2 
2 52 


53" 

— 52 

— 50 

— 50 



2d 
Es soll dieselbe in eine andere mit dem neuen Intervall co = — 

= 0,125^ umgewandelt werden. 

Um die gewünschte Vierteilung ausführen zu können, ist es 

notwendig, vorerst die in den Formeln (19) vorkommenden Grössen 

«i, «2, «4, ... zu berechnen. Wir finden 

für Aug. 1,25: «i = 3» 35' 52,5", «a = + 0,8" 
für Aug. 1,75: ch = 3» 35' 42,5", aa = — 5,8" 
für Aug. 2,25: «i = 3» 35' 6,5", «a = — 12,1" 
für Aug. 2,75: a^ = 3» 34' 5,5", (h = — 18,4". 

Mit diesen Grössen und den aus ihnen nach (19) berechneten 

neuen Differenzen bilden wir nachstehende Tafel: 







Tabelle XXXTI. 






Aug. 
1906 


1 

1 

I' 


TT' 


in' 


«1 


«2 


1,00. 
1,25. 
1,50. 
1,75. 
2,00. 
2,25. 
2,50. 
2,75. 
3,00. 


3« 35' 51,7" 
3 35 53,3 
3 35 48,3 
3 35 36,7 
3 35 18,6 
3 34 54,4 
3 34 23,9 
3 33 47,1 

i 


+ 1,6" 
5,0 

11,6 

— 18,1 

— 24,2 

— 30,5 

— 36,8 


6,6" 
6,6 

— 6,5 
-6,1 

— 6,3 

— 6,3 


3« 35' 52,5" 
3 35 42,5 
3 35 6,5 
3 34 5,5 


+ 0,8" 
— 5,8 
12,1 
—18,4 
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Die Richtigkeit der berechneten Differenzen, welche den Inhalt 
der mit I' überschriebenen Spalte bilden, erhellt offenbar daraus, 
dass die Summen 

3« 35' 61,7" + 3» 35' 53,3" = 7» 11' 45" 
3» 35' 48,3" + 3« 35' 36,7" ^ T> 11' 25" 
3» 35' 18,6" + 30 34' 54,4" = 7« 10' 13" 
30 34' 23,9" + 3» 33' 47,1" =- 7« 8' 11" 
mit den in der Tabelle XXXV notierten ersten Differenzen über- 
einstimmen. Der weitere Verlauf der Rechnung besteht jetzt 
darin, dass wir die in der Spalte I' des obigen Schemas tabu- 
lierten Differenzen in genau derselben Weise behandeln wie die 
ersten Differenzen der Tabelle XXXV. 

Statt der Formehi (19) bedienen wir uns der folgenden : 

/(a + |o) — (a, 0) = 01'—«,'+ «4' 

und (19') 

f(a+\a))—f(a-\-^w)=. «/ + «/ — «4' +•••• 
Die Tabelle XXXVI liefert nun 

fUr Aug. 1,625: a,' = 1« 47' 54,15", oa' = — 1,04" 
fttr Aug. 1,875: a/ = 1» 47' 48,35", «»' = — 1,86" 
für Aug. 2,125: Oi' = 1» 47' 39,30", «a' == — 2,64" 
für Aug. 2,375: a,' = lo 47' 27,20", a^' = — 3,42". 
Mit diesen Werten bilden wir die 

Tabelle XXXVII. 



August 
1906 



1,000. 
1,500. 
1,625. 
1,750. 
1,875. 
2,000. 
2,125. 
2,250. 
2,375. 
2,500. 
3,000. 



MonddMam: 
Spica 



66« 16' 66" 

72 28 41 

74 16 36 

76 4 29 

77 52 19 
79 40 6 
81 27 48 
83 15 25 

85 2 56 

86 60 19 
93 68 30 



I" 



10 


47' 


55,19" 




47 


53,11 




47 


50,21 




47 


46,49 




47 


41,94 




47 


36,66 




47 


30,62 




47 


23,38 



«l 



1« 47' 54,15" 



1 47 48,35 



1 47 39,30 



1 47 27,20 



«a 



— 1,04" 



— 1,86 



2,64 



3,42 









1 


Argument 

i 1 


1 

• 


4 


a 

9 


^a-/ - '■ 


1 
t 

a — 

t 


4- ' 


1 

i a — 


I-. 


1 

! a — 




1 

a — 


i-; 


a 

1 




1 a — 

1 


1- 


1 

a — 




a — 


3 


a — 


> 


i 


'-,' - 


1 

1 

1 < 


^- 
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Betrachtet man die unter I" verzeichneten DiflPerenzen erster 
Ordnung, so erkennt man, dass auch hier die oben gestellte 
Bedingung, welche dahin lautet, dass bei der Spaltung jedesmal 
die Summe der beiden neuen Differenzen genau gleich sein muss 
dem Wert der alten, erfüllt ist. 

§ 36. Wir haben bis jetzt angenommen, die zu verkleinernden 
Tafelintervalle seien um eine ganze Potenz von 2 von einander 
entfernt. Nachstehend soll jetzt die Methode besprochen werden, 
derer wir uns bei den ebenfalls häufig vorzunehmenden Drei- und 
Fünfteilungen zu bedienen gedenken. 

Wenn wir bei der Interpolation in die Mitte genötigt waren,. 

in den Interpolationsformeln ^ = - zu setzen, so wird in ent- 

sprechender Weise eine Drei- bezw. Fünfteilung offenbar dadurch 
erzielt werden können, dass man in denselben Formeln der Phase t 

12 1 2 S 4 

nacheinander die Werte ^, -, bezw. -, -^, -3-, - beilegt und aus den 

ö ö ü o o o 

so erhaltenen Ausdrücken die gewünschten Funktionswerte un- 
mittelbar berechnet. In diesem Sinne aufgefasst, würde uns die 
Teilung der Intervalle allerdings kaum etwas neues bringen. Daher 
wenden wir uns denn auch ohne weiteres derjenigen Methode zu,, 
die, wie die in § 35 besprochene, darin besteht, dass man statt 
der Funktionswerte ein Differenzenschema berechnet, aus welchem 
jene durch blosses Aufsummieren der Differenzen erster Ordnung 
hergeleitet werden können. 

Als Differenzenschema verwenden wir dasjenige der neben- 
stehenden Tabelle XXXVIII, welche in genau gleicher Weise 
gebildet ist, wie die früher besprochene Tabelle XVII und auch 
in diese dadurch übergeführt werden kann, dass man ä = 1 setzt. 

Im Gegensatz zu dem in § 35 auseinandergesetzten Verfahren 
muss nun aber bemerkt werden, dass von den in der Tabelle 
XXXVIII sich vorfindenden Differenzen nur diejenigen direkt^ 
d. h. mit Hilfe von Formeln, berechnet werden, welche daselbst 
durch fettere Schrift kenntlich gemacht sind. Die Lücken der 
übrigen Glieder werden dadurch ausgefüllt, dass man, mit der 
letzten Differenzenreihe beginnend,* reihenweise durch successivea 

* Ais letzte Differenzenreihe ist selbstverständlich immer diejenige aufzufassen^ 
deren unmittelbar berechnete Glieder derart wenig von einander abweichen, dass die^ 
übrigen leicht daraus abgeleitet werden können. 
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Addieren oder Subtrahieren die Glieder der um eine Ordnung 
niedrigeren Differenzenreihe herleitet und so zur Funktionenreihe 
«elbst ansteigt. Dabei geben die im voraus unmittelbar berech- 
neten Glieder den festen Rahmen, auf welchen die Rechnung 
allseitig abschliessen muss. Es müssen nämlich nachstehende 
Bedingungsgleichungen, in denen die rechten Seiten durch eben 
jene Glieder gebildet sind, genau erfüllt sein: 



k—l 



<20) 



k 

__fc— 1 



2 

i = k 




2i-l 



(a-f Wö>-] — 2^a>, 3) = (fl+/i-f lce>,2)— (fl4-/io>,2) 



fc— 1 



7 , (a + Wö>+^ö>,2) =(fl-f -^±-a>,|) — (fl-| — ^:p-a>,l) 



k 



2 




2i-l 



{a-\-no)-] — — — ö>, 1) = (fl -f- ^ + ^ ^> •) — (fl -f- i» ö>, 0) 

^=...,_3,_2, — l,0, + l, + 2, ... 

Bei der Herstellung der Tabelle XXXVIII gingen wir von der 
Annahme aus, es sei k eine ungerade Zahl Demnach kann der 
Übergang zur Drei-, bezw. Fünfbeilung ganz einfach dadurch 
bewerkstelUgt werden, dass man in der Tabelle sowohl, als auch 
in den obigen Gleichungen ä = 3, bezw. ä = 5 setzt. So erhält 
man z. B. aus (20) für k = S: 



(21) 



(a + nce> + ^ce>,4) =(flH — ^y^^,3) — (a-j — ~^^ö;,3) 
2i-l 





[aA-nco-A — cw, S) = {a-\-n-{-\o), 2) — (fl + iiö>, 2) 

o 
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f < = 4- 1 



(21) 




1 = — 1 

«==3 



/ I 1 ^* ct\ / 1 2/1 4-1 .V . ,2/1 — I ,. 




2i 1 

{a-\-nw-\ — w,l) = (a-r{-n-{-\co,0) — (a-\-noj,iy 



n 



6 



3, - 2, -- 1, 0, + 1, + 2, + 3, 



und fllr k = 5: 



(22) 



i = 4-2 




(a-\-nco-]-^co,4:) =(fl + -^^ö>,3)— (fl-| — ^- — 



f=r— 2 

1 = 5 



t 
5 



w,iy 





i = 5 




2i 1 

{a-\-na)-\ TTT— ö>, 3) = (fl + '» + ^^j 2) — (fl + '^^j 2> 



/ I I ^' o^ / I 2/1+1 -^ / I 2/1 — I ,. 

ia-{-noj-{--co,2) =(flH 2^ö>,l) — (flH — a>,l> 



2 ^' 1 

{a-j-no)-] — — a>, l) = (fl-]-/i-|-lcy, 0) — (fl + /i<w, 0)^ 



/i 



10 



., - 3, - 2, - 1, 0, + 1, + 2, + 3, 



§ 37. Übergehend zur Berechnung der in der Tabelle XXXVIII 
ausgezeichneten Differenzen, wird unsere Aufgabe zunächst wohl 
darin bestehen müssen, dass wir zwischen diesen und den Dif- 
ferenzen der ursprünglichen Tabelle (etwa XVII) der numerischen 
Eechnung leicht zugängliche Beziehungen aufzustellen suchen. 
Zur Herleitung derselben bedienen wir uns der Formel (14),. 
welche lautet: 

/•(a + ^a>)=(a + jö>,0) + (^ — ^)(a + |a>,l)+(02(a+^ö>,2)• 
+ (^)2 • "1 (t — ^) (a + |a>, 3) + (t+1), (a + |o., 4) 
+ (^ + l)4-^(^ — -|)(a + |ö>,5) + (^ + 2)6(a + jö>, 6> 



+ (^ + 2)7.^(a+|a>, 7) + 



110 

Ä — 1 *4-l 

Hierin setzen wir nun nacheinander t = ^ , , bezw. t = ^\ , 

1 1 ~~~" fi% 1 I wi 

t)der fttr Ä = — : ^ = — ?; — , bezw. t = ^ — und erhalten 

m 2 2 

= (a + ia,, 0)+|(a+|-a>, 1) 
+ 8 11^'" + ^' + ^^^"+^'"' ^^ 

l = — 1 

+ S TT ("» + 2 i + 1) (a + I ö,, 3) 



n 



48 

• = -2 

+ 384 IT ^"^ + ^ * + ^^ ^" + ^ *"' *^ 
+5^ TT (»« + 2 » + 1) (a + I «., 5) 



R 



3840 

t=— 3 



+ 46Ö8Ö n (»« + 2 i + 1) (a + I «*, 6) 



m 



2 
i=— 3 



+ '•"' 



-(a + io>, 0)+|(a+|a>, 1) 



^ JT (»» + 2 i + 1) (a + I a>, 2) 



i=— 1 



+ ä JT (»« + 2 t + 1) (a + 4 c«, 3) 



t — — 2 



i = — 2 

+ 0^^ TT (»« 4- 2 i + 1) (a + T *«, 5) 



3840^^ 



n. 
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i = — 3 

(w + 2 i + 1) (a + i <ü, 6) 



1080, IJ 



46080 

2 
« = -8 

+ 57^ TT (»» + 2i 4- 1) (« + T «». 7) 



645120 _ 



u. 



= w (a -}- i ö>, 1) 
1=— 1 

i = -2 

+ tS7^ TT (»» + 2t + 1) (a 4- t «', 5) 



U. 



1920 

i = — 3 



+ 82i6Ö,.IJ;^ + 2*4-l)(« + i«',7) 



+ 

Wir können hier a gegen a -\- n at vertauschen und erhalten 
allgemein : 

(a -I ^ ö>, 1) = m (a 4 ^^ — o>, 1) 



<23) 



+ T^'n/*^ + 2 i + 1) (a + ^-^ -, 5) 



i = — 2 

1920 

i = — 3 

322560 

2 

+ 

oder kürzer: 



+ 5^S^ 11/*^ + 2i + 1) (a + ^^ö^ -. 7) 



3) 



<24) 



/ I 2» + l .X , I 2n+l .. , ., , 2n+l 

+ 

n= 3, _ 2, — 1, 0, + 1, + 2, + 3, - 



mit 



(25) 
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i = — 1 



r = — 3 



^ = 32256Ö,.I|^(^+2^+l^ = 32266ö(*^"^> <^"^> ^^^^-^^^^ 

durch welche Beziehungen die in der Spalte V der Tabelle XXXVIII 
vorkommenden Differenzen 

..., (fl - 1 .., 1), (fl - i a>, 1), (fl + i a>, 1), (fl 4- » a>, 1), ... (26) 

offenbar eindeutig bestimmt sind. Um auch für die Differenzen 

..., (fl - «>, 2), (fl, 2), (fl + w, 3), ... (27) 

der Spalte 11' ähnliche Formeln zu erhalten, bilden wir das Dif- 
ferenzenschema der nebenstehenden Tabelle XXXIX in der Weise^ 
dass wir die nach (24) berechneten Grössen (26) als Funktions- 
werte betrachten. Dieses Schema hat nun für die Grössen (27) 
die gleiche Bedeutung, wie dasjenige der Tabelle XVII ftlr die 
Grössen (26). Nach (24) wird somit: 

(fl, 2) = m { m (a, 2) + ^ (a, 4) + B {a, 6) -| } 

+ Ä[m (a, 4) + -4 (a, 6) H j + B{m(a, 6) -\ }+ --^ 

= m^ (a, 2) + wi ^ (a, 4) + w -B (a, 6) -| 

-^ m Ä(a, 4:) -^ m B (a, 6) -\- "' 

+ ^Ma,6) + ..., 
oder kürzer: 

(fl, 2) = m2 (a, 2) + Ä' (a, 4) + B' (a, 6) -\ . 

Vertauscht man wiederum a gegen a-\-n w^ so wird 

{a-\-na),2)=^mHa + nw,2) + A'ia+no}^4:)+B'{a+noj,6)'{-;{28} 

worin zur Abkürzung 

Ä'=2mÄ und B' =2 7n B-{-Ä^ 

gesetzt wurde. In gleicher Weise findet man unter Benutzung 
der Tabelle XL und der Formel (24) die Beziehung: 

- I 

+ B" (« + ^^ «>, 7) +- \m 




;/ 



3. 






mia — et), 4: 

m (a, 4) -f- -^ 
w (a -(- w, 4 



( :«, 6) 







5. I 


■ 
1 
■ 






i CO, 6) 






> 


m2 (a - 


- V oj, 5) 



m^ (<^ H" 2 ^^' ^^ 



« 



1 
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]^ 1 *+l 

Hierin setzen wir nun nacheinander t = , , bezw. t = ^ . , 

t)der tiXr k = —: t = ~^ , bezw. i = — ^^^ — und erhalten 

w 2 ' 2 



(o + la>, 0)+|(a+|c«, 1) 



+ ^ JT (w + 2 i + 1) (o + 1 ü>, 2) 
i=— 1 



n 



+ 5 I (m + 2 t -H 1) (« + I «>, 3) 



m 

. = -2 



384 

i = — 2 

m 



,IJ, 



+ 3iö I ' (»« + 2i + l)(a + |ö>, 5) 



.=-3 

+ 46080.11/'" + 2i + 1) (a + 1 «,, 6) 



2 
i = — 3 



+ '•"•" 



m 



_(« + !«,, 0)+J(a+|ö;,l) 



— ^ JT (m + 2 i + 1) (a 4- I CO, 2) 



i=— 1 



+ ^ ]!"(»» + 2t + !)(« + i«*. 3) 

i ^ — 2 

~m 17^*" + ^' + ^^^" + ^"''*^ 

i = — 2 

+ 3S0 n/*" + 2i + 1) (a + 1 «;, 5) 



46080.11 
m 
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1 = — 3 

(m + 2 i + 1) (a + i ö>, 6) 



+ 645120.11/^ + 2i + 1) (« + T «>, 7) 



= m (a + i ^7 1) 
i==— 1 

+ 1^ n (»» + 2i + 1) (a + 4 o*, 5) 



,n. 



1920 

i = — 3 

m 



+ öööHä7^ (m + 2t + l)(a + i«;, 7) 



822660 II 



Wir können hier a gegen a -\- n w vertauschen und erhalten 
allgemein : 

(a -j J_ ü>, 1) = m (a -I J— «>, 1) 

+ S'n ^'^ + 2t + 1) (a + ^^«>,3) 



<23) 



+ liö.n/*" -f- 2i + 1) (a + ^^^ «*, 5) 



i = — 3 



+ 32i6Ö.n ^"^ + 2i + 1) (a + 2n^ ,^ ^^ 

+ ' ••' 

■oder kürzer: 



<24) 



/ I Zfl + I A / I 2«-|-l ,, , ,. , 2n+l -, 
(a-j -l_a>,i) = w(a-| J— ß>, 1) + ^ (a -| g^^'iS) 

+ 

n = 3, _ 2, — 1, 0, + 1, + 2, + 3, - 



mit 
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(25) 



l=— 1 






322560.1 I' ^ » ' 322560 



durch welche Beziehungen die in der Spalte F der Tabelle XXXVIII 
vorkommenden Differenzen 

. . . , (a - I c., 1), (a - 1 c., 1), (a + i a>, 1), (a + f a>, 1), . . - (26) 

offenbar eindeutig bestimmt sind. Um auch für die Differenzen 

. . ., (a - c, 2), (a, 2), {a + .e>, 2), ; - - (27> 

der Spalte 11' ähnliche Formeln zu erhalten, bilden wir das Dif- 
ferenzenschema der nebenstehenden Tabelle XXXIX in der Weise^ 
dass wir die nach (24) berechneten Grössen (26) als Funktions- 
werte betrachten. Dieses Schema hat nun flir die Grössen (27) 
die gleiche Bedeutung, wie dasjenige der Tabelle XVII filr die 
Grössen (26). Nach (24) wird somit: 

{a,2) =m{m (a, 2) -^ Ä (a, 4) -f 5 (a, 6) -j } 

+ ^{ m (a, 4) + ^ (a, 6) -j } + B{m{a, 6) + --•}+ •- 

= m^ (a, 2) + w ^ (a, 4) -f- m B {a, 6) -j 

-\-m Ä{a,4:)-{'mB(a,6)-\ 

+ ^Ma,6)+..., 
oder kürzer: 

(fl, 2) = m2 (a, 2) -f- Ä' (a, 4) + B' (a, 6) -| . 

Vertauscht man wiederum a gegen a-^n w, so wird 

{a + no),2)=^mHa + no),2) + Ä'ia+no),i)+B'{a+noj,6)+.-,{28} 

worin zur Abkürzung 

Ä' = 2mÄ und B' = 27n B-^ Ä'^ 

gesetzt wurde. In gleicher Weise findet man unter Benutzung 
der Tabelle XL und der Formel (24) die Beziehung: 

+ £"(« + ^^«',7)+- (29> 



/ 






m{a — o;, 4 

m (a, 4) -]- ^ 
m (a -f- ö;, 4 



o., 6) 



te :«>, 6) 



5. Ü 



m^ (a — ^ o^, 5) 
m^ {a -\- \ oj^ 5) 
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mit: Ä" -^Sm'^ Ä, B'' ~= S m'- B -^ S m Ä^, 
aus welcher nun auch die ausgezeichneten Differenzen der Spalte 
III', nämlich 

-, (a-^oj, 3), (a-^co, 3), (« + i co, 3), (a + | ro, 8), -. (30) 

ohne Mühe berechnet werden können. In der Praxis wird es 
wegen der Kleinheit der Koeffizienten, mit welchen die in den 
rechten Seiten der Formeln (24), (28) und (29) vorkommenden 
Differenzen der Ausgangstabelle (XVII) multipliziert erscheinen, 
nur ganz selten notwendig sein, ausser den in (26), (27) und 
(80) angegebenen auch noch diejenigen der Spalte IV' oder gar 
noch der Spalte V u. s. w. berechnen zu müssen. 

Wenn wir trotzdem nicht bei der Formel (29) stehen 
])leiben, sondern derselben auch noch die übrigen, welche zur 
Berechnung der vierten, fünften, sechsten und siebenten Diffe- 
renzen dienen, folgen lassen, so hat dies seinen Grund vor allem 
darin, dass eine Reihe von theoretischen Untersuchungen über die 
Interpolation bei den mathematischen Tabellen es als wünschens- 
wert erscheinen lässt, jene Formeln in ihrer Gesamtheit kennen 
zu lernen. Die Ableitung der zur Berechnung der Differenzen 

2/1 + 1 2/1 + 1 
(a, + /i ö>, 4), (a -] ^— oj^ 5), (a + /i oj, 6) und (a -\ ^^ — oj, 7) 

notwendigen Formeln gestaltet sich von nun an um so einfacher, 
als die B wegfallen und die Ä ein leicht zu behaltendes Bildungs- 
gesetz erraten lassen. Man findet: 
(fl + /i ö>, 4) = m* {a + n co, 4) + ^"' {a + n oj, 6) -| , (31) 

, , 2/1+1 j,. ^. , 2n+l ,. , .,y. . 2nf 1 ^. , ,^^. 
(«H 2^— ^^,5) ---mHa-i -^(o,ö) + A^'' {a-\ ^^, /)+-., (32) 

(fl + /i oj, 6) -- m^ {a + n ö>, 6) H , (33) 

(a + ^"^^ oj, 7) = m^ (a+ ^'' + ^ co, 7) +.--, (84) 

worin : ^"' ^^4 m^ A^ ^*^' _^ 5 7W* A', 

§ 38. Der Übergang zu den bei der Drei- bezw. Fünfteilung 
des Intervalls zu verwendenden Formeln wird jetzt offenbar 
in der Weise bewerkstelligt, dass wir in den im vorigen Para- 
graphen abgeleiteten Beziehungen für m die speziellen Werte 

,, bezw. - substituieren. Im ersten Falle erhalten wir mit Rück- 

sieht auf (25) aus (24) : 

8 
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§ 39. Die Formeln (35) bis (39) gestatten, vorausgesetzt dass 
die Differenzen von höherer als der siebenten Ordnung vernach- 
lässigt werden dürfen, noch folgende bemerkenswerte Umformung: 
Schreibt man z. B. (39) in der Form 



(«+-f-c,.,5)^243- • 2 
und beachtet (41), so erhält man : 



1 , . 2n4-l ^,5 1 , , 2n+l ^, 



3 2187 



(« 



2/.+ 1 



w 



,6) 



1 



2n+l 



2fl-|-l 



In gleicher Weise findet man statt (38), (37), (36) und (35) nach- 
einander : 



1 4 

(a + /? a>, 4) ^ — (rt + n oj, 4) — g (a + fl a>, 6), 



(43) 



. , 2/. + 1 „^ 1 . 2»+l „^ 



2 



27^- ' 2 

1 - . 2 /» -I- 1 „. 



/ I 2/14-1 -- 



(44) 



12 1 

(rt 4- w oj, 2) — g (a -f /i <w, 4) — -^ (a+fl «>, 6), (45) 



9 



(a -|- /i ö>, 2) = 

(«H ir— «» 1) =^ ö (« H ?^«',1) — q(«+ 



3 



2/1-1-1 



w,8)'.(46) 



2 ' "^ 3 '"'' ' 2 "' ^' 3 "•" ' 2 

Für den Fall, dass in der gegebenen Tabelle (XVII) keine hohem 
als fdnfte und in der gesuchten (XXXVin filr Ä = 3) keine höhern 
als dritte Differenzen in Rechnung gezogen werden müssen,** 
lassen sich die obigen Formeln in folgende einfachere überführen: 



(47) 



, . 2/1 fl „^ 1 r 



^ , 2«+l -^ 1 , 2n + l _. 
{a-\ 2 — ^^)3)- y(aH 2 — "*' ^^ 



(a + nw,2) .U{a + Hoj, 2) 



27 



(a f » cy, 4) 



. , 2n + l 1 

(«I— 3— ,1)^3 



. , 2W4-1 ., , , 2/. + 1 ,, 



§ 40. Übergehend zur Herleitung der Formeln für m = 
erhalten wir: 



1 

5' 



* Der rechnerische Vorteil dieser Formeln gegenüber den ursprünglichen besteht 
allein in der grossem Einfachheit der Koeffizienten. 

** Z. B. bei der Herstellung der astronomischen Jahrbücher (Ephemeriden- 
gammlungen). 
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+ 15625 ('^ +-2 '"'^^ 



52 , , 2n-j-l _., 



(a + r^ ^ , 7) + , m) 



390625 ' ' 2 

1 2 

(« + /i w, 2) = — (a 4- « cu, 2) — g25 (« + »» "*» *) 

33 
'+ 73^25 ia -\- n oj, 6) — , (49) 

, I g/i + 1 „, 1,1 2*1+1 „, 3,1 2n+l ., 
(«+^^«',3) = j25(a + -^-,3)-3^(a + -^«.,5) 

. 57 , , 2w4-l 

1 4 

(« + /JW, 4) = gggCa + w«', 4)— ^25(" + **''''6)+---' (51) 

, ,2/1+1 _, 1 , ,2« + l 

(« + — 2^«>,6) = 3^ (« +— f^a>, 5) 

-15625^«+— 2 "'^^+ ' ^^2) 

(a + « w, 6) = Ygggg (a + « w, 6) — (53) 

und 

,, 8/1 + 1 „ ■ 1 , , 2 « + 1 

{a-\-—f-a,, 7) == ^j2^ (a + — ^^ «, 7) - (o4) 

Zur bessern Übersicht ordnen wir diese Formeln, unter Beifügung 
der siebenstelligen Logarithmen der Zahlenfaktoren, wiederum in 
ein Schema (siehe Tabelle XLII). 

§ 41. Wie bei der Dreiteilung, so lassen sich unter derselben 
Voraussetzung auch die Formeln (48) bis (62) auf eine der Rech- 
nung leichter zugängliche Form bringen. Schreibt man z. ß. (52) 
in der Form 

, , 2/1+1 ., 1 , 2n+l ,., . 1 . , 2n + l ^, 
(« + -2-^,5) -3125(^ + -^^^'^)-^-7M25(^+-T- 

so erhält man unter Berücksichtigung von (64) 

, , 2/1 + 1 ., 1 . , 2n+l ^, ... 2/1+1 „, ,.., 



(Zu Seite 116.) 



{a^n CO, 6) ; (a H — ^^^f- w, 7) 



(«+ 



2/1 + 
2 



(a -f- " < 



, ,2/1 + 



(« + « 






. 33 

'"78125 
[6,6257239] 



15629 
[6,4082400] 



1 



l 



(« + « 



(^ I 2/.-hJ 



"T" 15625 
[5,8061800] 



52 



390625 
[6,1242432] 



57 



~^ 390625 
[6,1641148] 



1 



15625. 
[5,8061800] 



1 



"^78125 
[5,1072100] 
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Ähnliche Umformungen ergeben: 

(o + /i w, 4) = — (a -|- « w, 4) — 4 (« + /i w, 6), (56) 

, .2/1 + 1 _, 1 ^ 1 2w + l g. „, L*H-i R. 
(a -I J— «, 3) =- Y^g (a-| gL— w, 3) — 3 (a H g*— ö>,6) 

18. , 2/1—1-1 „. ,„_^ 

— y(«-| 2 ^'^)' (57) 

(a + /i«,2)==^(a + nß>,2)--2(« + /i«>,4) — ^(a + /ic<i,6), (58) 
, ,8/1 + 1 ,, 1 . 2n + l ,. , ,2/1 + 1 „, 

- -^ (« + — ^ «>, 6). (59) 

Unter der Voraussetzung, dass in der Ausgangstabelle keine höhern 
als die fttnften und in der neu zu erstellenden Tabelle keine höhern 
als die dritten Differenzen berücksichtigt werden müssen, ergeben 
sich die einfachen Fonneln: 

, , 2n+l _, 3 , , 2n+l _, 
(«H 2 ""'^"25^""^ 2 ^' ^ 



(60) 



, ,2/1 + 1 „^ 1 



(a + /i w, 2) = ^Ua -{- n oj, 2) — ^ {a -\- n (o, 4)1 



, . 2/1+1 ., 1 



. . 2«+l ., 1 r , 2n+l _^ 

(.a-\ 2""^' ^~25P"^ 2 '"'^ 

14 ^ , 2n+l _.\) 
-125^« + -^-"'^))) 



3 14 1 

Wenn man bedenkt, dass ^ und ^t^ sich nur um -z-^ = 0,008 

' 25 125 125 ' 

von einander unterscheiden, so kann die dritte Formel von (60) 
mit einer in den weitaus meisten Fällen genügend grossen Ge- 
nauigkeit auch wie folgt geschrieben werden: 

, , 2/i-hl . ■ 1 , , 2w+l .. . . 2/1-1-1 ^^ 

% 42. Die vorstehend entwickelten Formeln, die zum grossen 
Teil an Einfachheit und Bequemlichkeit nichts zu wünschen übrig 
lassen, arbeiten völüg sicher und lösen, wie die nachfolgenden 
Beispiele dartun sollen, unter Berücksichtigung der Bedmgungs- 
gleichungen (21) und (22) das Problem der Drei- bezw. Fünfteilung 
in einer filr die Praxis durchaus befriedigenden Weise. 
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Beispiel 4. Ausgehend von den Werten für 

log sin 18^ 30' -= 9,5014764 

log sin 19^ = 9,5126419 

log sin 19^ SO' =- 9,5234953 

log sin 20^ =- 9,5340517 

log sin 20^ 30' -- 9,5443253 

log sin 21^ -^ 9,5543292 

log sin 21^ 30' =- 9,5640754 

soll eine Tabelle mit einem Intervall von 10' hergestellt iverden. 

Zur Lösung dieser Aufgabe bedürfen wir nach obigem vor 
allem der durch die gegebenen Funktionswerte bestimmten Diffe- 
renzen der verschiedenen Ordnungen. Sind diese einmal berechnet 
(siehe Teil A der Tabelle XLIII), so findet man die ausgezeichneten 
Differenzen der Abteilung B derselben Tabelle durch Benutzung 
der Formehi (47). Das Ausfüllen der Lücken in I', II' und III' 
geschieht in der früher (§ 36) beschriebenen Weise unter steter 
Berücksichtigung der in (21) gegebenen Bedingungsgleichungen. 
Auf vorliegendes Beispiel angewendet, lauten diese letztern wie 
folgt : 

5,3 + 5,8 + 5,2 = (— 314,1) — (— 329,9) = + 15,8 

5,0 4- *,9 + 4,6 = (— 299,6) — (— 314,1) -- + 14,5 

— 319,3 — 314,1 — 309,1 ^ 34244 — 85186 = — 942 

35511 + 85186 + 34867 -- 9,5340517 — 9,5234958 = 105564 
34553 4- 34244 + 33939 = 9,5443258 — 9,5840517 = 102736. 

Nachdem dafür gesorgt ist, dass diese Gleichungen genau erfüllt 
sind, finden sich die gesuchten Funktions werte (Teil C) durch 
blosses Aufsummieren der . Differenzen der Spalte I'. Darnach 
wird z. B.: 

log sin 19« 40' = 9,5284958 + 0,0035511 = 9,5270464. 

Die Vergleichung der gefundenen Funktionswerte mit denjenigen 
einer Logarithmentafel zeigt vollständige Übereinstimmung. 

Beispiel 5. Der Nautical Almanax) für das Jahr 1906 gibt für 
dm Planeten Mars nachstehende Werte in Rektaszetision : 

Aug. 24. 9* 22'" 51,3V 
Aug. 29. 9^ 35^^' 27,17' 
Sept. 3. 9* ^?" 54,35' 



1 



119 

Sept, 8. 10^ 0^ 1S,60' 
SepL 13, 10^ 12^ 25,20' 
SepL 18, 10^ 24"^ 29,94' 
SepL 23. 10^ 36"' 28,10". 
Man berechne hieraus eine Tabelle^ deren Intervall einen Tag beträgt. 

Die gewünschte Tabelle (Teil C der Tabelle XLIV) wird in 
genau derselben Weise erhalten, wie diejenige des vorhergehenden 
Beispiels. Die Berechnung der ausgezeichneten Differenzen erfolgt 
mit Hilfe der Formeln (60). Die Bedingungsgleichungen (22) gehen 
über in: 

0,0050« + 0,0048« + 0,0046« + 0,0045« -f 0,0044^ 

= (— 0,2977«) — (— 0,3210«) = + 0,0233« 

0,0043« + 0,0042« + 0,0039« + 0,0038« + 0,0035« 

= (— 0,2780«) — (— 0,2977«) = + 0,0197« 

— 0,3066« — 0,3021« — 0,2977« — 0,2934« — 0,2892« 

= 2"^ 26,336« — 2°» 27,826« = — 1,489« 

2°» 28,452« + 2r 28,136« + 2°^ 27,826« + 2r 27,518« + 2"* 27,216« 

= 10»^ 0"» 13,60« — 9»^ 47°^ 64,36« -- 12^ 19,15« 
2"* 26,918« + 2«^ 26,625« + 2°^ 26,336« -f- 2"* 26,051« -f 2"^ 25,770*» 
= IQ^ 12"^ 26,20« — 10»^ 0°^ 13,60« = 12"^ 11,70«. 

Die gesuchten Funktionswerte ergeben sich wiederum durch 
Aufsummieren der Differenzen I'.* 

§ 43. Laut der der Tabelle XXXVIII zu Grunde gelegten 
Voraussetzung kann für k jede beliebige ungerade, positive Zahl 
gesetzt werden, mithin ausser 3 und 5 auch 7, 9, 11 u. s. w. 
Nun lehrt aber die Erfahrung, dass der mühsamste Teil der bei 
der Verwendung des vorstehend besprochenen Teilungsverfahrens 
erforderlichen Arbeit gerade im Ausfüllen der zwischen den direkt 
berechneten Differenzen sich befindlichen Lücken besteht. Dieser 
Tatsache eingedenk, wird der Berechner von Tabellen gut tun, 
das rohe Intervall, wo immer es angeht, so zu wählen, dass er 
sich bei der Verkleinerung desselben auf die Drei- resp. Fünf- 
teilung beschränken kann. 

Wünscht man z. B. eine von 10 zu 10' fortschreitende Sinus- 
tafel (vergl. Tabelle III der Schlömilch'schen Logarithmentafel), so 
ist es vorteilhaft, die Rechnung etwa folgendermassen anzulegen : 



Vergl. Nautical Almanac 1906, Seite 239. 
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Man verschafft sich die Werte für sin 4», sin 8^, sin 12» u. s. w. 
Hierauf verkürzt man das Intervall durch wiederholtes Inter- 
polieren in die Mitte nach und nach auf -, -, - des ursprüng- 

J 4 o 

hchen Intervalls, also auf 2% 1», 30'. Einmal so weit, erhält 
man dann die verlangte Tabelle mittelst einer Dreiteilung, welche 
hier umso leichter durchgeführt werden kann, als durch die vor- 
ausgehenden Verkleinerungen die Ordnung der noch in Betracht 
kommenden Differenzen bedeutend herabgemindert wurde. 

§ 44. Der Umstand, dass die Stufenzahlen der beiden heute 
bei uns im Gebrauche stehenden Zahlensysteme Vielfache von 
3 und 5 smd, enthebt den Rechner von vornherein der Not- 
wendigkeit, dem k höhere Werte als die eben genannten beilegen 
zu müssen.* Wenn wir nachstehend trotzdem ftlr die Zehn- 
teilu7ig ein selbständiges, somit von der Fünfteilung unabhängiges 
Verfahren vorschlagen, so geschieht dies namentlich auf Grund 
der Tatsache, dass sie sowohl bei Tabellen mit sexagesimaler, 
als auch bei solchen mit dezimaler Grundlage gelegentlich mit 
grossem Vorteil vorgenommen werden kann. Die hierbei erfor- 
derlichen Rechnungen sind im wesentlichen dieselben, wie wir 
sie oben für ungerades k auszuführen genötigt waren (vergleiche 
§ 36 u. ff.); auch hier berechnen wir mittelst eigens zu diesem 
Zwecke herzuleitender Formeln eine Reihe von Differenzen, auf 
welche die ganze Rechnung abschliessen muss. 

Mit Rücksicht darauf, dass die Zehnteilung nur dann mit 
Erfolg angewendet wird, wenn höhere als dritte Differenzen auf 
das Resultat keinen Einfluss mehr auszuüben vermögen, be- 
schränken wir uns bei der Ableitung der Formeln auf dritte 
Differenzen in der geauchtep und auf vierte Differenzen in der 
gegebenen. Tabelle.** Im übrigen sollen auch jetzt die direkt zu 
berechnenden Differenzen durch fettere Schrift kenntlich gemacht 
sein. 



* So lassen sich z. B. die häufig notw endiget Sechs- und Zehnteilungen durch 
«Ine vorausgehende Interpolation in die Mitte auf. eine Drei-, bezw. !^iinfteilung 
zurückführen. 

** Die Ausdehnung der Formeln auf höhere Differenzenordnungen bietet den 
Studierenden nach Kenntnisnahme der in § 45 gegebenen Anleitung nicht die 
geringsten Schwierigkeiten. 
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§ 45. Da im vorliegenden Falle die Anzahl der Teile, in 
welche das ursprüngliche Intervall zerlegt werden muss, eine 
gerade ist, so ist ohne weiteres klar, dass die zu berechnenden 
Differenzen ungerader Ordnung der gesuchten Tabelle nicht, wie 
es für ungerade k der Fall war, auf den zwischen je zwei auf- 
einander folgenden Funktionswerten der gegebenen Tabelle gezogen 
gedachten Horizontalen hegen können. Diesem Umstand ist es 
denn auch zuzuschreiben, dass die in § 37 zur Berechnung der- 
selben gefundenen Beziehungen hier nicht zur Verwendung ge- 
langen dürfen. 

Indem wir mit k wiederum allgemein die Anzahl der Teile 
bezeichnen, in welche das ursprüngüche Intervall geteilt werden 

soll, erhalten wir, die Funktionswerte fia-{-^oj)un6f(a-\ — ^~- cd) 
als Ausgangswerte wählend, 

oder, wenn wir wie in § 37 Ä = — setzen, 

. . l-\-m ,l4-2m. ^, .i . 
(fl H ^ — a>, 1) = /-(a -] L_ aj)—f{a-}-^ w). 

Unter Benutzung der Formel (12) ergibt sich hieraus: 



(«, 0) + {^±^\ (« + >,!) + {"^X ia, 2) 

+ (^^+ 1)3 (« + i «>, 3) + (i^+ 1\ ia, 4) 
- («' 0) - (|)x (« + T «>, 1) - {^\ ia, 2) - (^+ 1)^ (a + I «,, 3) 



■2 

.2 /vyiS fiiTii 



^m(a + ^£^, 1+-^- <a, 2) + (g- -^j — 24)^^ + IT^' 3) 

, m* 5'wa\ 

+ V2i-^(^'*^ 

oder kürzer: 
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(« + 



mit 



1 -f zu 



a>, l)=m (a + -i-ö>, l) + ^(a, 2) + jB(a + |co, 3) 

+ C(a, 4) (61) 



(62) 



Ä=^ 



B 



6" ' T ~24 



(7 = 



m. 



öm' 



24 48 

Die Berechnung der Differenzen gerader Ordnung erfolgt nun 
wieder nach den in § 37 gefundenen Formeln. Darnach findet man : 

(fl, 2) == m^ (a, 2) + A' (a, 4) (63) 



mit A' 



12 



(m^ — 1). 



Um endUch noch eine Formel zur Berechnung der dritten 
Differenzen zu erhalten, erinnern wir daran, dass dieselben aus 
den Differenzen zweiter Ordnung in genau gleicher Weise abge- 
leitet werden, wie die ersten Differenzen aus den Funktions- 
werten. Bilden wir demgemäss die 

Tabelle XLV. 



Funktiofi 


i. Differenz 


2, Biffereyiz 


• 

m2 (a — a>, 2) + A' {a — ö>, 4) ' 

m2 (a, 2) + A' (a, 4) 

1 

m^ (a 4" ö>, 2) + A' (a -\- w^ 4) 

• 

1 


1 

1 

mM« — 1 w, 3) 
m-^ {a + Y ö>, 3) 

• 


• 
^2 (a — ö>, 4) 

^2 ((7,^ 4) 

771 2 (a -|- ö>, 4) 

• 



und berücksichtigen die Beziehungen (61) und (62), so wird 

= m^ {ß -{- \ CO, ?,) -\- — (a, 4). 



(« + ^-.3) 



(64) 
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Substituiert man jetzt in den drei gefundenen Formeln (61), (63> 
und (64) für m den Wert t^* und vertauscht . a gegen a -\- n o> 
(n = ..., — 3, — 2, — 1, 0, -|- 1, -|" 2, + 3, ---), so finden wir: 

3 , . 2»t+l -^ 83 ■ . 



(65) 



2000' ' 2 ' ' 80000 

1 33 

(a + /iß>,2) = __ (ör -|- „ w, 2) — ^^^^ (a + n «, 4) 

1 , , 2w+l _, , 1 



Mit Hilfe dieser Formeln gelingt es nun ohne Mühe, die 
ausgezeichneten Differenzen zu berechnen und so das Gerippe 
des Differenzenschemas der gewünschten Tabelle zu bilden. Das 
Ausfüllen der Lücken geschieht in der früher beschriebenen Weise 
unter steter Beachtung nachstehender Beziehungen: 

I 2 % — — 1 

(a-\-n (o-\ __ö>,3) = (a4-/i+lc/J,3)— (a-|-/iö>,2> 



(66) 



2J<- I -" . 20 
» = i 

«~5 



V'(a4-iiö> + -^ö>,2)=(a4-/i+0,65o>,l)— (a+/i— 0,46o>,l> 



y, , , 2.-1 



(a + ^ oj -^ 20^^"^'^^ ^ (fl + /i + 1 ^t>, 0) — [a-{-nw, 0> 



»=.--,- 3, — 2, - 1, 0, + 1, + 2, + 3, ■-. 

Baispiel 6. Der „Nautical Almanac" für dns Jahr 1906 gibt 
für den Planeten Saturn nachfolgende Werte in Rektaszension : 

Mai 7. 23'' 0>» 4,92' 

12. 23 1 27,92 

17. 23 2 43,64 

22. 23 3 51,72 

27. 23 4 51,75 

Juni 1. 23 5 48,42 

6. 23 6 26,52 

11. 23 7 0,84. 

Man berechne hieraus eine Tabelle, deren Intervall 0,5 Tag beträgt^ 



* Es versteht sich von selbst, dass m nicht aoE den Wert -^ beschränkt ist ; im 
G^enteil kann dafür der reziproke Wert jeder beliebigen geraden Zahl gesetzt werden.. 
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Die Herstellung der gewünschten Tabelle (XL VI) erfolgt 
in derselben Weise wie in den Beispielen 4 und 5 mit Hilfe 
der Formeln (65) und unter Berücksichtigung der Bedingungs- 
gleichungen (66). 

Alle bis jetzt beschriebenen Verfahren, betreffend die Teilung 
•der Intervalle, wickeln sich, wie eine mehrfache Erprobung bei 
der Verfeinerung von Tafeln gezeigt hat, äusserst rasch und glatt 
ab. In nicht ganz demselben günstigen Lichte dagegen zeigen 
sie sich, wenn man die Funktionswerte der durch wiederholte 
Teilung erhaltenen Tabellen auf ihre Genauigkeit prüft. Gele- 
gentlich vorgenommene Proben beweisen, dass die Sicherheit der 
letzten Stelle sehr häufig zu wünschen übrig lässt. Diesem Übel- 
stand kann indessen leicht dadurch abgeholfen werden, dass man 
die direkt berechneten Funktionswerte mit einer grösseren Anzahl 
von Dezimalen tabuliert, als schliesslich beibehalten werden. 

§ 46. Die vorstehend angestellten Betrachtungen über die 
Herstellung mathematischer Tabellen gelten allgemein, mithin 
auch für den Fall, dass die zu tabulierenden Funktionen ganze 
rationale sind. Dessen ungeachtet geben wir für solche ein 
anderes Verfahren und zwar deshalb, weil das Konstantwerden 
der Differenzen einer bestimmten Ordnung wesentliche Verein- 
fachungen gestattet. 

Um die Vorstellung zu fixieren, denken wir uns eine ganze 
Funktion F {x) vom dritten Grade, so dass nach § 5 die Grössen 
d der Tabelle I durchweg denselben Wert di besitzen. Dann lässt 
sich, wenn ausserdem noch die Anfangsglieder ai, &i und Ci bekannt 
sind, die Tabelle durch blosses Aufsummieren berechnen. Aus der 
Entstehungsart des Differenzenschemas ergeben sich nämlich die 
Beziehungen : 

<?2 = Ci -f- rflj 03=^^2+ ^2, Ö4 ===^3 + ds, Co = C4 + ^*7 ••• 
62 = fti -j- Gli &3 ^^ &2 4" ^2, Ö4 =' b'i -{- C3, 65 == 64 -[- C4, •-• 
a2 = Ol + &i, a-i = (h + ^2, «4 = 0^3 + Ö3, «6 = »4 + 64, ••• 

Als numerisches Beispiel diene die Berechnung der den Argu- 
menten X* = 0, 1, 2, 3, ••• entsprechenden Werte der Funktion 

F {X) = l()x^— 101 x^ —imx-[- 1799. 

Js[ach Tabelle X ist: 

. a,=-\- 1799, 61 = — 200, Ci -- — 142 und dk = +.60. 
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Mai 17. 

17,5. 
18. 
18,5. 
19. 
19,5. 
20. 
20,5. 
21. 
21,5. 
33. 
22,5. 
28. 
28,5. 
24. 
24,5. 
25. 
25,5. 
26. 
26,5. 
37. 
27,5. 
28. 
28,5. 
29. 
29,5. 
30. 
30,5. 
31. 
31,5. 
Juni 1. 



'. I 
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Das Aufsummieren geht nun in der Weise vor sich, dass man 
in die Spalte der di, ck, c?3, ^4, •-• die Zahl 60 und an die Stelle 
von Ci den Wert — 142 einträgt. Die übrigen c berechnet man 
dann mit Hilfe der oben gegebenen Beziehungen. In gleicher 
Weise ergeben sich hierauf auch die b und die a. 

Das hier dargelegte Verfahren ist an sich durchaus brauchbar^ 
es hat jedoch den Nachteil, dass die zu addierenden Grössen nicht 
senkrecht über einander stehen. Um demselben wirksam entgegen- 
zutreten, verfahren wir folgendermassen : Wir bilden ein von 
hnks nach rechts sich ausdehnendes Rechenschema mit vier- 
zeiligen Spalten. In die erste Zeile kommen die Zahlen d, die 
nach der Voraussetzung sämtlich den Wert di == 60 besitzen. 
In die zweite, dritte und vierte Zeile der ersten Spalte (von links) 
setzt man der Reihe nach die zum voraus berechneten Grössen 
Ci = — 142, bi = — 200 und ai = 1799. Ist dies geschehen,, 
so stellt sich die Rechnung wie folgt: 



60 


60 


60 


60 

-|-38 

— 446 


60 


60 


60 


60 


• 


142 
200 


82 


22 


-|-98 


4-158 


-|-218 


-1-278 


• 


342 


— 424 


— 408 
-1-387 


310 


152 


+ 66 


• 


-1-1799 


-{-1599 


4-1257 


-1-833 


21 


331 


483 


• 



Hierin ist z. B. : 

— 424 = — 82 + (— 342), — 446 -- — 22 + (— 424) 
-f- 387 == — 446 + 833, — 21 -= — 408 + 387. 

Es bedarf wohl keiner besonderen Auseinandersetzungen, wie die 
Rechnung anzulegen ist, wenn statt der dritten Diflferenzenreihe 
eine andere als Ausgangspunkt zu dienen hat. 

§ 47. Unsere nächste Aufgabe wird nun darin zu bestehen 
haben, zu zeigen, wie man die Gheder 6i, Ci, di u. s. w. der von 
Ol ausgehenden absteigenden Treppe des Differenzenschemas I 
berechnen kann. 

Bei der Lösung dieser Aufgabe gehen wir von der Annahme 
aus. die zu tabulierende Funktion sei vom n^" Grade. Damit 
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liefern uns die in § 6 hergeleiteten Beziehungen für J F {x\ 
J^F(x)j J^Fix), '", J''F{x) unmittelbar die gewünschten Grössen, 
bedeuten doch die hnken Seiten jener Gleichungen nichts anderes, 
^Is die mit &i. Ci, di, ••- bezeichneten Glieder des nach rechts bis 
zur n^"^ Differenzenreihe erweiterten Schemas der Tabelle I. 

Um mit der bisher üblichen Bezeichnungsweise Übereinstim- 
mung zu schaffen, ersetzen wir 

Fia), J F{a), J^ F (a), J^ F (a), .-., J« F (a) 

durch die in § 18 eingeführten Symbole 

<a, 0), (a + ^ io, 1), (a+ I w, 2), (a+ |\.>, 3), .-., (a +|- w, n). 

So gehen die erwähnten Beziehungen, wenn man in denselben 

"^ F 
noch allgemein fiir F^ (x) das Symbol ^— ^ substituiert, über i 

die folgenden: 



m 



<67) 



|«.+i-.,i,=4'>(|f).» + ^S'(P.<-+^"(||).«"+- 

(»+|,..2)=4»(||).-».+4»(|^).,..+...+^"(0)...-, 

^"(II).'-+-+^?'(ID.'-. 



(a + -|^e>,3) 



dx^ 



d^* 



(a-\^^w,n) 



"^'^ Vax«/''' ' 



in welchen die Koeffizienten Ä durch die in § 6 gegebene Be- 
ziehung (2) vollständig bestimmt sind. 

§ 48. Wenn wir trotzdem bei der numerischen Berechnung 
der in Frage stehenden Koeffizienten nicht jene Beziehung ver- 
wenden, so geschieht dies lediglich deshalb, weil es möglich ist, 
einfache Formeln abzuleiten, welche jedes beliebige Ä direkt zu 
liefern im stände sind. Um diese Formeln aufzufinden, verfahren 
wir folgendermassen : 

Nach obigem erhält man z. B. für die Differenz der q^'^ Ord- 
nung (q < 71) : 
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Vertauscht man hierin a gegen a -\~ w, so wird 



+ • 

Andrerseits findet man, da 

(a-{- io-^ ^(0^ q) = (a-{-±w, g) + (a + ^ eo, g + 1) ist, 

+-!*"(P).--+^'*"(P)/»-+' 
=^"(||).'»'+(4'i-+A^-")(|S^)/-- 

+(41.+4«')(|S?)/"'"'+- 

Durch Vergleichung der Koeffizienten der beiden für {(i-\-<o-{-^(ü^q) 
gefundenen Ausdrücke ergeben sich nachstehende Beziehungen: 
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_ y<(9) 



iW) 



Aq) I Aq+1) 
^-j-2 "T ^q-i-2 



Aq) Aq) 



Ol 1^ 



2! 



1! 



oder 






.4 



(q) 
q ? 

(?) 

9 I >«(9> 



= ^ + A+1» 



Aus der ersten dieser beiden Beziehungen folgt unmittelbar, dass 
Mit Rücksicht hierauf erhält man endlich : 

^«•j-.^ — ^ -f- -TTT" ~r A4-27 



(68) 



lq+3 



J^' 



?+4 



3! 
4! 



2! 

3 j I Ol I ^«+3? 



1(9) 



2! 

1(9) 



(9) 



A+5 — 6!-t--4r+'-3r+"ir'+"^+'' 



Ersetzen wir jetzt hier q nach und nach durch alle zulässigen 
Werte und beachten, dass nach (1) in § 6 



4^> 



aW— ^ aW 



3r 



41 



,A'' 



5!' 



i(i) 






(69) 



2!'""° 3!''*** 4!''*" 5!' ' "" nl 

ist, so können alle übrigen A ohne Schwierigkeit berechnet 
werden. So erhält man z. B. für g = 1 : 



^f=^+4'' 



-+-=1 

2! '2! ' 



A?- 



3! 



4' 



a!^' 



2! 

id) 



+ 4" 



1 . .J_ I A_l 
3!^2!2!^3! 12' 



KD 



i-a 



mT- qi -r 2! "r-^4 4!~'~3! 2! '"2! 3!"f'4! 4' 



4! 



3! 



u. s. w. 



Vergleiche die Anmerkung * zu (3) in § 6. 
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Mit Hilfe der in (68) zusammengestellten Recursionsformeln und 
der in (69) gegebenen Werte für Ä^\ J^\ ••-, A^P gelingt es nun 
auch ohne Mühe, die eingangs dieses Paragraphen erwähnten 
Beziehungen aufzufinden. Aus der ersten der Gleichungen (68) 
ergeben sich die folgenden: 



^fit;>- 


- 2t~r A+i) 


4^-- 


1 j(v 1) 

"2! "^^ ' 


j<r'^ 


. 1 1 Ai-p 


A^^ - 


-2, 1 A,, 


4^P 


1 



2!' 
somit durch Addition: 

oder, wenn ^ -|- 1 durch q ersetzt wird, 

4'it - |. (70) 

Aus der zweiten Gleichung von (68) erhält man ferner: 

(9) 



^7+3 — ^ -r ivr I A+2, 



3! ' 2! 



.4,/+2 — ^ -| 5^ h A/-(-i , 



1 j(9-2) 



A+1 — Tjri — ^5T~+^ 



<.'/-2) 



-^4 — ^f+ ITT +^3 , 



3! ' 2! 



3!' 

-- ^■'=^ + il'^*, = ^ + |iy'i Ode.- 
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^^3 , 3! ' 2! 4^*^ ^ 24 • 

Ersetzt man wiederum q-\- l durch q, so wird 

In gleicher Weise ergibt sich aus 

4^3 = ^^^; (72) 

aus 



5! • 4!Z-J "^^ ' 3!^ '"^' • 2! 

.(.) ^ g(15gB-f-30g^ + 5g-2) , 

u. s. w. 

Die hier gefundenen Ausdrücke für ^g?fi, ^?^2, -4j$^3? ^^^, u. s. w. 
lassen wohl zur Genüge erkennen, dass diese letztem hein ein- 
fdches Bildungsgesetz befolgen, dass also die Kenntnis der all- 
gemeinen Form irgend eines Koejßizienten nicht auf diejenige des 
nächstfolgenden schliessen lässt. 

In nebenstehender Tabelle XLVII geben wir die Zahlenwerte 
derjenigen A, welche bei der Tabulierung einer Funktion zehnten 
Grades in Rechnung gezogen werden müssen. Die Anordnung 
geschah in der Weise, dass die Koeffizienten ohne Änderung 
der Reihenfolge gleich in die zur Berechnung der Differenzen 
(a + Y ^'^ l)j (öt + f <f^, 2), (ö^ + f ^'>, 3), u. s. w. in (67) gegebenen 
Reihen übertragen werden können. 



« 


AI) 
-^1 






1 


1 


1 




1 


1 






1 


1 


2 




1 

l 

2 


1 


1 








3 




3 

t 


1 


2 




4 

i 


1 


2 








5 




5 1 


1 


2 




6 


1 1 

1 
t 


3 






i 


7 




7 


1 


2 




1 

8 


1 

1 

1 

1 


4 


: 




1 


9 


' 


9 


1 






1 

1 ( 




2 




r t 

10 

; 1 


1 
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Als numerisches Beispiel diene die Tabulierung der Funktion 

iP (X) = 3 a;* — 5 a;2 -f- 7 

für die Argumente 0, 1, 2, 3, 4, -••. 

Um die in Betracht kommenden Differenzen zu erhalten, 
setzen wir das in den Formeln (67) auftretende a gleich und finden : 

J' (0) = 7 



(1^) = (12 x^ — 10 x)o = 



d'F 



(S6x^— 10)o = — 10 



(Lr)o - (72 .)„ = 



^2 ;c* ^ 







72. 



Es wird somit: 

(a + I w, 1) -= — 5 -f- 3 -= — 2 
(a + I <M, 2) = — 10 + 42 = -f- 32 
(a + 1 (0, 3) = + 108 

(a + j w, 4) = + 72. 

Mit diesen Werten bilden wir jetzt das in § 46 besprochene 
Eechenschema : 



+ 72 


.72 


72 


72 


72 


72 


72 


72 


• 


+ 108 


180 


252 


324 


396 


468 


540 


612 


• 


+ 32 


140 


320 


572 


896 


1292 


1760 


2300 


• 


2 


30 


170 


490 


1062 


1958 


3250 


5010 


• 


+ 7 


5 


35 


205 


695 


1757 


3715 


6965 


• 



Hierin enthält die letzte Zeile die gesuchten Funktionswerte. 

§ 49. Ausser dem vorstehend auseinander gesetzten und bei 
der .Tabulierung von F (x) = S x* — b x^ -{-7 angewandten Ver- 



J^2 

fahren stehen uns noch zwei andere zur Verfügung. Das eine 
besteht darin, dass man so viele aufeinander folgende Funktions- 
werte direkt berechnet, als notwendig sind, um die übrigen durch 
Aufsummieren der aus jenen abgeleiteten Differenzen erhalten 
zu können. Da es nach § 5 möghch ist, die konstante Differenz 
unmittelbar hinzuschreiben, so ist die Zahl der zu berechnenden 
Funktionswerte gleich dem Grad der zu tabulierenden Funktion. 

Das andere Verfahren ergibt sich auf Grund der in § 5 an- 
gestellten Betrachtungen über die Art dei' Differenzenbildung 
einer vorgelegten Funktion und wird folgendermassen gehandhabt: 
Ausgehend von der gegebenen Funktion F (;v) bildet man, indem 
man in derselben statt x successive die Werte x -j- <w, x -{- 2 ov, 
X -f- 3 CO, X -{- 4: o)^ u. s. w. substituiert, die aufeinander folgenden 
FunktJons werte i^ (x -f- fo), F ix -|- 2 «>), F (x -{- 3 co), F (x-^-^t co), 
u. s. w. und bildet hieraus nach der in § 4 und ff. gegebenen 
Anleitung die Differenzen der verschiedenen Ordnungen. Alsdann 
ersetzt man in F {x) und den Gliedern der in F {x) beginnenden 
absteigenden Treppe x durch irgend einen Ausgangswert und 
berechnet die ihnen entsprechenden Zahlenw^erte. Der weitere 
Verlauf der Rechnung ist der gleiche, wie der in § 46 beschriebene. 

Zur Erläuterung des Verfahrens diene folgendes Beispiel: 
Man wünscht mie Tabelle zu konstruieren, ivelche die Kreisflächen 
für die beständig um ^ cm wachsenden Halbmesser zu liefern im 
Stande ist,* 

Bezeichnen wir den Halbmesser irgend eines der Kreise mit x^ 
so sind wegen lo =^ j die Halbmesser der nachfolgenden Kreise 

bezw. .'^ + T' '^ ~f~ Y' ^ H" T' ■^' + 1? ^- s- ^M ^^^ demnach die 
Flächeninhalte x'^ :r, (x + |)2 ;:, {x-{-~y' r., {x + |)2 ;:, {x -\- 1)^ ;:, 
u. s. w. Mit diesen Grössen bilden wir jetzt nachstehendes Schema: 



* Vergleiche : Die arithmetischen Keihen höherer Ordnungen und die figurierten 
Zahlen, von Dr. Anton (PrOgrammschrift des Oelsmischen Gymnasiums, 1850). 
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Arg. 



Funktion 




X 



/V» 2 T 






4 

2 



6 



(X + ^-)ä r -- (.c^ + g- •* + 16^ '^ 
2 1 

(U; 4- -)2 TT = u-ä + .t + t) " 

4 4 

(X -I- ^)ä TT = {.c-^ _|- -_ ,c + — ) ;: 

ix + l)ä - =- U'^ + 2 .« + 1) r 

5 ö 25 

(o; + j)^ - -= (.c^ + k •*' + 16^ '^ 



6 



9 



^ + J I' (a; + j)ä5r = (:«2 + 3.c + -j-);: 



1 , 1 
2^+16 

1 , 3 

2^ + 16 
1 , 5 

2*"^"16 

1 , 7 

2*^~l6 

1 , 9 

2^"f"l6 

1 ,11 

2 •''+16 



•iZ 



TT 



7t 



n 



TT 



welches für x -=0 in folgende Tabelle übergeht : 

Tabelle XLVIII. 



Hälbni, 



Fläche 



n 



■ ii 

0,25 i| 0,19634954 

0,5 ;! 0,78539816 

0,75 I 1,76714586 

1 :j 3,14159264 
1,25 ; 4,90873850 

I 

1,50 I 7,06858344 

1,75 i 9,62112746 

2 12,56637056 
2,25 15,90431274 



0,19634954 

0,58904862 

0,98174770 

1,37444678 

1,76714586 

2,15984494 \ 

2,55254402 

2,94524310 

3,33794218 



0,39369908 

0,39269908 
0,39269908 
0,39269908 
0,39269908 
0,39269908 
0,39269908 
0,39269908 
0,39269908 
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Der Umstand, dass die Ausgangswerte der vorstehenden 
Tabelle infolge der Irrationalität von ;r mit Abrundungsfehlern 
behaftet sind, lassen es begreiflich erscheinen, dass die durch 
das Aufsummieren jener fehlerhaften Differenzen gefundenen 
Funktionswerte um so mehr von den direkt berechneten ab- 
weichen müssen, je weiter man einerseits mit der Ordnungsniimmer 
der Smnmenreihen und andererseits mit der Ausdehnung der Tabelle 
vorwärts geht. Um dies' einzusehen, wollen wir der Einfachheit 
halber vorerst annehmen, es seien nur die Differenzen der höchsten 
Ordnung mit einem Abrundungsfehler so behaftet. Unter dieser 
Voraussetzung und mit Rücksicht auf den Satz, dass die Fehler 
der Differenzen zugleich die Differenzen der Fehler sind (vergL 
§ 12), betrachten wir die 5o als die konstanten Differenzen, aus 
welchen durch successives Aufsummieren die Differenzen nullter 
Ordnung — die Funktionswerte — hergeleitet werden und erhalten 
das Schema 



«0 


So 


0»> 

»0 


«0 


So 


«0 


«0 


«0 


So 


• 





«0 


2 s« 


3 CO 


4£o 


5 CO 


6 So 


7 s« 


8 s« 


• 








«0 


3 So 


6 So 


10 So 


15 So 


21 So 


28 So 


• 











So 


4£o 


10 CO 


20 So 


35 So 


56 s« 


• 














So 


5^0 

^0 


15 So 


35 So 


70 s« 


• 

















6 s« 


21s« 


56 So 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 



durch welches die Richtigkeit des oben ausgesprochenen Satzes 
offenbar dargetan ist. Ein Bück auf vorstehendes Täfelchen lasst 
ausserdem erkennen, dass die darin vorkommenden Grössen aus 

n {n—V) (w — 2) - - - (n — Ä + 1) 



{n)k • So = 



^0 



1.2-3 k 

dadurch abgeleitet werden können, dass man sowohl dem k als 
dem n alle möglichen ganzen, positiven Zahlenwerte beilegt. So 
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findet man z. B. die Glieder der dritten Zeile, d. i. der zweiten 
Summenreihe, für /c = 2 und w = 0, 1, 2, 3, 4, u. s. w. Entgegen 
der oben gemachten Annahme wollen wir jetzt voraussetzen, es 
seien auch die übrigen, durch Rechnung gefundenen Differenzen 
mit Abrundungsfehlem in das Rechenschema eingetragen worden : 
Bezeichnen wir diese der Reihe nach mit si, ^2, «3, u. s. w., so 
zeigen ähnliche Untersuchungen, wie die soeben angestellten, dass 
,-das (n -f- 1)^ Ghed der (k -\- 1)'^" Zeile des Fehlerschemas durch 

{n)k • So + (w)jt_i • ei -f (n)A-_2 • h + (^)^-3 • «3 + (^)i-4 • -* H 

bestimmt ist. 

Aus dem gesagten ergibt sich für den Rechner folgende wichtige 
Verhaltungsmassregel : Wenn die Benutzung völlig scharfer Werte 
aus irgend welchen Gründen als nicht enipfehlenstvert erscheint, 
so ist bei der Anwendung des vorstehend beschriebenen Summen- 
verfahrens die Anzahl der mitzunehmenden Dezimaien derart zu 
bemessen, dass der Einfluss der Abrundungsfehler unschädlich bleibt. 

§ 50. Wir machen nun im fernem noch die Annahme, es 
hege eine fertige Tabelle irgend einer ganzen, rationalen Funktion 
vor uns und es handle sich darum, dieselbe in eine andere, mit 
einem /c-mal kleineren Intervall zu verwandeln. Es versteht sich 
von selbst, dass die in den §§ 36 bis 45 besprochenen Teilungs- 
verfahren auch hier zur Verwendung kommen können und zwar 
um so mehr, als die Konstanz der GUeder der letzten Differenzen- 
reihe wesenthche Vereinfachungen sowohl in Bezug auf die 
vorbereitenden Rechnungen (Berechnung der ausgezeichneten Dif- 
ferenzen) als auch hinsichtiich des Ausfüllens der Differenzen- 
spalten gestattet. Zur Erläuterung möge die 

Tabelle XLIX. 



X 


F(x)=x^ 


I 


n 


ni 


1 

2 
3 
4 
5 
6 

• 


1 

8 

27 

64 

125 

216 

• 


7 
19 
37 
61 
91 

• 


6 
12 
18 
24 
30 
86 

• 


6 
6 
6 
6 
6 

• 
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enthaltend die dritten Potenzen der Zahlen 1, 2, 3, 4, u. s. w. in 
eine Tabelle mit 10-mal kleinerem Intervall verwandelt werden, 
d. h. in eine solche, welche die dritten Potenzen von 1,1, 1,2, 1,3, 
1,4, u. s. w. zu liefern im stände ist. 

Indem wir Fil)--^ 1 als Ausgangswert betrachten und ihn 
demnach, gemäss unserer bisherigen Uebung, mit (o, 0) bezeichnen 
können, berechnen wir mit Hilfe der in § 45 gegebenen Formeln 
(66) die Differenzen 

(o + 0,56 w, 1), (o, 2) und (o -|- 0,66 w, 3). 
Man findet: 



(a + 0,66 <«, 1) 

(«,2) 

(o + 0,56 <y, 8) 



- 1 



10 
1 

100 ■ 

1 
1000 



7+T 



6 



200 
0,06 



6 



3 



2000 



• 6 -= 0,721 



• H 



0,006. 



Die übrigen Zahlen der neuen Tabelle ergeben sich jetzt in der früher 
beschriebenen Weise successive mittels Addition und Subtraktion: 

Tabelle L. 



X 



Fix) 



-■- X» 



V 



n' 



1 

1,1 

1,2 
1,3 

1,4 
1,5 

1,6 

1,7 

1,8 

1,9 

2 

2,1 

2,2 

2,3 

2,4 



1 

1,331 
1,728 
2,197 
2,744 
3,375 
4,096 
4,913 
5,832 
6,859 
8 

9,261 
10,648 
12,167 
13,824 



0,331 
0,397 
0,469 
0,547 
0,631 
0,731 
0,817 
0,919 
1,027 
1,141 
1,261 
1,387 
1,519 
1,657 



0,06 

0,066 
0,072 
0,078 
0,084 
0,090 
0,096 
0,102 
0,108 
0,114 
0,120 
0,126 
0,132 
0,138 
0,144 



in' 



0,006 
0,006 
0,006 
0,006 
0,006 
0,006 
0,006 
0,006 
0,006 
0,006 
0,006 
0,006 
0,006 
0,006 
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Während wir bei der Herstellung der TabeUe XL VI genötigt 
waren, dem n der Formeln (65) so viele Werte beizulegen, als 
gegebene Funktionswerte vorlagen, genügte es in vorstehendem 
Beispiel, für jede Diflferenzenspalte je ein Glied zu berechnen — 
ein Vorteil, der bei der Herstellung von ausgedehnten Tafeln 
nicht hoch genug angeschlagen werden kann. 

§ 51. Ein zweites, ausschliesslich bei der Tabulierung von 
ganzen rationalen Funktionen anwendbares Verfahren besteht 
darin, dass man die Glieder der mit (a, 0) beginnenden abstei- 
genden Treppe der neuen Tabelle als Funktionen derselben Glieder 
der gegebenen Tabelle darstellt. Zur Herleitung der allgemeinen 
Form diese!' Funktionen bedienen wir uns der 



Tabelle LI. 



Arg. 



Funktion 



a 

a 

a 
a 



k 
2 



w 



k 
3 



(O 



k 



CO 



,4 



a 



k 



(0 



L 6 
' Ar 



(«,0) 

, 2 
(« + ^ «, 0) 

(« + 1 ^"^ 0) 

4 
(a + v^ (0, 0) 

■ 5 
(« + i. '") 0) 

(a-\-~oj,0} 



n 



3 
5 

(«+äfc-,l) 



III 



(« + §^oj, 2) 
/ I 10 o. 



■ 6 „ 
(a + ä^t^iS) 

9 
(a + äÄ'"'^^ 



in welcher /c die Anzahl der Teile bedeutet, in welche das ursprüng- 
liche Intervall geteilt werden soll. Danach erhält man zunächst: 

(a + ~ CO, 1) -=(a + l <o, 0) — (a, 0), 



(74) 
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= («+ |.o,0) — 2(a+^«>,0) + (a,0) 

=^ ^"^ + ä '"' ^^ "" ^ ^"^ "l" ä '"' ^^ + ^^ + T '"' *^ 
^(a + |w,0) — 3(a+|«>,0) + 3(a + iw,0)-(a,O), 
oder allgemein: 

(« + 1^ '0, 9) = (« + 1 '"> 0) - (3), (a I- ^ w, 0) -t {qh (« 1" ^ "^ «) 

- T'(— D' (q)i (a + ^^ o,, 0). 

Hierin ist nun nach Formel (10): 
(«)o(a+|w,0)-(g)«(a,0) + (3)o(|)j (a f >,l) + («)«(|),(a+|«,2) |--- 

|-(3)o(|)^,,(a + '^'<',«-l) f (9)o(|\ (a4-|'",3)H ■■• 

+ <«)• (|)„ ^"^ -^-"^ '"' "^' 
(3)i(«+'T<'>,0) = (g)i (a,0)f(g)i (^)j(a f }«>,l)4-(2)i(^)^(a-f |<«,2) [- 

+ (3)i ( V"),-! ^" + '^ '"' *~^^ + ^*^' (^),/ ^" + ?' '"' ^^ ^' • • ■ 

{3)»(a+^^",0) = (2)2 (a,0)+(g)a (^)^(a+|a>,l) \-(q),(^)^(a -\fco,2) f. 

+ (3)» (^),.i (a -h T '". ^1) ^ (3)ä (^), (« + lOJ,q)-\-.. 
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(?Vi {a+jw, 0) = («Vi (a, 0)+(gV,(^)j (a+T«*' D + (9Vi(^)2(«+|«', 2) -|- - 

+ (9Vi (^),.i (« + ^ '0, 3-1) + (gVi (^\ (a + 1 <-, ?)+■•■ 
+ («Vi (r)„ (a f 1 <". «)» 



(«), («, 0) = (g), (a, 0). 
Man erhält somit durch Addition: 



2k 






(« + ^'o, 9) = > .(- D' (9). (« + T '". 0) 




(a, 0) V(- 1)' (g), + (a + 1 w, D X^" 1^' (3).- (^), 



^ 



+ (a + 1 «, 2) Jj- D' (g), (^X 
+ 



7^0 




+ (a + ^<«,g-l) > ^(-iv(g),(0)^^ 



1 = 



+ 



3)^(-l)'(3),(V), 



+ (a + ^ «>, «) y (- 1)' (?), (^)^, 
oder, wenn allgemein 




(-1)'(3).(V)=5 






(75) 



gesetzt wird, 

(«+|^«>,?)=4'^a,0) + 4'^a + |<",l) + 5^2'^ (« + !«', 2) f • 
(76)1 +5(?)(a + ?^a>,g_l)4-£(')(a + |«>,3)+... 

+ 4'^a + |o>,M). 
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Nnn ist nach (75): 

4" = «.. (I). - <* ( ti). +iM^)--± (,v, (1) 



k'r 



= (</)o 



— (?) 






1 • 



+ ((/) 



¥(^-') (^-^) (2=?'-3)''-(^-i>+l) 



2 • 



1 



P 



(?Vi • 



(i-l)(i-2)(i-3) (i-.+l) 



k\k 



+ , . Ä ^Ä *'' ^A "'' vä 



(9)o 



(i)'-^;:(ir+^r;(r- 



k 

±4"(i) 



p! 



f 



— ig) 



{m~^:m'+ö:;(^r- ...±ö-'> (^) 



1 • 



+ (q) 



(¥)'- ^H c-rrv ^: (¥)'- ■ ■ ic-" (2^) 



2 • 



i?! 



a)'-ö:;'(r+öinir--±qir9 



_ (9)g-l • 






H vft/ * 



oder nach Potenzen von k geordnet : 



(77) 



B 



(y) 



iL 

kP 






1)' iq)i (q—iy 



Op-l 

kP-^ 




(— 1)' {q)i (q—ir' 



i ---^ 



* Der Strich ( — ) über dem C soll andeuten, dass die mit den Zahlen 1 bis p — 1 
gebildeten Kombinationen ohne Wiederholung als Produkte aufzufassen sind. 
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(77) 



+ fe2j^-l)'(9M9-0''-- 



— k 




(— D' (q)i iq—iT^ 




(— D' (0, {q—t\ 



l-:{) 



worin : 



Vi 



77(1) 



(2) 



1 



(78)* 



bp.2 = 



. _p± 

(3) 



Ö/j-3 -^ 



c 



p-1 



PI 



f 




^ 



. (i)-i-l)(2p~i-l) 



. (p-'i-l)2 {2p-i-l) {p-i-2) 



C 



(p-i) 



&i = 



p-i 



p! 



1 • 2 - 3 - 4 ••• (i>-^l) 
p\ 



1 
P 



Berücksichtigt man noch, dass nach anderwärts angestellten Unter* 
suchungen die auf der rechten Seite von (77) stehenden Summen 
immer dann verschwinden, wenn p > g ist, so geht (76) über in: 



* Die Auffindung des den Zähler irgend eines der Koeffizienten b bildenden 
Ausdrucks mag an einem Beispiel erläutert werden: Es ist 



'-^ 1.3 
1-4 
1-5 



l'p—1 



f 2-3 
2.4 
2.5 

m m • • 

2'p—l 



• f 3.4 
3.5 



3-7^—1 



f ... \.(p-3)(p-2) -\-(p^2){p-l) 
(p-3) (p-l) 



^^ >-2)(P+l) ).g >^)(H-2) ^^_3 /i>--4)(H-3) ^ . . ._^ (^_3) . L(2^=?) 

b u Ct u 

\~ iP — 2) ^— ^- ~, mit dem allgemeinen Gliede : 

. (£_i _2H£+J0 oder ,: . (p-i--^)i2p-i-l\ 
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<« +|l«', ?) = ^^{a + \w, q) + 5$!^,(a+?±i o>, g + l) 



2k 



q+2 



(?) 



+ B\%2 (o + T "'. «' + 2) H h B^'^a + Io*,«), 



<79) 

mit welcher Beziehung die allgemeiae Lösung der eingangs dieses 
Paragraphen gestellten Aufgabe als abgeschlossen betrachtet 
werden kann, 

§ 52. Zur Erläuterung der vorstehend gefundenen Besultate 
machen wir die Annahme, die der gegebenen Tabelle zu Grunde 
gelegte Funktion sei eine solche vom fünften Grade. Damit 
«rhält man aus (77) unter Beachtung von (78) für 9 = 1 und 
p = l, 2, 3, 4, 5: 



<80) 



B^j» 



1 

k 



(1) 



B 






B 



(X) 

5 



2Ää 
1 

6Ä;ä 
1 
24Ä;* 
1 



2k 
1 



1 



2Ä2 ' SÄ 



4ä;8 



■ 11 1 
"•" 24*« 4Ä 



1 



i + 



120Ä;« 12Ä;* ' 24^» 12A2 ' bV 



+ 



für ? = 2 und p = 2, 8, 4, 5: 



<81) 



B 



für ^ = 8 und p 



B 



(2) 



B 



(2) 



k^ 



k^ 



B 



(2) 



3 , 11 



(2) 



12** 2*3 ' 12Ää 

1 Z__i_ J 

Rift ' ±1(3 



<82) 



ß 



(3) 



5 



(3) 



4*0 e/c' 

^ 3, 4, 5: 
3 3 



4*3 6*!»' 



B 



(3) 



2** 
5 



2A« 



4*0 ** ' 4*3' 



\+ 
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für ^7 = 4 und ^ = 4, 5 : 



(83) 



B 



(4) 



1 



5 Ä5 A* 



und endlich für ^ = 5, p = 5 : 



(84) 



ü(5)_ A 
X>- TT« 



Der Umstand, dass sich die Formel (79) für A = 1 auf 

reduziert, hat offenbar zur Folge, dass für diesen speziellen Wert 
von k alle B ausser B^i\ Bf\ Bf\ :^\ Bf^ (welche sämtUch = 1 
werden) verschwinden. Diese Tatsache liefert uns eine höchst 
einfache Kontrolle für die Richtigkeit der in (80) bis (84) an- 
gesetzten Formeln. So wird z. B. 

1 7.7 5 3—14421 — 10 



B 



<2) 



o 



7 . 7 5 3— 144-21 

4 6 '4 6"" 12 



= 0. 



Nachstehend geben wir jetzt noch die aus (79) abzuleitenden 
Beziehungen für die häufig notwendigen Fünf- und Zehnteilungm. 
Sie lauten für A = 5 : 



<86) 



21 



(a + 4 «'' *) + 15625 ^^'^ + ^ **' ^)' 
(a+^«;,2) = ^(aH-|w,2) — -^(aH-|w,3) 



625 
25 



125 



16 . . 4 .- 66 



+ 625 ^'*+-^"^'*^- 



3125 



(a 4- f w, 5), 



s 



(«+5i«>.3) 



1 - . 3 o. 6 ... , 4 



125 



(a + |,^,3)- — (a + 4.«,4) 



625 



6 . . 5 



+ ßö^(« + l""^'ö)' 



(''+^<'^>*) 



625 

6 , , 4 .- 8 



e 



625 






3125 



(a + 1 (0, 5), 



(a +^w,6) = ^— ^ (a + y w, 5), 



3125 
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und für A = 10: 

1 9 

I 57 ^ , 3 .-,. 1653 , , 4 ,^ 

, 64467 , .5 ., 
+ 4ÖÖÖÖÖÖ^"+^'"''^^' 



(« + a '0, 2) = Tnn (« + 1 '0. 2) - ^r^rpr^ (a + 1 «, 3) 



(86) 



309 .4 2679 

40000 ^" "•" 2 '"' ^ 400000 



, 0\jy , 1 4 -s iiO<W , I 5 ex 
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+ 80000^"+^^^'^^' 

(a -I- |ß>, 5) =. — — — (a + I oj, 5). 



Als numerisches Beispiel diene die in § 50 gegebene Tabelle 
XLIX. Für dieselbe wird: 

(«+4-,l) = Iö-7-2Öo -12 + ^.6^ 0,831 



(« +|->,2) = i7v^ • 12 - :,^ . 6 :^- 0,066 



i- 12--^ 
100 ' 1000 

(«+1^,3)==^. 6-^0,006, 

welche Werte mit den in der Tabelle L gefundenen in der Tat 
übereinstimmen. 

Anhangsweise mag darauf hingewiesen werden, dass die^ 
soeben gefundenen Zahlen selbstverständlich auch mit Hilfe der 
Formeln (67) direkt aus F (x) berechnet werden können und zwar 
in der Weise, dass in denselben co durch ^ oj ersetzt wird. 
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